Gonzalo Neira
Eva Vergara
David S. Holmes- 1Center for Bioinformatics and Genome Biology, Centro Ciencia & Vida, Fundación Ciencia & Vida, Santiago, Chile
- 2Facultad de Medicina y Ciencia, Universidad San Sebastián, Santiago, Chile
Verrucomicrobia are a group of microorganisms that have been proposed to be deeply rooted in the Tree of Life. Some are methanotrophs that oxidize the potent greenhouse gas methane and are thus important in decreasing atmospheric concentrations of the gas, potentially ameliorating climate change. They are widespread in various environments including soil and fresh or marine waters. Recently, a clade of extremely acidophilic Verrucomicrobia, flourishing at pH < 3, were described from high-temperature geothermal ecosystems. This novel group could be of interest for studies about the emergence of life on Earth and to astrobiologists as homologs for possible extraterrestrial life. In this paper, we describe predicted mechanisms for survival of this clade at low pH and suggest its possible evolutionary trajectory from an inferred neutrophilic ancestor. Extreme acidophiles are defined as organisms that thrive in extremely low pH environments (≤ pH 3). Many are polyextremophiles facing high temperatures and high salt as well as low pH. They are important to study for both providing fundamental insights into biological mechanisms of survival and evolution in such extreme environments and for understanding their roles in biotechnological applications such as industrial mineral recovery (bioleaching) and mitigation of acid mine drainage. They are also, potentially, a rich source of novel genes and pathways for the genetic engineering of microbial strains. Acidophiles of the Verrucomicrobia phylum are unique as they are the only known aerobic methanotrophs that can grow optimally under acidic (pH 2–3) and moderately thermophilic conditions (50–60°C). Three moderately thermophilic genera, namely Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus, have been described in geothermal environments. Most of the investigations of these organisms have focused on their methane oxidizing capabilities (methanotrophy) and use of lanthanides as a protein cofactor, with no extensive study that sheds light on the mechanisms that they use to flourish at extremely low pH. In this paper, we extend the phylogenetic description of this group of acidophiles using whole genome information and we identify several mechanisms, potentially involved in acid resistance, including “first line of defense” mechanisms that impede the entry of protons into the cell. These include the presence of membrane-associated hopanoids, multiple copies of the outer membrane protein (Slp), and inner membrane potassium channels (kup, kdp) that generate a reversed membrane potential repelling the intrusion of protons. Acidophilic Verrucomicrobia also display a wide array of proteins potentially involved in the “second line of defense” where protons that evaded the first line of defense and entered the cell are expelled or neutralized, such as the glutamate decarboxylation (gadAB) and phosphate-uptake systems. An exclusive N-type ATPase F0-F1 was identified only in acidophiles of Verrucomicrobia and is predicted to be a specific adaptation in these organisms. Phylogenetic analyses suggest that many predicted mechanisms are evolutionarily conserved and most likely entered the acidophilic lineage of Verrucomicrobia by vertical descent from a common ancestor. However, it is likely that some defense mechanisms such as gadA and kup entered the acidophilic Verrucomicrobia lineage by horizontal gene transfer.
Introduction
Acidophiles can be defined as organisms that grow optimally at pH <5 (Dopson, 2016). A common further classification defines moderate acidophiles as organisms that grow optimally between pH 3 and 5 (Foster, 2004; Benison et al., 2021) and extreme acidophiles that grow optimally at pH ≤ 3 (Johnson and Aguilera, 2015). The latter are particularly challenged for survival as they tolerate proton gradients across their membrane several orders of magnitudes higher than neutrophiles or moderate acidophiles (Baker-Austin and Dopson, 2007; Slonczewski et al., 2009; Quatrini and Johnson, 2018; Hu et al., 2020). Acidophiles are found in anthropogenic environments such as industrial mineral recovery (bioleaching) sites and acid mine drainage (AMD) and in natural environments including acidic hot springs, volcanic mud pools, and acid rock drainage (ARD) (Baker-Austin and Dopson, 2007; Huang et al., 2013; Cárdenas et al., 2016; Schmitz et al., 2021). Acidophilic organisms are important in the biotechnology industry as they play key roles in biomining processes and bioremediation of AMD (Johnson and Hallberg, 2005; Schippers et al., 2013). They are also a source of acid-stable enzymes (Baker-Austin and Dopson, 2007; Sharma et al., 2012; Dopson and Holmes, 2014; Liu et al., 2015). Furthermore, they are important from an evolutionary standpoint as some members inhabit what are thought to be early Earth-like environments (González-Toril et al., 2005; Baker-Austin and Dopson, 2007; Mirete et al., 2017). Extreme acidophiles are found across all the Tree of Life (ToL) with members in Bacteria, Archaea, and Eukarya domains (Hedrich and Schippers, 2021).
The phylum Verrucomicrobia is composed mainly of neutrophiles found commonly in soil and fresh or marine waters (Bergmann et al., 2011; Freitas et al., 2012; Cabello-Yeves et al., 2017) and also as part of the intestinal microbiome (de Vos, 2017). Exceptionally, a clade of extremely acidophilic aerobic methanotrophs has been described in geothermal environments (Dunfield et al., 2007; Pol et al., 2007; Islam et al., 2008; van Teeseling et al., 2014). This clade includes organisms from three different genera: Methylacidiphilum, identified in hot acidic geothermal environments (optimal growth pH 2 - 3 and optimal temperature growth 50–60 °C) (Hou et al., 2008; Islam et al., 2008; Anvar et al., 2014); Methylacidimicrobium, from the same geothermal location but living at lower temperatures (optimal growth temperature 30 to 44°C and optimal growth pH between 1 and 3) (van Teeseling et al., 2014) with the exception of Methylacidimicrobium thermophilum AP8 which grows optimally at 50 °C and pH 3 to 5 (Picone et al., 2021a). A third organism, Ca. Methylacidithermus pantelleriae, from a metagenome-assembled genome (MAG) represents the only known member of this proposed genus (Picone et al., 2021b). All three genera are characterized by being the only bacteria capable of oxidizing methane (methanotrophy) in extremely acidic environments (van Teeseling et al., 2014). Recently, it has also been discovered that organisms from these genera can grow as mixotrophs by oxidation of H2, in addition to methane, as an adaptation to the fluctuations of methane in the geothermal environments they inhabit (Carere et al., 2017; Mohammadi et al., 2017, 2019; Schmitz et al., 2021). They also use lanthanides instead of other metals as co-factors in some enzymes (Pol et al., 2014) and are being exploited for the development of industrial applications of rare earth elements including lanthanide (Cotruvo, 2019; Featherston and Cotruvo, 2021).
Phylogenetic profiles of Verrucomicrobia acidophiles have been based on analyses using 16S rRNA (Dunfield et al., 2007; Pol et al., 2007; van Teeseling et al., 2014; Picone et al., 2021b; Schmitz et al., 2021), pmoA (Anvar et al., 2014; Picone et al., 2021a), or seven housekeeping genes (Erikstad et al., 2019). In this paper, we extend the phylogenetic analyses of this group using fully sequenced genomic data, including all sequenced genomes of Verrucomicrobia methanotrophs.
To date, most of the extensive studies of extreme acid resistance and its evolution have been made in organisms that aerobically oxidize reduced iron and/or sulfur in biomining operations and similar acidic econiches (e.g., AMD) (Baker and Banfield, 2003; Valdés et al., 2008; Cárdenas et al., 2016; Christel et al., 2018; Vergara et al., 2020; Zhang et al., 2020b; González-Rosales et al., 2022). These organisms maintain a near neutral pH inside the cell while confronting an external pH that can be several orders of magnitude lower (e.g., pH ≤ 3) (Baker-Austin and Dopson, 2007). Acid resistance mechanisms have been classified for convenience of discussion into two different classes, a “first line” of defense, which encompasses membrane changes designed to impede the entry of protons into the cell, and a “second line” of defense, which either consumes or expels protons that have penetrated the cytoplasmic membrane (Vergara et al., 2020).
Previous studies have identified genes predicted to form part of the “first line” of defense in the moderate thermophilic group Methylacidiphilum (Hou et al., 2008; Kruse et al., 2019). Six different copies of the starvation lipoprotein (Slp), a protein with no proven function but that is co-expressed in an acid resistance island in E. coli (Mates et al., 2007) were identified in M. infernorum V4 and five different orthologs also identified in both other complete genomes M. fumariolicum SoIV and M. kamchatkense Kam1. Also, two potassium transporters, Kch and Kdp, were described which have been shown in other microorganisms to generate an inner membrane positive charge that repels protons (Buetti-Dinh et al., 2016; Christel et al., 2018). In addition, modification of the fluidity of the membrane via the incorporation of hopanoids has also been proposed to be present in a wide range of acidophiles and deletion of the hopanoid biosynthesis pathway in acid-tolerant bacteria reduces the growth at low pH (Welander et al., 2009; Mangold et al., 2013; Sohlenkamp and Geiger, 2016), but this has not yet been investigated in the acidophilic Verrucomicrobia.
Mechanisms, potentially forming part of the second line of defense against acid resistance, have been reported previously in some strains of Methylacidiphilum but are unexplored in Methylacidimicrobium and Methylacidithermus (Hou et al., 2008; Kruse et al., 2019). Amino acid decarboxylation systems, such as arginine and glutamate decarboxylation, include enzymes that work in pairs with the corresponding transporter by binding protons and then secreting them outside the membrane (Foster, 2004). Another second line of defense mechanism involves proton/cation antiporters (e.g., sodium or calcium cations) that belong to a large family of transporters that can regulate the intracellular pH by the expulsion of protons while maintaining the cytoplasmic charge unchanged by the entrance of a positively charged cation (Padan et al., 2004). The ATPase F0-F1 is a common mechanism that can work either by importing protons to generate ATP or in a reverse fashion resulting in the expulsion of the protons (Foster, 2004).
In this study, we expand our knowledge of potential acid resistance systems in strains of the extremely acidophilic genera, Methylacidiphilum, Methylacidimicrobium, and Ca. Methylacidithermus. Predicted acid resistance genes and pathways were then used to build a global model of pH homeostasis mechanisms in these genera that allowed us to hypothesize how acid resistance was acquired during the evolution of acidophilic methanotrophs. It is hypothesized that vertical inheritance from a common ancestor that is predicted to be a neutrophile and Horizontal Gene Transfer (HGT) from other acidophiles are the main mechanisms for the acquisition of acid resistance in acidophilic methanotrophs.
Materials and methods
Retrieval of genome sequences and quality assessment
A dataset of 20 genome sequences from acidophilic methanotrophs of the Verrucomicrobia phylum were obtained from AciDB (Neira et al., 2020) and retrieved using the National Center for Biotechnology Information (NCBI) assembly ID through the Batch Entrez tool (Kans, 2020). In addition, 49 genomes of neutrophiles classified as representative species in the genome taxonomy database (GTDB) (Parks et al., 2022) and NCBI were retrieved using the Batch Entrez tool. From these genomes, 24 are associated with the phylum Planctomycetes and were used exclusively for phylogenetic analysis and 25 from the Verrucomicrobia phylum were used for both phylogenetic and comparative genomics analysis. Quality assessment of all genomes was carried out by CheckM v1.1.2 (Parks et al., 2015), defining >90% completeness and <10% contamination as high-quality genomes according to the proposed standards for high-quality draft genomes (Chun et al., 2018). Detailed information on the nomenclature, taxonomy, and source of the genomes in our study is provided in Supplementary Table 1.
Genome characterization
Average nucleotide identity (ANI) between genomes was calculated with FastANI v1.1 using default parameters, where organisms with ANI values >95% were considered to belong to the same species (Jain et al., 2018). In a group of species clustered together with ANI > 99%, only the highest quality genome was selected for subsequent analysis. Digital DNA-DNA hybridization (dDDH) was carried out using the genome-to-genome distance calculator v2.1 (GGDC) (Meier-Kolthoff et al., 2013) with Formula 2 as recommended for datasets with draft genomes. Average amino acid identity (AAI) was calculated between each pair of genomes using CompareM (https://github.com/dparks1134/CompareM), using default parameters.
Phylogenomic analyses
Phylogenetic reconstruction was inferred using the PhyloPhlAn3 pipeline using a set of 400 universal markers (Asnicar et al., 2020). The parameters utilized in each of the steps of the PhyloPhlAn3 pipeline are the following: Diamond v2.0.11 (Buchfink et al., 2015) was used for mapping each marker in the different genomes; MAFFT (Katoh et al., 2005) with the L-INS-I algorithm for the multiple sequence alignment (MSA) of identified markers; IQTREE v1.6.12 (Nguyen et al., 2015) was used for the construction of maximum-likelihood trees from the concatenated MSA; Modelfinder (Kalyaanamoorthy et al., 2017) was used to select the best-suited evolutionary models according to the Bayesian (BIC) and Akaike information criterion (AIC); and UFBoot (Minh et al., 2013) was used to infer the robustness of the tree with 1,000 replicates. Domain composition of proteins was predicted using used PFAM domains in TREND (Gumerov and Zhulin, 2020, 2022). The consensus tree was visualized using iTOL v6 (Letunic and Bork, 2019) and edited with Inkscape (https://inkscape.org/).
Identification of orthologous groups of genes
Orthogroups are defined as proteins that have evolved from a common ancestor and include both orthologs (homologous proteins from two species) and paralogs (set of proteins that have a common origin in the same genome) and were identified using Orthofinder v2.3.3 (Emms and Kelly, 2018). Predicted protein sequences were compared using Diamond v2.0.11 (Buchfink et al., 2015) in an all-versus-all sequence similarity search. DendroBLAST (Kelly and Maini, 2013) was used to generate unrooted protein trees for each of the identified orthogroups. Protein alignment for all orthogroups was performed using MAFFT (Katoh et al., 2005). FastTree (Price et al., 2010) was used for tree inference in each orthogroup.
Acquisition and analyses of acid resistance genes
Potential genes and pathways involved in microbial resistance against low pH were identified through an extensive literature search and were functionally characterized using Pfam (Finn et al., 2014) and InterPro families (Blum et al., 2021), providing a curated list of proteins associated with acid resistance (Supplementary Table 2). All genomes were functionally annotated using the Interproscan pipeline (Jones et al., 2014). InterPro and Pfam were then used to identify similar proteins in the genomes of acidophiles and neutrophiles of Verrucomicrobia using a custom python script. All proteins predicted to be involved in acid resistance were identified in the previously defined orthogroups in order to analyze the conservation across the organisms in the study (Supplementary Table 3). The conservation of genetic context between Verrucomicrobia genomes and manually curated mis-annotations for the genes of interest were determined by MAUVE (Darling et al., 2004). Phylogenetic tree reconstructions were made as described above for all orthogroups involved in acid resistance.
Evolutionary inference of acid resistance genes
Inference of branch-site-specific events was made using the tree topology of the 400 conserved markers species tree obtained in the phylogenomic analysis as described above. The presence and absence of genes related to acid resistance and associated genes (same genomic context) were mapped onto each branch of the phylogenetic tree. The inference of evolutionary events was made using maximum parsimony criteria as incorporated in the COUNT software (Csuös, 2010). HGT predictions were made using HGTector (Zhu et al., 2014). To study in more detail, possible HGT events phylogenetic trees were constructed for each predicted protein to be acquired by these mechanisms, and in these cases, the 100 best hits from BLASTP against the RefSeq protein database followed by removal of duplicate sequences using SeqKit (Shen et al., 2016) were added to the sequences in the corresponding orthogroup. Phylogenetic tree construction was performed with the same parameters and software detailed in the previous section with the additional step for the estimation of the root position for each orthologous unrooted tree with minimal ancestor deviation (MAD) method applied with MAD phylogenetic rooting algorithm (Tria et al., 2017).
Results and discussion
Genomic index and features
Twenty genome sequences from acidophiles of the bacterial phylum Verrucomicrobia were downloaded, of which 7 are complete genomes and 13 are in draft state. All genomes have a contamination index lower than 10% and a completeness index higher than 90% as defined by CheckM. Acidophiles and neutrophiles genome sizes showed differences; the acidophilic clade has an average genome size of 2.37 Mb, in contrast to neutrophiles clades of Opitutales, Chtoniobacterales, and Verrucomicrobiales having an average size of 5.34, 6.3, and 4.28 Mb, respectively (Cortez et al., 2022) (Supplementary Table 1). Genomic indices were calculated to identify the taxonomic rank and species for the genomes under study. ANI and dDDH values were used to evaluate the species clusters in the dataset; in addition, ANI values over 99% were used to de-replicate the genomes for the subsequent analysis (values summarized in Supplementary Tables 4, 5, respectively). Species clusters not previously reported were identified in the case of Methylacidiphilum sp. IT5 with Methylacidiphilum infernorum (ANI of 96% and dDDH of 68.7 using formula 2 with a 75% probability of both genomes being the same species) and Methylacidiphilum sp. IT6 with Methylacidiphilum sp. Phi (ANI of 98% and dDDH of 86). All previously reported species were also supported by this analysis. Five genomes from Methylacidiphilum fumariolicum and one genome of Methylacidiphilum kamchatkense were excluded for further analysis by de-replication criteria, and in both cases, the complete genome of each species was kept. AAI values were used to analyze the genus-level taxonomic association of the genomes (Supplementary Table 6). The values divide the different species into three different genera: (i) Ca. Methylacidithermus pantelleriae is the only representative for Methylacidithermus; (ii) Methylacidimicrobium with six different species, and (iii) Methylacidiphilum with five different species. An overview of the final 14 genomes with their genomic characteristics and geographic origin is presented in Table 1.
Table 1. Genomic features of the organisms used in this study.
Phylogenomic analysis
A phylogenomic tree was reconstructed based on the 400 conserved markers identified with PhyloPhlAn3 resulting in an alignment of 9,640 positions (Figure 1). The tree was rooted using genomes from the Planctomycetes phylum as an outgroup. Two main clades are distinguished, the first one formed by acidophilic genomes from all three genera and Limisphaera ngatamarikiensis (only representative genome of the Pedosphaerales order) and the second one with all neutrophiles from the orders Opitutales, Chthonobacteriales, and Verrucomicrobiales. The branching order in the acidophilic clade suggests that Ca. Methylacidithermus pantelleriae is the first to branch followed by Methylacidimicrobium and Methylacidiphilum. This genome-driven phylogeny supports the branching order reported for a 16S rRNA tree (Schmitz et al., 2021) and consolidates and expands previous phylogenomic studies that used gene markers (Cai et al., 2022; Kim et al., 2022).
Figure 1. Phylogenetic distribution of Verrucomicrobia. Maximum Likelihood tree obtained with PhyloPhlAn3 using 400 conserved markers. Planctomycetes genomes were collapsed and used as an outgroup. Support values of the ultrafast bootstrap repetitions are shown in each corresponding branch. The taxonomic rank of each clade is shown at the right of the tree, genus level for all the acidophilic genomes, and order level for the neutrophiles (in yellow shadow).
Acid resistance: First line of defense
Proton influx into the cell is thought to be reduced in acidophiles by the development of an inside positive potential which is different from neutrophiles that have inside negative potential (Richard and Foster, 2004; Baker-Austin and Dopson, 2007). One of the ways that acidophiles generate this positive potential is via the use of membrane cation transporters, especially K+ transporters. Acidophiles also have a highly impermeable cell membrane to restrict proton influx into the cytoplasm (Konings et al., 2002) that can be generated by altering the lipid and fatty acid composition of the membrane (Feng et al., 2021). Hopanoid incorporation into the membrane has also been implicated in enhancing membrane stability (Welander et al., 2009).
Potassium uptake systems
Three different potassium uptake systems were identified in the methanotrophs:
(i) a voltage-associated potassium channel (vCh) was predicted in all acidophiles and in 8 of 25 neutrophiles. A phylogenetic tree was constructed with vCh and rooted with MAD (as described in Methods). Three main clades of vCh (clades I, II, and III, Figure 2) were identified, all sharing a common architecture with similarity to the potassium transporter KcsA; the latter could be a potential starting point in the evolution of potassium channels (Diskowski et al., 2015). Domain composition analysis (Zhang et al., 2020a) of clade I that contains sequences from Verrucomicrobiales and Opitutales suggests that it evolved as a fusion between TrkA and KcsA. Clades II and III share a similar architecture with an N-terminal KcsA-like domain and a C-terminal KtrA-like domain (Figure 2). Clade II contains Verrucomicrobiales that branch close to Planctomycetes, Cyanobacteria, and Alphaproteobacteria. Clade III contains all the acidophilic methanotrophs in a monophyletic branch that branches as the final “twig” from a group that contains Deltaproteobacteria, Cyanobacteria, and Chlorobiaceae. The diversity in phylum and families forming the clade III, out from Verrucomicrobia, strength the HGT possibility for this gene in the acidophiles. The last is supported by additional positive results for the Methylacidiphilum from HGTector that suggests this gene could be related to an HGT event (Supplementary Table 7).
(ii) A second potassium uptake system identified is the high-affinity system Kdp (Epstein, 2003). The classic composition of the Kdp cluster is characterized by the presence of the effector subunits KdpFABC, in charge of forming the inner membrane potassium transporter and the two-component regulatory system KdpDE. The gene cluster formed by kdpABC was identified in Methylacidimicrobium and Methylacidiphilum genera, with kdpF detected exclusively in these acidophiles. The subunit kdpF is known to be involved in the stabilization of the potassium transporter (Dorus et al., 2001), suggesting that its presence in the acidophiles may be an adaptation to low pH. In addition, the genomic context of the regulatory genes kdpDE in the acidophilic genera is different from the classic organization detected in neutrophiles (Figure 3) (Epstein, 2003). Prediction of HGT events coupled with parsimony analysis shows that the kdpDE cluster in these genera is most likely inherited vertically (Figure 8). Interestingly, genome rearrangement events have been linked to adaptation in different ecological niches (Darling et al., 2008; Qian and Zhang, 2014; Yan et al., 2018) and we hypothesize that the potential translocation of kdpD/kdpE in the acidophiles plays a role in low pH adaptation possibly by modifying the expression of the effector proteins KdpFABC in acidic conditions (Ballouz et al., 2010; Cholo et al., 2015). This remains to be tested experimentally.
(iii) We detected a third potassium uptake system (Kup) in acidophiles of Methylacidimicrobium and also in several related neutrophiles. However, phylogenetic tree reconstruction shows that Kup of the acidophilic Verrucomicrobia clusters close to Kup of other extreme acidophiles such as Acidithiobacillus and Acidiferrobacter rather than Kup of neutrophilic relatives (Figure 4), strongly suggesting that it was acquired by HGT, which is supported by HGTector results predicting Acidithiobacillus as a putative donor (Supplementary Table 7). As Kup is embedded in the membrane directly facing the physicochemical constraints of the low pH environment, this HGT event could enable the acidophilic methanotrophs to respond to the challenge of low pH. This is consistent with the suggestion that acidophiles are prone to share genes that carry adaptation to their specific econiche (Beard et al., 2021), and supports the proposed role of Kup in some acidic environments (Trchounian and Kobayashi, 1999).
Figure 2. Phylogenetic tree and domain composition of potassium transporter vCh. In the blue background is shown clade I, yellow = clade II, and purple = clade III. All sequences were clustered to phylum or class level when possible. Acidophiles are indicated in red letters and the Verrucomicrobia neutrophiles are in bold letters. On the right side is shown the domain composition of proteins in each clade as predicted by PFAM in TREND, where ion_trans_2 = domain responsible for recognition of K+. At the right bottom of a box are the Model Domains for the different potassium channels (Uniprot accession as reference); TrkA from Escherichia coli (acc: P0AGI8) and Salmonella typhimurium (acc: P0A2J9), Ktr from Bacillus subtilis (acc: O32080), KcsA from Streptomyces lividans (acc: P0A334), and Kch from E. coli (acc: P31069). Baln, Balneolaceae; Planc, Planctomycetes; Chlor, Chloroflexi; Bact, Bacteroidetes; Abd, Abditibacteriota; Deltapr, Deltaproteobacteria; Thdesbact, Thermodesulfobacteriales; Thdesvi, Thermodesulfovibrio; Spir, Spirochetes; Firm, Firmicutes.
Figure 3. Cluster organization for kdp genes. Collapsed cladogram of Verrucomicrobia is represented on the left. Functional genes are shown in purple and regulator subunits are shown in orange, with synteny represented by yellow connecting background. *Insertion of hypothetical protein identified in Methylacidiphilum kamchatkense. **Transposase identified in Phragmitibacter flavus and Brevifollis gellanilyticus.
Figure 4. Rooted phylogenetic tree of potassium transporter Kup. Acidophiles clades are shown in red with violet background for acidophiles from Verrucomicrobia. Yellow background and bold letters for neutrophilic Verrucomicrobia, green background for neutrophiles from Proteobacteria, and other phyla are without color background.
Membrane modifications
The starvation lipoprotein (Slp) is one of the proteins identified involved in acid resistance (Alexander and St John, 1994). Slp was detected exclusively in the acidophilic and not the neutrophilic Verrucomicrobia. The number of copies ranged from a minimum of five in Methylacidiphilum fumariolicum to seven in the Methylacidimicrobium genus. A rooted phylogenetic tree with the best hits from NCBI for the Slp family was done (Figure 5), showing a close relation with sequences from Gammaproteobacteria, including sequences from the acidophiles Acidihalobacter. The sequences from Verrucomicrobia form different clades suggesting that there are possible different origins for the multiple Slp gained. The lipobox characteristic for this protein was present in 55 of 79 of the orthologs. However, the classic Asp was not found in the +2 position suggesting that Slp is targeted to the outer membrane (Zückert, 2014; Vergara et al., 2020). Extensive gene duplication of slp in extreme acidophiles has been noted previously in the genera Leptospirillum (Vergara et al., 2020) and Acidihalobacter (Boase et al., 2022). Currently, little is known about the mechanism of action of Slp. Transcription analysis in E. coli reveals that a gene cluster denominated acid fitness island (AFI) that includes slp which is overexpressed under acidic conditions (Tucker et al., 2002) and it is also known that slp mutants have a diminished capability to grow in extremely acidophilic pH (Mates et al., 2007), arguing for a link between the presence of slp and the capability to grow in low pH conditions.
Figure 5. Rooted phylogenetic tree for Slp in Verrucomicrobia acidophiles (violet background) and their best Blast hits. In red are the Slp predicted proteins from acidophiles and in red asterisks are marked the collapsed clades where an acidophile protein sequence is present. The clades were collapsed to class or phylum level when possible, indicating the number of leaves in parenthesis.
One of the main mechanisms identified in acidophilic bacteria to participate in membrane modifications is hopanoids biosynthesis (Belin et al., 2018; Vergara et al., 2020). The genes for squalene biosynthesis, which is the precursor of hopanoid, are present in all acidophiles and 15 of 25 neutrophiles, however, with different genes to complete the process. For instance, in neutrophiles, the squalene synthase (sps) gene was identified, while in acidophiles, hpnCDE were identified as the enzymes involved in the production of squalene, following a canonical synthesis pathway (Kannenberg and Poralla, 1999). The squalene production genes in several Verrucomicrobia genomes could be explained by their multiple roles as the precursor for sterols and pentacyclic triterpenes in addition to hopanoids, molecules that perform a range of functions on the stabilization of membranes not only, particularly, associated with low pH adaptation (Belin et al., 2018). In contrast, hpnF codifying for squalene hopane cyclase and hpnG producing the first branched hopanoid molecules were identified in all acidophiles but only in three neutrophiles (Chthoniobacter flavus Ellin428, Terrimicrobium sacchariphilum, and Limisphaera ngatamarikiensis). The hopanoid transporter hpnN required for the transport of the cytoplasm synthesized hopanoid molecules to the membranes was identified in the organisms that have hpnFG with the addition of two other Verrucomicrobiales genomes. The predicted final products in the pathway for acidophiles and Limisphaera ngatamarikiensis could be the branched hopene molecules bacteriohopanetetrol (BHT) and aminobacteriohopanetriol (ABH), in contrast to most neutrophiles in Verrucomicrobia, reaching until squalene (Figure 6). HpnO and HpnA were identified in the majority of organisms. The differences in the hopanoid biosynthesis pathway could be key in the development of acid resistance as the branching hopene molecules have been related to adaptation to extremely low pH providing the necessary membrane modifications to reduce the entrance of protons to the cytoplasmic space (Mangold et al., 2013; Belin et al., 2018).
Figure 6. Predicted hopanoid biosynthesis in Verrucomicrobia. The percentage of organisms with each of the proteins is presented at the side of each specific reaction arrow. The percentage in acidophiles is presented in red and neutrophiles in blue. Gray dotted arrow represents a not verified step mediated by HpnA, double arrow indicates spontaneous reactions. FPP, Farnesyl diphosphate; BHT, bacteriohopanetetrol; ABH, Aminobacteriohopanetriol.
Acid resistance: Second line of defense
Cation antiporters
The sodium–proton exchanger (NaH, nhe2) is a membrane protein that transports Na+ into the cell and H+ out of the cell (antiport). It has been implicated in pH homeostasis in a wide range of species (Masrati et al., 2018). NaH, encoded by nhe2, is predicted to occur in the acidophilic methanotrophs Methylacidiphilum and Methylacidimicrobium where it is present in multiple copies except for Methylacidiphilum sp. IT5 which has only one copy (Supplementary Table 3). Phylogenetic analysis indicates that the acidophilic methanotrophs have three distinct clades of NaH suggesting that it may have different origins/functions. The first clade (termed NaH-1) contains sequences from both of the acidophilic genera and is closely related to sequences from Firmicutes and Actinobacteria (Supplementary Figure 1). The second clade (NaH-2) is associated only with Methylacidiphilum, branches close to thermophilic microorganisms such as Thermus. There is a third clade (NaH-3) that branches close to Proteobacteria, and we speculate that it is not necessarily involved in acid resistance (Supplementary Figure 1). We hypothesize that the NaH restricted to the acidophiles helps in proton export in acidic environments. It remains to be deduced whether it arose by HGT from other acidophiles or by gene duplication from the neutrophilic NaH representative that was present in the ancestral line.
The calcium/proton antiporter (Cah, chaA) was identified in all species of the Methylacidiphilum genus and is predicted to have been gained by HGT. To study further the distribution of chaA, a phylogenetic tree analysis was performed. Two clearly distinct clades could be identified (Supplementary Figure 2). Clade one contains sequences from Methylacidiphilum which branch close to, but distinct from, Proteobacteria and Nitrospira lenta and also contains the sequences from Opitutales and one Chtoniobacterales. Clade two does not have any Methylacidiphilum representatives but contains sequences from Chtoniobacterales and Verrucomicrobiales. The distance between Cah from the acidophiles and other members of the Verrucomicrobia is consistent with the idea that it was acquired by HGT supported by HGTector results (Supplementary Table 7). Alternatively, chaA in the acidophiles has been subjected to evolutionary pressure to adapt to dealing with low pH environments.
Decarboxylation systems
A GadA (glutamate decarboxylase) decarboxylation system, potentially involved in acid resistance (Bearson et al., 1997; Foster, 2004), was identified in both acidophilic and neutrophilic Verrucomicrobia. GadA was identified in all acidophilic species and only in two neutrophiles (Supplementary Table 3). Phylogenetic tree analysis of GadA (Figure 7) indicates that its closest relatives include multiple examples from other acidophiles of the genus Acidithiobacillus, suggesting that it was obtained by HGT. GadA has been shown to enhance the ability of bacteria to cope with extremely low pH (Damiano et al., 2015) and to function in E. coli when stressed at low pH (Foster, 2004). GadA works in conjunction with the antiporter system for the export of protonated products. Two conserved potential antiporter systems, GadC and AdiC (Hou et al., 2008; Kruse et al., 2019), were detected in methanotrophs but it remains to be demonstrated experimentally that they function in collaboration with the decarboxylation systems.
Figure 7. Rooted phylogenetic tree for glutamate decarboxylase GadA. In violet and yellow background are highlighted the acidophiles and neutrophiles from Verrucomicrobia, respectively, with bold letters for the ones present in this study. Collapsed clades to phylum or class level, when possible, with number of leaves between parentheses. Alphaprot, Alphaproteobacteria; Betaprot, Betaproteobacteria; Gammaprot, Gammaproteobacteria; Deltaprot, Deltaproteobacteria; Planct, Planctomycetes; Acctinob, Actinobacteria; Chlor, Chloroflexi; Cyanob, Cyanobacteria.
ATPase F0-F1
ATPase F0-F1 is a system used for both the generation of ATP and the extrusion of protons in highly acidic conditions (Sun, 2016). Two different ATPase F0-F1 clusters were identified in the Verrucomicrobia genomes. The first, conserved across all organisms exhibit a genetic structure that follows the canonical organization atpBEFHAGDC of an F-type ATPase involved in ATP synthesis. The second cluster is only present in Methylacidiphilum and Methylacidimicrobium and has the organization atpDCZBEFAG, resembling that of an N-type ATPase (Supplementary Figure 3A) (Schulz et al., 2017). The motif ESLxxY was detected in subunit atpE (Supplementary Figure 3B), which is highly conserved in acidophiles and has been suggested to be involved in proton extrusion as a mechanism for acid resistance (Foster, 2004; Sun et al., 2012).
Intracellular buffering
Intracellular buffering is another mechanism commonly used by cells for internal pH homeostasis (Slonczewski et al., 2009). One of the pathways predicted in both acidophilic and neutrophilic Verrucomicrobia is the arginine/agmatine deiminase pathway. This pathway includes the participation of the arginine decarboxylase (SpeA or AdiA) and two additional proteins: agmatine deiminase (AguA) and N-carbamoylputrescine amidohydrolase (AguB). The synteny of aguAB is highly conserved in Verrucomicrobia but their genomic context differs between acidophilic and neutrophilic representatives. For example, all genomes from the Methylacidimicrobium clade have a syntenic downstream region that contains the hopanoid transporter hpnN and the Methylacidimicrobium and Methylacidiphilum genera have an upstream slp (Supplementary Figure 4). Since we have detected a possible relationship between slp and hpn with acid resistance, it is possible that aguAB are also associated with pH homeostasis.
Another mechanism commonly used for intracellular buffering is the phosphate uptake system (pstABCS) (Hutkins and Nannen, 1993; Peng et al., 2017). PstABC was identified in conserved orthogroups across acidophilic and neutrophilic Verrucomicrobia and we suggest that this gene cluster was acquired by vertical inheritance. In contrast, the gene potentially encoding the phosphate binding subunit PstS was identified only in three neutrophiles but was conserved across all acidophiles. To explore this in further detail, a phylogenetic tree was constructed. This tree shows a clear separation of pstS between the acidic and neutrophilic clades (Supplementary Figure 5) suggesting that its acquisition was by HGT.
Carbonic anhydrase (Can) has also been implicated in intracellular buffering by a proton-consuming reaction and is upregulated under acidic conditions in Helicobacter pylori (Bury-Moné et al., 2008). Ca was predicted in both acidophilic and neutrophilic Verrucomicrobia but could be separated by phylogenetic analysis into two orthogroups corresponding, respectively, to acidophilic and neutrophilic Verrucomicrobia, with the caveat that Ca. Methylacidithermus and two organisms from Methylacidimicrobium have copies of both orthogroups (Supplementary Table 3). Using bioinformatic approaches, we could not detect functional differences between the two orthogroups and their potential role in pH homeostasis in the Verrucomicrobia remains to be elucidated.
Evolutionary trajectory of acid resistance genes
To gain insight into the evolutionary history of the acidophilic Verrucomicrobia, genes predicted to be involved in acid resistance were mapped onto a cladogram derived from the maximum likelihood tree of Verrucomicrobia. Gene gains (HGT and gene duplications) were annotated using Dollo parsimony analysis and HGT prediction (Figure 8). The cladogram is consistent with the idea that the acidophiles (blue background, Figure 8) evolved from a common ancestor shared with the neutrophiles (yellow background, Figure 8). Interestingly, the early evolution of the acidophiles appears to be concurrent with the development of moderate thermophilia (growth ~50°C) (red branches, Figure 8). The inferred neutrophilic ancestor already had several genes involved in acid resistance including kdpABCD and various Hpn modifications (black genes, Figure 8). We speculate that these could have provided an advantage for growth in moderately low pH conditions (pH 4-6). However, for growth in extremely low pH, various genes were gained such as slp, N-atp, hpnCDE, and nhe (green highlighted gene, Figure 8). It is proposed that several of these genes were gained by HGT from other acidophilic lineages as discussed above; other genes such as slp are predicted to have been gained by HGT followed by rampant gene duplication.
Figure 8. Hypothesized evolutionary trajectories of genes predicted to be involved in acid resistance in Verrucomicrobia (cladogram). Acidophiles = blue background with purple numbered nodes; neutrophiles = yellow background with unnumbered orange nodes. Clades with red branches are moderate thermophiles and black branches are mesophiles. Genes in black are inferred to be inherited by vertical descent, genes in green are predicted to be inherited by 1HGT or 2gene duplications. Gene events were predicted by Count with Dollo parsimony, HGTector, and phylogenetic analyses.
Reduction of genome size was also observed in the proposed evolutionary trajectory with an average genome size of the neutrophiles of 5.3 Mb, whereas it was 2.37 Mb in the acidophilic clades (Supplementary Table 1). Genome reduction has been proposed to be an adaptive mechanism to reduce energy costs involved in replication and translation in nutrient-scarce environments (Sowell et al., 2009) and under acid stress (Cortez et al., 2022).
Model
A model of genes and mechanisms predicted to be involved in acid resistance in extreme acidophiles belonging to the Verrucomicrobia group is shown in Figure 9. The model incorporates the concept of first and second lines of defense, where the first line of defense consists of mechanisms that are used to impede the entry of protons into the cell. These are principally, but not exclusively, located in the cell membrane. Whereas second line of defense mechanisms are mainly involved in the neutralization or expulsion of protons that have either escaped the first line of defense and have entered the cell or have been produced within the cell.
Figure 9. Model of the proposed acid resistance mechanisms in Verrucomicrobia acidophiles. Top section shows the predicted first line of defense mechanisms: blue = potassium/H+ exchange systems involved in the generation of a reversed membrane potential and the membrane systems Slp and Hpn; green, Slp; yellow, proteins associated with hopanoid biosynthesis. Bottom section shows the proposed second line of defense mechanisms: purple, decarboxylation and antiporter systems; green, N-type ATP synthase; pink, Pst transporter and carbonic anhydrase.
Data availability statement
The original contributions presented in the study are included in the article/Supplementary material, further inquiries can be directed to the corresponding author/s.
Author contributions
GN and DH conceived and designed the research. GN performed the research. GN, EV, and DH analyzed the data and wrote the manuscript. All authors participated in the writing and approval of the final manuscript.
Funding
This work was supported by FONDECYT 1181717 (DH), Centro Ciencia & Vida, and FB210008 Financiamiento Basal para Centros Científicos y Tecnológicos de Excelencia de ANID.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2022.900531/full#supplementary-material
References
Alexander, D. M., and St John, A. C. (1994). Characterization of the carbon starvation-inducible and stationary phase-inducible gene slp encoding an outer membrane lipoprotein in Escherichia coli. Molec. Microbiol. 11, 1059–1071. doi: 10.1111/j.1365-2958.1994.tb00383.x
Anvar, S. Y., Frank, J., Pol, A., Schmitz, A., Kraaijeveld, K., den Dunnen, J. T., et al. (2014). The genomic landscape of the verrucomicrobial methanotroph Methylacidiphilum fumariolicum SolV. BMC Genom. 15, 1–12. doi: 10.1186/1471-2164-15-914
Asnicar, F., Thomas, A. M., Beghini, F., Mengoni, C., Manara, S., Manghi, P., et al. (2020). Precise phylogenetic analysis of microbial isolates and genomes from metagenomes using PhyloPhlAn 3.0. Nat. Commun. 11, 1–10. doi: 10.1038/s41467-020-16366-7
Awala, S. I., Gwak, J. H., Kim, Y. M., Kim, S. J., Strazzulli, A., Dunfield, P. F., et al. (2021). Verrucomicrobial methanotrophs grow on diverse C3 compounds and use a homolog of particulate methane monooxygenase to oxidize acetone. ISME J. 15, 3636–3647. doi: 10.1038/s41396-021-01037-2
Baker, B. J., and Banfield, J. F. (2003). Microbial communities in acid mine drainage. FEMS Microbiol. Ecol. 44, 139–152. doi: 10.1016/S0168-6496(03)00028-X
Baker-Austin, C., and Dopson, M. (2007). Life in acid: pH homeostasis in acidophiles. Trends Microbiol. 15, 165–171. doi: 10.1016/j.tim.2007.02.005
Ballouz, S., Francis, A. R., Lan, R., and Tanaka, M. M. (2010). Conditions for the evolution of gene clusters in bacterial genomes. PLoS Computat. Biol. 6, e1000672. doi: 10.1371/journal.pcbi.1000672
Beard, S., Ossandon, F. J., Rawlings, D. E., and Quatrini, R. (2021). The flexible genome of acidophilic prokaryotes. Curr. Issues Molec. Biol. 40, 231–266. doi: 10.21775/cimb.040.231
Bearson, S., Bearson, B., and Foster, J. W. (1997). Acid stress responses in enterobacteria. FEMS Microbiol. Lett. 147, 173–180. doi: 10.1111/j.1574-6968.1997.tb10238.x
Belin, B. J., Busset, N., Giraud, E., Molinaro, A., Silipo, A., and Newman, D. K. (2018). Hopanoid lipids: from membranes to plant–bacteria interactions. Nat. Rev. Microbiol. 16, 304–315. doi: 10.1038/nrmicro.2017.173
Benison, K. C., O'Neill, W. K., Blain, D., and Hallsworth, J. E. (2021). Water activities of acid brine lakes approach the limit for life. Astrobiology 21, 729–740. doi: 10.1089/ast.2020.2334
Bergmann, G. T., Bates, S. T., Eilers, K. G., Lauber, C. L., Caporaso, J. G., Walters, W. A., et al. (2011). The under-recognized dominance of Verrucomicrobia in soil bacterial communities. Soil Biol. Biochem. 43, 1450–1455. doi: 10.1016/j.soilbio.2011.03.012
Blum, M., Chang, H. Y., Chuguransky, S., Grego, T., Kandasaamy, S., Mitchell, A., et al. (2021). The InterPro protein families and domains database: 20 years on. Nucleic Acids Res. 49, D344–D354. doi: 10.1093/nar/gkaa977
Boase, K., González, C., Vergara, E., Neira, G., Holmes, D., and Watkin, E. (2022). Prediction and inferred evolution of acid tolerance genes in the biotechnologically important acidihalobacter genus. Front. Microbiol. 13. doi: 10.3389/fmicb.2022.848410
Buchfink, B., Xie, C., and Huson, D. H. (2015). Fast and sensitive protein alignment using DIAMOND. Nat. Methods 12, 59–60. doi: 10.1038/nmeth.3176
Buetti-Dinh, A., Dethlefsen, O., Friedman, R., and Dopson, M. (2016). Transcriptomic analysis reveals how a lack of potassium ions increases Sulfolobus acidocaldarius sensitivity to pH changes. Microbiol. 162, 1422–1434. doi: 10.1099/mic.0.000314
Bury-Moné, S., Mendz, G. L., Ball, G. E., Thibonnier, M., Stingl, K., Ecobichon, C., et al. (2008). Roles of α and β carbonic anhydrases of Helicobacter pylori in the urease-dependent response to acidity and in colonization of the murine gastric mucosa. Infect. Immun. 76, 497–509. doi: 10.1128/IAI.00993-07
Cabello-Yeves, P. J., Ghai, R., Mehrshad, M., Picazo, A., Camacho, A., and Rodriguez-Valera, F. (2017). Reconstruction of diverse verrucomicrobial genomes from metagenome datasets of freshwater reservoirs. Front. Microbiol. 8, 2131. doi: 10.3389/fmicb.2017.02131
Cai, Y., Yun, J., and Jia, Z. (2022). Phylogeny and metabolic potential of the methanotrophic lineage MO3 in beijerinckiaceae from the paddy soil through metagenome-assembled genome reconstruction. Microorganisms 10, 955. doi: 10.3390/microorganisms10050955
Cárdenas, J. P., Quatrini, R., and Holmes, D. S. (2016). Genomic and metagenomic challenges and opportunities for bioleaching: a mini review. Res. Microbiol. 167, 529–538. doi: 10.1016/j.resmic.2016.06.007
Carere, C. R., Hards, K., Houghton, K. M., Power, J. F., McDonald, B., Collet, C., et al. (2017). Mixotrophy drives niche expansion of verrucomicrobial methanotrophs. ISME J. 11, 2599–2610. doi: 10.1038/ismej.2017.112
Cholo, M. C., van Rensburg, E. J., Osman, A. G., and Anderson, R. (2015). Expression of the genes encoding the Trk and Kdp potassium transport systems of Mycobacterium tuberculosis during growth in vitro. BioMed Res. Int. 2015, 1–11. doi: 10.1155/2015/608682
Christel, S., Herold, M., Bellenberg, S., Buetti-Dinh, A., El Hajjami, M., Pivkin, I. V., et al. (2018). Weak iron oxidation by Sulfobacillus thermosulfidooxidans maintains a favorable redox potential for chalcopyrite bioleaching. Front. Microbiol. 9, 3059. doi: 10.3389/fmicb.2018.03059
Chun, J., Oren, A., Ventosa, A., Christensen, H., Arahal, D. R., da Costa, M. S., et al. (2018). Proposed minimal standards for the use of genome data for the taxonomy of prokaryotes. Int. J. System. Evolution. Microbiol. 68, 461–466. doi: 10.1099/ijsem.0.002516
Cortez, D., Neira, G., González, C., Vergara, E., and Holmes, D. S. (2022). A large-scale genome-based survey of acidophilic bacteria suggests that genome streamlining is an adaption for life at low pH. Front. Microbiology. 13:803241. doi: 10.3389/fmicb.2022.803241
Cotruvo, J. A. (2019). The chemistry of lanthanides in biology: recent discoveries, emerging principles, and technological applications. ACS Central Sci. 5, 1496–1506. doi: 10.1021/acscentsci.9b00642
Csuös, M. (2010). Count: evolutionary analysis of phylogenetic profiles with parsimony and likelihood. Bioinformatics 26, 1910–1912. doi: 10.1093/bioinformatics/btq315
Damiano, M. A., Bastianelli, D., Al Dahouk, S., Köhler, S., Cloeckaert, A., De Biase, D., et al. (2015). Glutamate decarboxylase-dependent acid resistance in Brucella spp.: distribution and contribution to fitness under extremely acidic conditions. Appl. Environ. Microbiol. 81, 578–586. doi: 10.1128/AEM.02928-14
Darling, A. C., Mau, B., Blattner, F. R., and Perna, N. T. (2004). Mauve: multiple alignment of conserved genomic sequence with rearrangements. Genome Res. 14, 1394–1403. doi: 10.1101/gr.2289704
Darling, A. E., Miklós, I., and Ragan, M. A. (2008). Dynamics of genome rearrangement in bacterial populations. PLoS Genet. 4, e1000128. doi: 10.1371/journal.pgen.1000128
de Vos, W. M. (2017). Microbe Profile: Akkermansia muciniphila: a conserved intestinal symbiont that acts as the gatekeeper of our mucosa. Microbiol. 163, 646–648. doi: 10.1099/mic.0.000444
Diskowski, M., Mikusevic, V., Stock, C., and Hänelt, I. (2015). Functional diversity of the superfamily of K+ transporters to meet various requirements. Biol. Chem. 396, 1003–1014. doi: 10.1515/hsz-2015-0123
Dopson, M. (2016). “Physiological and phylogenetic diversity of acidophilic bacteria,” in Acidophiles: Life in Extremely Acidic Environments. 79–92. doi: 10.21775/9781910190333.05
Dopson, M., and Holmes, D. S. (2014). Metal resistance in acidophilic microorganisms and its significance for biotechnologies. Appl. Microbiol. Biotechnol. 98, 8133–8144. doi: 10.1007/s00253-014-5982-2
Dorus, S., Mimura, H., and Epstein, W. (2001). Substrate-binding clusters of the K+-transporting Kdp ATPase of Escherichia coli investigated by amber suppression scanning mutagenesis. J. Biol. Chem. 276, 9590–9598. doi: 10.1074/jbc.M009365200
Dunfield, P. F., Yuryev, A., Senin, P., Smirnova, A. V., Stott, M. B., Hou, S., et al. (2007). Methane oxidation by an extremely acidophilic bacterium of the phylum Verrucomicrobia. Nature 450, 879–882. doi: 10.1038/nature06411
Emms, D. M., and Kelly, S. (2018). OrthoFinder2: fast and accurate phylogenomic orthology analysis from gene sequences. BioRxiv 466201.
Epstein, W. (2003). The roles and regulation of potassium in bacteria. Progr. Nucleic Acid Res. Molec. Biol. 75, 293–320. doi: 10.1016/S0079-6603(03)75008-9
Erikstad, H. A., Ceballos, R. M., Smestad, N. B., and Birkeland, N. K. (2019). Global biogeographic distribution patterns of thermoacidophilic Verrucomicrobia methanotrophs suggest allopatric evolution. Front. Microbiol. 10, 1129. doi: 10.3389/fmicb.2019.01129
Featherston, E. R., and Cotruvo, J. A. (2021). The biochemistry of lanthanide acquisition, trafficking, and utilization. Biochimica et Biophysica Acta (BBA)-Molec. Cell Res. 1868, 118864. doi: 10.1016/j.bbamcr.2020.118864
Feng, S., Qiu, Y., Huang, Z., Yin, Y., Zhang, H., Zhu, D., et al. (2021). The adaptation mechanisms of Acidithiobacillus caldus CCTCC M 2018054 to extreme acid stress: bioleaching performance, physiology, and transcriptomics. Environ. Res. 199, 111341. doi: 10.1016/j.envres.2021.111341
Finn, R. D., Bateman, A., Clements, J., Coggill, P., Eberhardt, R. Y., Eddy, S. R., et al. (2014). Pfam: the protein families database. Nucleic Acids Res. 42, D222–D230. doi: 10.1093/nar/gkt1223
Foster, J. W. (2004). Escherichia coli acid resistance: tales of an amateur acidophile. Nat. Rev. Microbiol. 2, 898–907. doi: 10.1038/nrmicro1021
Freitas, S., Hatosy, S., Fuhrman, J. A., Huse, S. M., Welch, D. B. M., Sogin, M. L., et al. (2012). Global distribution and diversity of marine Verrucomicrobia. The ISME J. 6, 1499–1505. doi: 10.1038/ismej.2012.3
González-Rosales, C., Vergara, E., Dopson, M., Valdés, J. H., and Holmes, D. S. (2022). Integrative genomics sheds light on evolutionary forces shaping the acidithiobacillia class acidophilic lifestyle. Front. Microbiol. 12:822229. doi: 10.3389/fmicb.2021.822229
González-Toril, E., Martínez-Frías, J., Gómez Gómez, J. M., Rull, F., and Amils, R. (2005). Iron meteorites can support the growth of acidophilic chemolithoautotrophic microorganisms. Astrobiology 5, 406–414. doi: 10.1089/ast.2005.5.406
Gumerov, V. M., and Zhulin, I. B. (2020). TREND: a platform for exploring protein function in prokaryotes based on phylogenetic, domain architecture and gene neighborhood analyses. Nucl. Acids Res. 48, W72–W76. doi: 10.1093/nar/gkaa243
Gumerov, V. M., and Zhulin, I. B. (2022). Correction to ‘TREND: a platform for exploring protein function in prokaryotes based on phylogenetic, domain architecture and gene neighborhood analyses'. Nucl. Acids Res. 50, 1795. doi: 10.1093/nar/gkac034
Hedrich, S., and Schippers, A. (2021). Distribution of acidophilic microorganisms in natural and man-made acidic environments. Current Issues Molec. Biol. 40, 25–48. doi: 10.21775/cimb.040.025
Hou, S., Makarova, K. S., Saw, J. H., Senin, P., Ly, B. V., Zhou, Z., et al. (2008). Complete genome sequence of the extremely acidophilic methanotroph isolate V4, Methylacidiphilum infernorum, a representative of the bacterial phylum Verrucomicrobia. Biol. Direct 3, 26. doi: 10.1186/1745-6150-3-26
Hu, W., Feng, S., Tong, Y., Zhang, H., and Yang, H. (2020). Adaptive defensive mechanism of bioleaching microorganisms under extremely environmental acid stress: Advances and perspectives. Biotechnol. Adv. 42, 107580. doi: 10.1016/j.biotechadv.2020.107580
Huang, Q., Jiang, H., Briggs, B. R., Wang, S., Hou, W., Li, G., et al. (2013). Archaeal and bacterial diversity in acidic to circumneutral hot springs in the Philippines. FEMS Microbiol. Ecol. 85, 452–464. doi: 10.1111/1574-6941.12134
Hutkins, R. W., and Nannen, N. L. (1993). Ph homeostasis in lactic acid bacteria1. J. Dairy Sci. 76, 2354–2365. doi: 10.3168/jds.S0022-0302(93)77573-6
Islam, T., Jensen, S., Reigstad, L. J., Larsen, Ø., and Birkeland, N. K. (2008). Methane oxidation at 55 C and pH 2 by a thermoacidophilic bacterium belonging to the Verrucomicrobia phylum. Proc. Nat. Acad. Sci. 105, 300–304. doi: 10.1073/pnas.0704162105
Jain, C., Rodriguez, -R L. M, Phillippy, A. M., Konstantinidis, K. T., and Aluru, S. (2018). High throughput ANI analysis of 90K prokaryotic genomes reveals clear species boundaries. Nat. Commun. 9, 1–8. doi: 10.1038/s41467-018-07641-9
Johnson, D. B., and Aguilera, A. (2015). “The microbiology of extremely acidic environments,” in Manual of Environmental Microbiology 4th Edition, chapter 4.3.1. M. V. Yates, H. Cindy, C. H. Nakatsu, R. V. Miller and S. D. Pillai, eds. New York: ASM Press. doi: 10.1128/9781555818821.ch4.3.1
Johnson, D. B., and Hallberg, K. B. (2005). Acid mine drainage remediation options: a review. Sci Total Environ. 338, 3–14. doi: 10.1016/j.scitotenv.2004.09.002
Jones, P., Binns, D., Chang, H. Y., Fraser, M., Li, W., McAnulla, C., et al. (2014). InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236–1240. doi: 10.1093/bioinformatics/btu031
Kalyaanamoorthy, S., Minh, B. Q., Wong, T. K., Von Haeseler, A., and Jermiin, L. S. (2017). ModelFinder: fast model selection for accurate phylogenetic estimates. Nat. Methods 14, 587–589. doi: 10.1038/nmeth.4285
Kannenberg, E. L., and Poralla, K. (1999). Hopanoid biosynthesis and function in bacteria. Naturwissenschaften 86, 168–176. doi: 10.1007/s001140050592
Kans, J. (2020). “Entrez direct: E-utilities on the UNIX command line,” in Entrez Programming Utilities Help. National Center for Biotechnology Information (US).
Katoh, K., Kuma, K. I., Toh, H., and Miyata, T. (2005). MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 33, 511–518. doi: 10.1093/nar/gki198
Kelly, S., and Maini, P. K. (2013). DendroBLAST: approximate phylogenetic trees in the absence of multiple sequence alignments. PLoS ONE 8, e58537. doi: 10.1371/journal.pone.0058537
Kim, H. W., Kim, N. K., Phillips, A. P., Parker, D. A., Liu, P., Whitaker, R. J., et al. (2022). Genome sequence of a thermoacidophilic methanotroph belonging to the verrucomicrobiota phylum from geothermal hot springs in yellowstone national park: a metagenomic assembly and reconstruction. Microorganisms 10, 142. doi: 10.3390/microorganisms10010142
Konings, W. N., Albers, S. V., Koning, S., and Driessen, A. J. (2002). The cell membrane plays a crucial role in survival of bacteria and archaea in extreme environments. Antonie Van Leeuwenhoek 81, 61–72. doi: 10.1023/A:1020573408652
Kruse, T., Ratnadevi, C. M., Erikstad, H. A., and Birkeland, N. K. (2019). Complete genome sequence analysis of the thermoacidophilic verrucomicrobial methanotroph “Candidatus Methylacidiphilum kamchatkense” strain Kam1 and comparison with its closest relatives. BMC Genom 20, 1–15. doi: 10.1186/s12864-019-5995-4
Letunic, I., and Bork, P. (2019). Interactive Tree Of Life (iTOL) v4: recent updates and new developments. Nucleic Acids Res. 47, W256–W259. doi: 10.1093/nar/gkz239
Liu, Y., Tang, H., Lin, Z., and Xu, P. (2015). Mechanisms of acid tolerance in bacteria and prospects in biotechnology and bioremediation. Biotechnol. Adv. 33, 1484–1492. doi: 10.1016/j.biotechadv.2015.06.001
Mangold, S., Rao Jonna, V., and Dopson, M. (2013). Response of Acidithiobacillus caldus toward suboptimal pH conditions. Extremophiles 17, 689–696. doi: 10.1007/s00792-013-0553-5
Masrati, G., Dwivedi, M., Rimon, A., Gluck-Margolin, Y., Kessel, A., Ashkenazy, H., et al. (2018). Broad phylogenetic analysis of cation/proton antiporters reveals transport determinants. Nat. Commun. 9, 1–14. doi: 10.1038/s41467-018-06770-5
Mates, A. K., Sayed, A. K., and Foster, J. W. (2007). Products of the Escherichia coli acid fitness island attenuate metabolite stress at extremely low pH and mediate a cell density-dependent acid resistance. J. Bacteriol. 189, 2759–2768. doi: 10.1128/JB.01490-06
Meier-Kolthoff, J. P., Auch, A. F., Klenk, H. P., and Göker, M. (2013). Genome sequence-based species delimitation with confidence intervals and improved distance functions. BMC Bioinform.14, 1–14. doi: 10.1186/1471-2105-14-60
Minh, B. Q., Nguyen, M. A. T., and von Haeseler, A. (2013). Ultrafast approximation for phylogenetic bootstrap. Molec. Biol. Evolut. 30, 1188–1195. doi: 10.1093/molbev/mst024
Mirete, S., Morgante, V., and González-Pastor, J. E. (2017). “Acidophiles: diversity and mechanisms of adaptation to acidic environments,” in Adaption of Microbial Life to Environmental Extremes. p. 227–251. doi: 10.1007/978-3-319-48327-6_9
Mohammadi, S. S., Pol, A., van Alen, T., Jetten, M. S., and Op den Camp, H. J. (2017). Ammonia oxidation and nitrite reduction in the verrucomicrobial methanotroph Methylacidiphilum fumariolicum SolV. Front. Microbiol. 8, 1901. doi: 10.3389/fmicb.2017.01901
Mohammadi, S. S., Schmitz, R. A., Pol, A., Berben, T., Jetten, M. S., and Op den Camp, H. J. (2019). The acidophilic methanotroph Methylacidimicrobium tartarophylax 4AC grows as autotroph on H2 under microoxic conditions. Front. Microbiol. 10, 2352. doi: 10.3389/fmicb.2019.02352
Neira, G., Cortez, D., Jil, J., and Holmes, D. S. (2020). AciDB 1.0: a database of acidophilic organisms, their genomic information and associated metadata. Bioinform. 36, 4970–4971. doi: 10.1093/bioinformatics/btaa638
Nguyen, L. T., Schmidt, H. A., Von Haeseler, A., and Minh, B. Q. (2015). IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies. Molec. Biol. Evolut. 32, 268–274. doi: 10.1093/molbev/msu300
Padan, E., Tzubery, T., Herz, K., Kozachkov, L., Rimon, A., and Galili, L. (2004). NhaA of Escherichia coli, as a model of a pH-regulated Na+/H+ antiporter. Biochim. Biophys. Acta (BBA)-Bioenerget. 1658, 2–13. doi: 10.1016/j.bbabio.2004.04.018
Parks, D. H., Chuvochina, M., Rinke, C., Mussig, A. J., Chaumeil, P. A., and Hugenholtz, P. (2022). GTDB: an ongoing census of bacterial and archaeal diversity through a phylogenetically consistent, rank normalized and complete genome-based taxonomy. Nucleic Acids Res. 50, D785–D794. doi: 10.1093/nar/gkab776
Parks, D. H., Imelfort, M., Skennerton, C. T., Hugenholtz, P., and Tyson, G. W. (2015). CheckM: assessing the quality of microbial genomes recovered from isolates, single cells, and metagenomes. Genome Res. 25, 1043–1055. doi: 10.1101/gr.186072.114
Peng, T., Ma, L., Feng, X., Tao, J., Nan, M., Liu, Y., et al. (2017). Genomic and transcriptomic analyses reveal adaptation mechanisms of an Acidithiobacillus ferrivorans strain YL15 to alpine acid mine drainage. PLoS One 12, e0178008. doi: 10.1371/journal.pone.0178008
Picone, N., Blom, P., Hogendoorn, C., Frank, J., van Alen, T., Pol, A., et al. (2021b). Metagenome assembled genome of a novel verrucomicrobial methanotroph from Pantelleria Island. Front. Microbiol. 12, 1176. doi: 10.3389/fmicb.2021.666929
Picone, N., Blom, P., Wallenius, A. J., Hogendoorn, C., Mesman, R., Cremers, G., et al. (2021a). Methylacidimicrobium thermophilum AP8, a novel methane-and hydrogen-oxidizing bacterium isolated from volcanic soil on pantelleria Island, Italy. Front. Microbiol. 12, 225. doi: 10.3389/fmicb.2021.637762
Pol, A., Barends, T. R., Dietl, A., Khadem, A. F., Eygensteyn, J., Jetten, M. S., et al. (2014). Rare earth metals are essential for methanotrophic life in volcanic mudpots. Environ. Microbiol. 16, 255–264. doi: 10.1111/1462-2920.12249
Pol, A., Heijmans, K., Harhangi, H. R., Tedesco, D., Jetten, M. S., and Den Camp, H. J. O. (2007). Methanotrophy below pH 1 by a new Verrucomicrobia species. Nature 450, 874–878. doi: 10.1038/nature06222
Price, M. N., Dehal, P. S., and Arkin, A. P. (2010). FastTree 2–approximately maximum-likelihood trees for large alignments. PLoS ONE 5, e9490. doi: 10.1371/journal.pone.0009490
Qian, W., and Zhang, J. (2014). Genomic evidence for adaptation by gene duplication. Genome Res. 24, 1356–1362. doi: 10.1101/gr.172098.114
Quatrini, R., and Johnson, D. B. (2018). Microbiomes in extremely acidic environments: functionalities and interactions that allow survival and growth of prokaryotes at low pH. Current Opin. Microbiol. 43, 139–147. doi: 10.1016/j.mib.2018.01.011
Richard, H., and Foster, J. W. (2004). Escherichia coli glutamate-and arginine-dependent acid resistance systems increase internal pH and reverse transmembrane potential. J. Bacteriol. 186, 6032–6041. doi: 10.1128/JB.186.18.6032-6041.2004
Schippers, A., Hedrich, S., Vasters, J., Drobe, M., Sand, W., and Willscher, S. (2013). Biomining: metal recovery from ores with microorganisms. Geobiotechnology I 141, 1–47. doi: 10.1007/10_2013_216
Schmitz, R. A., Peeters, S. H., Versantvoort, W., Picone, N., Pol, A., Jetten, M. S., et al. (2021). Verrucomicrobial methanotrophs: ecophysiology of metabolically versatile acidophiles. FEMS Microbiol. Rev. 45, fuab007. doi: 10.1093/femsre/fuab007
Schulz, S., Wilkes, M., Mills, D. J., Kühlbrandt, W., and Meier, T. (2017). Molecular architecture of the N-type ATP ase rotor ring from Burkholderia pseudomallei. EMBO Rep. 18, 526–535. doi: 10.15252/embr.201643374
Sharma, A., Kawarabayasi, Y., and Satyanarayana, T. (2012). Acidophilic bacteria and archaea: acid stable biocatalysts and their potential applications. Extremophiles 16, 1–19. doi: 10.1007/s00792-011-0402-3
Sharp, C. E., Smirnova, A. V., Graham, J. M., Stott, M. B., Khadka, R., Moore, T. R., et al. (2014). Distribution and diversity of Verrucomicrobia methanotrophs in geothermal and acidic environments. Environ. Microbiol. 16, 1867–1878. doi: 10.1111/1462-2920.12454
Shen, W., Le, S., Li, Y., and Hu, F. (2016). SeqKit: a cross-platform and ultrafast toolkit for FASTA/Q file manipulation. PLoS ONE 11, e0163962. doi: 10.1371/journal.pone.0163962
Slonczewski, J. L., Fujisawa, M., Dopson, M., and Krulwich, T. A. (2009). Cytoplasmic pH measurement and homeostasis in bacteria and archaea. Adv. Microb. Physiol. 55, 1–317. doi: 10.1016/S0065-2911(09)05501-5
Sohlenkamp, C., and Geiger, O. (2016). Bacterial membrane lipids: diversity in structures and pathways. FEMS Microbiol. Rev. 40, 133–159. doi: 10.1093/femsre/fuv008
Sowell, S. M., Wilhelm, L. J., Norbeck, A. D., Lipton, M. S., Nicora, C. D., Barofsky, D. F., et al. (2009). Transport functions dominate the SAR11 metaproteome at low-nutrient extremes in the Sargasso Sea. ISME J. 3, 93–105. doi: 10.1038/ismej.2008.83
Sun, Y. (2016). “F 1 F 0-ATPase functions under markedly acidic conditions in bacteria,” in Regulation of Ca2+-ATPases, V-ATPases and F-ATPases. p. 459–468. Cham: Springer. doi: 10.1007/978-3-319-24780-9_22
Sun, Y., Fukamachi, T., Saito, H., and Kobayashi, H. (2012). Respiration and the F1Fo-ATPase enhance survival under acidic conditions in Escherichia coli. PLoS ONE 7, e52577. doi: 10.1371/journal.pone.0052577
Trchounian, A., and Kobayashi, H. (1999). Kup is the major K+ uptake system in Escherichia coli upon hyper-osmotic stress at a low pH. FEBS Lett. 447, 144–148. doi: 10.1016/S0014-5793(99)00288-4
Tria, F. D. K., Landan, G., and Dagan, T. (2017). Phylogenetic rooting using minimal ancestor deviation. Nat. Ecol. Evolut. 1, 1–7. doi: 10.1038/s41559-017-0193
Tucker, D. L., Tucker, N., and Conway, T. (2002). Gene expression profiling of the pH response in Escherichia coli. J. Bacteriol. 184, 6551–6558. doi: 10.1128/JB.184.23.6551-6558.2002
Valdés, J., Pedroso, I., Quatrini, R., and Holmes, D. S. (2008). Comparative genome analysis of Acidithiobacillus ferrooxidans, A. thiooxidans and A. caldus: insights into their metabolism and ecophysiology. Hydrometallurgy 94, 180–184. doi: 10.1016/j.hydromet.2008.05.039
van Teeseling, M. C., Pol, A., Harhangi, H. R., van der Zwart, S., Jetten, M. S., Op den Camp, H. J., et al. (2014). Expanding the verrucomicrobial methanotrophic world: description of three novel species of Methylacidimicrobium gen. nov. Appl. Environ. Microbiol. 80, 6782–6791. doi: 10.1128/AEM.01838-14
Vergara, E., Neira, G., González, C., Cortez, D., Dopson, M., and Holmes, D. S. (2020). Evolution of predicted acid resistance mechanisms in the extremely acidophilic Leptospirillum genus. Genes 11, 389. doi: 10.3390/genes11040389
Welander, P. V., Hunter, R. C., Zhang, L., Sessions, A. L., Summons, R. E., and Newman, D. K. (2009). Hopanoids play a role in membrane integrity and pH homeostasis in Rhodopseudomonas palustris TIE-1. J. Bacteriol. 191, 6145–6156. doi: 10.1128/JB.00460-09
Yan, Z., Li, X., Li, M., Zuo, W., and Shan, S. (2018). “Shift-net: Image inpainting via deep feature rearrangement,” in Proceedings of the European Conference on Computer Vision (ECCV). p. 1–17. doi: 10.1007/978-3-030-01264-9_1
Zhang, H., Pan, Y., Hu, L., Hudson, M. A., Hofstetter, K. S., Xu, Z., et al. (2020a). TrkA undergoes a tetramer-to-dimer conversion to open TrkH which enables changes in membrane potential. Nat. Commun. 11, 1–11. doi: 10.1038/s41467-019-14240-9
Zhang, Y., Zhang, S., Zhao, D., Ni, Y., Wang, W., and Yan, L. (2020b). Complete genome sequence of Acidithiobacillus ferrooxidans YNTRS-40, a strain of the ferrous iron-and sulfur-oxidizing acidophile. Microorganisms 8, 2. doi: 10.3390/microorganisms8010002
Zhu, Q., Kosoy, M., and Dittmar, K. (2014). HGTector: an automated method facilitating genome-wide discovery of putative horizontal gene transfers. BMC Genom. 15, 1–18. doi: 10.1186/1471-2164-15-717
Keywords: geothermal ecosystems, acid resistance mechanisms in hyperacidic environments, polyextremophiles, Verrucomicrobia phylogenomics, methanotrophy, lanthanides, horizontal gene transfer (HGT), biomining
Citation: Neira G, Vergara E and Holmes DS (2022) Genome-guided prediction of acid resistance mechanisms in acidophilic methanotrophs of phylogenetically deep-rooted Verrucomicrobia isolated from geothermal environments. Front. Microbiol. 13:900531. doi: 10.3389/fmicb.2022.900531
Received: 20 March 2022; Accepted: 22 July 2022;
Published: 23 September 2022.
Edited by:
María Sofía Urbieta, CONICET - UNLP, ArgentinaReviewed by:
Irene Sánchez-Andrea, Wageningen University and Research, NetherlandsGloria Paz Levicán, University of Santiago, Chile
Copyright © 2022 Neira, Vergara and Holmes. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: David S. Holmes, ZHNob2xtZXMyMDAwJiN4MDAwNDA7eWFob28uY29t
†These authors have contributed equally to this work