The instant application contains a Sequence Listing which has been submitted electronically in ASCII format and is hereby incorporated by reference in its entirety. Said ASCII copy, created on Feb. 17, 2022, is named 018617_01341_SL.txt and is 979,969 bytes in size.
Rare earth elements (REE) are essential for the manufacturing of modern electronics, sustainable energy technologies including electric motors and wind turbine generators; solid state lighting; battery anodes; high-temperature superconductors; and high-strength lightweight alloys. All of these applications place increasing demands on the global REE supply chain. As the world demand for sustainable energy grows, finding a reliable and sustainable source of REE is critical.
Current methods for refining REE often involve harsh chemicals, high temperatures, high pressures and generate a considerable amount of toxic waste. These processes give sustainable energy technologies reliant on REE a high environmental and carbon footprint. As a consequence, due to its high environmental standards, the United States has no capacity to produce purified REE.
There is growing interest in biological methods to supplement, if not completely replace traditional REE extraction and purification methods. Bioleaching is used to extract 5% of the world's gold, and ≈15% of the world's copper supply, and biomining in Chile alone accounts for 10% of the world's Cu supply.
The performance of REE-bioleaching lags behind thermochemical processes. For example, while thermochemical methods have 89-98% REE extraction efficiency from monazite ore, Aspergillus species can only achieve≈3-5%. The acid-producing microbe Gluconobacter oxydans B58 can recover≈50% of REE from FCC catalysts. However, techno-economic analysis indicates that even this extraction efficiency is still not high enough for commercial viability.
Recent efforts to improve bioleaching have focused exclusively on process optimization. It is believed that no previous genetic approaches have yet been taken for any bioleaching microbe. With recent advances in tools for reading and writing genomes, genetic engineering is an attractive solution for enhancing bioleaching. However, applying these tools to non-model microorganisms like G. oxydans can be a significant challenge. While there have been some advances for editing the genome of G. oxydans it has remained unknown where the genome can be edited to improve bioleaching results. Thus, there is an ongoing and unmet need for improved compositions, engineered organisms, and methods for separating REEs from compositions that contain them. The present disclosure is pertinent to this need.
The present disclosure provides a description of a whole genome knockout collection for Gluconobacter oxydans B58, and use of it to comprehensively characterize the genomics of rare earth elements (REEs) bioleaching. In total, 304 genes that notably alter production of G. oxydans' acidic bio-lixiviant, including 165 that make statistically significant changes, were identified. Based in part on this analysis, the present disclosure provides modified bacteria for use in bioleaching REEs. The modified bacteria comprise at least one engineered genetic change that is correlated with improved bioleaching of the REEs, relative to REE bioleaching by unmodified bacteria of the same species as the modified bacteria. The at least one genetic change results in decreased expression, or increased expression, of at least one gene. In non-limiting embodiments, at least one gene for which expression is modified encodes a protein that participates in phosphate-specific transport system signaling, or encodes a protein that participates in pyrroloquinoline quinone (PQQ) synthesis. In non-limiting examples, expression of a gene that encodes a protein that participates in the phosphate-specific transport system signaling is suppressed. In certain embodiments, the suppressed gene is pstS, pstB or pstC. In certain embodiments, a gene that encodes a protein that participates in the PQQ synthesis is increased. In non-limiting embodiments, the expression of at least one of the genes pqqA, pqqB, pqqC, pqqD, pqqE, tldD and tldE, is increased. In addition to these and other genetic modifications described herein, the modified bacteria exhibit increase expression of mgdh relative to expression of mgdh by unmodified bacteria. In certain embodiments, expression of pstS, pstB, pstC, or a combination thereof is reduced, or expression of pqqA, pqqB, pqqC, pqqD, pqqE, tldD, tldE, or a combination thereof is increased. In these contexts expression of mgdh may also be increased.
In another aspect, the disclosure provides for contacting a composition comprising the REEs with a composition produced by the described modified bacteria. The composition produced by the bacteria may be considered a lixiviant, or a biolixiviant because it is produced by the described bacteria. The disclosure provides separating REEs from the composition after contacting the composition with the biolixiviant. The separated REEs are suitable for use in a wide range of applications that will be apparent to those skilled in the art.
In another aspect, the disclosure provides kits that contain one or more sealable containers in which the described modified bacteria are held. The kits may further comprise printed material, such as instructions for use of the modified bacteria to form a biolixiviant, and/or to extract REEs from a composition where they are present.
Unless defined otherwise herein, all technical and scientific terms used in this disclosure have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure pertains.
Every numerical range given throughout this specification includes its upper and lower values, as well as every narrower numerical range that falls within it, as if such narrower numerical ranges were all expressly written herein.
The disclosure includes all polynucleotide and amino acid sequences described herein. Each RNA sequence includes its DNA equivalent, and each DNA sequence includes its RNA equivalent. Complementary and anti-parallel polynucleotide sequences are included. Every DNA and RNA sequence encoding polypeptides disclosed herein is encompassed by this disclosure. Amino acids of all protein sequences and all polynucleotide sequences encoding them are also included, including but not limited to sequences included by way of sequence alignments. Sequences of from 80.00%-99.99% identical to any sequence (amino acids and nucleotide sequences) of this disclosure are included.
The disclosure includes all polynucleotide and all amino acid sequences that are identified herein by way of a database entry. Such sequences are incorporated herein as they exist in the database on the effective filing date of this application or patent.
The disclosure includes modified microorganisms having any modified single gene, and modifications of all combinations of genes described herein in the text, figures, figure legends, and tables of this disclosure.
Any gene and any combination of the genes that are described herein may be excluded from the claims of this disclosure. In embodiments, a modified microorganism of the disclosure may comprise or consist of only one modification of a single gene. In embodiments, a modified microorganism of the disclosure may comprise or consist of any combination of gene modifications described herein. In embodiments, only one or only a combination of genes that influence bioleaching of REEs are modified.
In non-limiting embodiments, the disclosure provides modified bacteria in which the expression of at least one of the genes pqqA, pqqB, pqqC, pqqD, pqqE, tldD and tldE, is increased. In certain embodiments, expression of pstS, pstB, pstC, or a combination thereof is reduced, or expression of pqqA, pqqB, pqqC, pqqD, pqqE, tldD, tldE, or a combination thereof is increased. In certain embodiments, the modified bacteria exhibit increased expression of mgdh relative to expression of mgdh by unmodified bacteria, wherein the increased expression of mgdh is in the context of at least one other described genetic modification.
In embodiments, the modified bacteria comprises or consist of mutations that are selected from mutations in all of the genes listed in Table A, and including all numbers and ranges of numbers of genes between 1 gene and the total genes in Table A.
The disclosure includes modifications that disrupt one or a combination of genes, modifications that increase expression of one or a combination of genes, or a combination of modifications that decrease expression of one or more genes and modifications that increase expression of one or more genes. Thus, the modifications involve altering the expression of one or more genes. Increasing, e.g., overexpressing a gene, can be achieved using various techniques that will be apparent to those skilled in the art when given the benefit of the present disclosure. In embodiments, increasing expression of a gene is achieved by substituting an endogenous promoter with a promoter that increases expression of the gene, relative to expression of the gene that is produced by the endogenous promoter. By “substituting” a promoter it is meant that the endogenous promoter (e.g., the promoter that is ordinarily operatively linked to the gene of interest without genetic engineering) has been changed so that is does not drive expression of the gene in the modified bacteria, and therefore the substituted promoter drives gene expression. By making this change, more mRNA is transcribed, thus facilitating production of more protein encoded by the pertinent gene that is operatively linked to the promoter. Substituting a promoter can include inserting a new promoter, while leaving the endogenous promoter in place, or inserting the new promoter in place of the endogenous promoter. The promoter that is inserted so that it is operably linked to and therefore drives expression of the described gene(s) can be heterologous to the bacteria, meaning it is taken or derived from a different organism, or it may be endogenous to the organism but has been introduced into a new location such that it can drive expression of the described gene(s). Various prokaryotic promoters that are suitable for this purpose are known in the art and include, for example, tufa and tufB. The substituted promoter (e.g., the promoter that is introduced into the bacteria) may be a constitutive or inducible promoter. The substituted promoter may be a core promoter, a proximal promoter, or a distal promoter.
Representative and non-limiting embodiments of promoters that can be used to increase expression of one or more genes as described herein include:
As an alternative or in addition to promoter modification, the disclosure includes addition of and/or repositioning of enhancer elements to increase expression of the described gene(s).
As an alternative or in addition to changing promoters, the disclosure includes increasing copy number of the gene that is to be overexpressed. In embodiments, one or more copies of the gene can be inserted into a bacterial chromosome, or can be introduced into bacteria using a plasmid. A list of genes for which overexpression is encompassed by the disclosure is provided on Table A. The additional copies of the gene may be in tandem, such as in a polycistronic configuration, or may be separated by segments of the bacterial chromosome or plasmid. In embodiments, a composition comprising the described bacteria are modified by transformation using one or more plasmids, which may be configured to be replicated and transferred to other bacteria in a bacterial population, such as by horizontal transfer.
In another embodiment, the disclosure comprises decreasing expression of genes. Decreasing expression can be achieved using any suitable approach. In embodiments, decreasing expression comprises disrupting the gene such that the protein encoded by the gene is not produced, or a protein produced by the gene does not function in the same way as if it had not been modified. In embodiments, a protein that is encoded by a modified gene of this disclosure is produced but does not function to impede bioleaching of REEs from a composition comprising them. In embodiments, a modification of a gene comprises a knock-out of some or all of the gene. Modifications of the genes can be achieved using any suitable genetic engineering techniques. In non-limiting embodiments, the modification comprises an insertion, a deletion, or a combination thereof. The disclosure includes insertion within, or a deletion of any segment of a gene, including but not limited to a insertion or deletion of a single nucleotide, such that the encoded protein is not produced or its function is eliminated or reduced. In embodiments, an insertion replaces some or all of the described gene(s). In a non-limiting embodiment, the described gene(s) is modified by insertion of a transposable element. In non-limiting embodiments, the genes are modified using compositions and methods described in U.S. Pat. No. 11,053,493, from which the entire description is incorporated herein by reference. In embodiments, a modification of a gene comprises an insertion as described in Anzai, Isao A., et al. “Rapid curation of gene disruption collections using Knockout Sudoku.” Nature Protocols 12.10: 2110-2137 (2017), from which the entire disclosure is incorporated herein by reference. In alternative approaches, site specific nuclease, such as Cas nucleases, can be used to modify any of the described genes. In embodiments, a type I, type II or type III CRISPR system can be used. Thus, in embodiments, a guide-RNA directed nuclease can make any of the described modifications. In embodiments, recombination of a chromosome or plasmid can be used, such as by introducing a recombination template comprising additional copies of a gene, and/or a promoter, to facilitate recombination of the recombination template into a desired location. In an embodiment, homologous recombination is used, and as such, the recombination template includes left and right homology arms to specify the location of recombination. In embodiments, a transposon system can be used to interrupt a gene sequence, such as the Sleeping Beauty transposon system.
In embodiments, the modified bacteria comprise a modification of at least one gene described in
In embodiments, the disclosure comprises increasing expression of at least one gene described in Table A. In embodiments, the disclosure comprises decreasing expression of at least one gene described in Table A. In embodiments, the disclosure comprises increasing expression of at least one gene and decreasing expression of at least one gene described in Table A. In non-limiting embodiments, modified bacteria of this disclosure are modified such that they exhibit decreased expression of at least one of the following genes: GO_1415, pstA, pstB, pstC, pstS, ggtl, surA, petP, ykoH, speC, and tonB. In non-limiting embodiments, modified bacteria of this disclosure are modified such that they exhibit increased expression of at least mgdh, and/or genes involved in PQQ synthesis (e.g., pqqA, pqqB, pqqC, pqqD, pqqE, and tldD, also referred to as pqqABCDE as an operon, and tldE), or a combination thereof. In embodiments, any one or any combination of proteins expressed by the pqqA, pqqB, pqqC, pqqD, pqqE, and tldD genes can be modified to increase their activity, such as by modifying amino acids in an active site, or amino acids that improve structural stability, and the like.
Combinations of modifications that increase and decrease expression of genes are included in the disclosure. The disclosure also includes mixed populations of bacteria, wherein some of the members of the population have different genetic modifications than other members of the population.
In embodiments, the modified bacteria have at least one engineered genetic change that is correlated with improved bioleaching of REEs, relative to REE bioleaching by unmodified bacteria of the same species as the modified bacteria. The disclosure includes the proviso that the set of modified genes may exclude a disruption of membrane bound glucose dehydrogenase (mgdh) gene as the only modification of the described bacteria. However, this gene may also be disrupted, provided it is in the context of at least one other gene modification that is described herein. In non-limiting examples, at least one genetic change increases acidification of a medium in which the modified bacteria are present. In a non-limiting example, the at least one genetic change is in a gene that is part of a phosphate transport system. In embodiments, the bacteria are modified such that they comprise a mutated gene that comprises or consists of at least one of: GO_1415, pstA, pstB, pstC, pstS, ggtl, surA, petP, ykoH, speC, and tonB. In embodiments, the modification comprises a disruption of at least GO_1415, or pstC, or a combination thereof.
The disclosure includes compositions comprising one or more REEs and modified bacteria of the disclosure. The disclosure includes a biolixiviant produced by the modified bacteria and one or more REEs. In embodiments, the disclosure relates to separating combinations of REEs. In embodiments, the disclosure relates to separating any one or combination of lanthanum, cerium, praseodymium, neodymium, promethium, samarium, europium, gadolinium, terbium, dysprosium, holmium, erbium, thulium, ytterbium, lutetium, scandium, and yttrium, from a composition comprising one or more of the REEs. The composition comprising the REEs may be any composition of matter, including but not limited to solids, semi-solids, and liquids. In embodiments, the REEs are present in a feedstock. In non-limiting embodiments, the REEs are present in coal fly ash, virgin ore, electronic waste, fluid cracking catalysts, and the like.
The disclosure includes a method comprising contacting a composition comprising one or more types of REEs with a biolixiviant produced by modified bacteria of this disclosure. In an embodiment, the method further comprising separating and optionally purifying one or more types of REEs from the composition comprising the REEs and the biolixiviant.
The disclosure comprises isolated modified bacteria, cell cultures comprising the modified bacteria, and kits comprising the modified bacteria. In an embodiment, a kit comprises one or more sealed containers comprising the modified bacteria, which can be used in REE bioleaching approaches.
The disclosure includes media in which the bacteria are cultured, and bacterial secretions. In an embodiment, the disclosure provides a biolixiviant produced by the described bacteria. In embodiments, the kit contains a sealable or sealed container that contains a biolixiviant produced by the described bacteria. The disclosure also includes modifying bacteria so that they comprise at least one of the described gene modifications.
The disclosure includes all modified microorganisms described herein. The described approaches may be used to engineer any type of bacteria. In embodiments, the bacteria are Gram-negative bacteria. In embodiments, the bacteria are obligate aerobes. In embodiments, the bacteria modified as described herein comprise any member of the bacteria family Acetobacteraceae. In an embodiment, the bacteria is a type of Gluconobacter. In an embodiment, the modified bacteria are Gluconobacter oxydans. In this regard, G. oxydans secretes a biolixiviant rich in gluconic acid. This is produced by periplasmic glucose oxidation by the pyrroloquinoline quinone (PQQ)-dependent membrane-bound glucose dehydrogenase (mGDH). The final pH of the biolixiviant is a major factor in REE bioleaching. But, gluconic acid alone fails to explain bioleaching by G. oxydans: pure gluconic acid is far less effective at bioleaching than the biolixiviant produced by G. oxydans. This means that even the most previous successful efforts to up-regulate mGDH activity and gluconic acid production are unlikely to take full advantage of G. oxydans' biolixiviant production capabilities. Thus, the present disclosure reveals a curated set of genes that can be modified to improve REE extraction, as demonstrated in the following Examples.
To characterize the genome of G. oxydans and identify a comprehensive set of genes underlying its bioleaching capabilities, we built a carefully curated whole-genome knockout collection of single-gene transposon disruption mutants using Knockout Sudoku (
Development of a Knockout Collection for G. oxydans Covering 2,733 Genes
We built a saturating coverage transposon insertion mutant collection for G. oxydans B58 and catalogued and condensed it with the Knockout Sudoku combinatorial pooling method (
The progenitor collection catalog indicates that we were able to generate at least one disruption mutant for almost every non-essential gene in the G. oxydans genome. In total, we identified disruption strains for 2,733 genes out of the 3,283 genes in the G. oxydans B58 genome. Since every predicted gene contains at least seven AT or TA transposon insertion sites, the remaining 550 non-disrupted genes are likely to be essential. A Fisher's Exact Test for gene ontology (GO) enrichment representing 268 of the non-disrupted genes demonstrated significant enrichment in several essential ontologies, with the greatest enrichment in those relating to the ribosome and translation (
The progenitor collection catalog was used to create a condensed G. oxydans disruption collection with at least one representative per non-essential gene. 47 progenitor strains were verified by Sanger sequencing prior to condensing, of which 43 (92%) were confirmed to have the predicted transposon coordinate. We selected one mutant for all 2,733 disrupted genes, a second mutant for 2,354 genes, and a third mutant for 50 genes where mutant location information was poor. All mutants were struck out for single colonies, and 2-10 colonies per mutant were picked, depending on the predicted number of cross-contaminating disruption strains in the originating well. This condensed collection contains 17,706 mutants in 185 96-well plates.
The condensed collection catalog was validated by a second round of combinatorial pooling and sequencing. Of the 17,706 wells in the condensed collection, we were able to confirm the identity for 15,257. We confirmed 25 of these wells by Sanger sequencing, and 100% have the predicted transposon coordinate. Among these wells, we were able to verify the identity and location of 4,419 independent transposon insertion sites, representing a disruption mutant for 2,556 unique genes (
We screened the new G. oxydans B58 whole genome knockout collection for disruption mutants with differential acidification capability (
In total, we observed 304 genes that apparently controlled acidification (
We used gene ontology enrichment to determine which biological processes, metabolic functions, and cellular components the most significant gene disruption mutants are involved in (
Among the disrupted genes that led to a weaker acidity, several enrichment groups are related to the synthesis or use of the redox cofactor, PQQ, represented by pqqB, pqqC, pqqE, tldD, and mgdh (
Acidification rate is controlled by carbohydrate metabolism and respiration. Disruptions in the pentose phosphate pathway increase acidification rate (
We selected 14 strains with some of the most significantly increased or decreased acidification for further testing (
Additional disruptions that led to a more acidic biolixiviant included those in a hypothetical protein with no similarity to anything previously characterized (δGO_1415S); a gamma-glutamyltranspeptidase (δggtl); a periplasmic chaperone (δsurA); an HTH-type transcriptional regulator (δpetP); a two-component system sensor histidine kinase (δykoH); a Pyridoxal 5-phosphate (PLP)-dependent ornithine decarboxylase (δspeC); and a TPR domain protein that is a putative component of the TonB iron uptake system (δtonB).
9 of the tested strains produced biolixiviant significantly higher in pH than pWT (
We tested if 10 of the mutants that produced a more acidic biolixiviant could bioleach REE from retorted phosphor powder (RPP) from spent fluorescent lightbulbs more efficiently than pWT (
The δpstC mutant produced the most acidic biolixivant, and extracted the most REE from RPP: 5.5% total extraction efficiency as compared with pWT's 4.7%. Stated differently, δpstC removed 18% more REE from RPP than pWT. This increase in REE extraction remains significant even under a Bonferonni correction, the most stringent statistical test for significance. Without the adjustment, six of the better acidifiers were also better bioleachers than pWT (
Without intending to be bound by any particular theory, it is considered that disrupting the phosphate transport system de-represses acid production in G. oxydans. Six of the disruption strains that resulted in a lower biolixiviant pH (δpstC, δpstB, δggtl, δpstA, δpstS, and δykoH), including three that increased bioleaching (δpstC, δpstB, δggtl), are involved in phosphate transport, sensing and signaling.
In its natural environment, G. oxydans produces biolixiviants to liberate phosphate from minerals, not metals. Under phosphate-limiting conditions, the PstSCAB phosphate transporter will activate the histidine kinase, PhoR, which in turn phosphorylates the transcription factor PhoB, and activates the pho regulon, enabling phosphate assimilation and uptake. Under sufficient phosphate conditions, PhoB is deactivated by PhoR, which in turn inhibits expression of these genes. Without intending to be constrained by any particular view, it is considered that disrupting any of these genes prevents G. oxydans from sensing when it has released adequate phosphate and when to stop producing biolixiviants.
We also tested REE extraction by 4 mutants that produce a less acidic biolixiviant than pWT. Even under the most stringent statistical test, the Bonferonni correction, they were all worse bioleachers than pWT (
The δmgdh mutant was the worst bioleacher of all tested, considering its lack of gluconic-acid production. δmgdh reduced bioleaching by 97%. Disruption mutants that knocked out synthesis of mGDH's essential redox cofactor, PQQ, also produced significant reductions in biolixiviant acidity. δpqqC reduced bioleaching by ≈94%. While bioleaching by δmgdh and δpqqC was negligible compared to pWT, they were able to bioleach a statistically significant amount of REE compared to glucose alone. This indicates, that a bioleaching mechanism independent of mGDH exists in G. oxydans (
Disruption mutants in tldD and tldE were also much worse at bioleaching than pWT. δtldD reduces bioleaching by 92%, while δtldE reduces it by 63% (
It will be recognized from the foregoing description that bioleaching has the potential to revolutionize the environmental impact of REE production, and dramatically increase access to these critical ingredients for sustainable energy technology. The present disclosure related to this potential by providing for improved bioleaching by genetic engineering.
By constructing a whole genome knockout collection for G. oxydans, we are able to characterize the genetics of this process with high sensitivity and high completeness. In total we identified 165 gene disruption mutants that significantly change the acidity of its biolixiviant, rate of production, or both.
REE bioleaching by G. oxydans is predominantly controlled by two well-characterized systems: phosphate signaling and glucose oxidation that is supported by production of the redox cofactor PQQ. Interrupting phosphate signaling control of biolixiviant production by disrupting a single gene (pstC) can increase REE extraction by 18%. Disrupting the supply of the PQQ cofactor to the membrane bound glucose dehydrogenase reduces REE extraction by up to 92%.
Comprehensive screening of the G. oxydans genome also discovers completely new targets that contribute as much to REE bioleaching as previously characterized ones. For example, disrupting GO_1415, which encodes a protein of completely unknown function, increases REE bioleaching by 15%. Additionally, these results highlight the potential for a previously uncharacterized role of TldE in PQQ synthesis.
The discovery of the potential contribution of TldE to PQQ biosynthesis may allow for marked enhancement of the cofactor production through the additional over overexpression of this gene, and a consequent increase in dehydrogenase activity, including production of gluconic acid by mGDH and any useful downstream products. PQQ is an essential cofactor important for several other industrial applications of G. oxydans, including production of L-sorbose. Furthermore, PQQ alone has many applications across many biological processes from plant protection to neuron regeneration.
Without intending to be bound by any particular theory, it is believed that the present disclosure provides the first demonstration of improvement of bioleaching through genetic engineering. Furthermore, the creation of a whole-genome knockout collection in G. oxydans can facilitate its use as a model species for further studies in REE bioleaching and other industrially important applications of similar acetic acid bacteria. The findings of the two major systems contributing to acidification in G. oxydans according to this disclosure show that, for greatly improving bioleaching: reduce inhibition of regulation of acid production by disabling the phosphate-specific transport system, while over-expressing mgdh along with the expanded synthesis pathway for its cofactor PQQ.
Gluconobacter oxydans B58 Genome Sequencing
Gluconobacter oxydans strain NRRL B-58 (GoB58) was obtained from the American Type Culture Collection (ATTC), Manassas, VA. In all experiments, G. oxydans was cultured in yeast peptone mannitol media (YPM; 5 g L−1 yeast extract, 3 g L−1 peptone, 25 g L−1 mannitol), with or without antibiotic, as specified.
Genomic DNA was extracted from saturated culture using a Quick-DNA Miniprep kit from Zymo Research (Part number D3024, Irvine, CA). Genomic DNA library was prepared and sequenced using a TruSeq DNA PCR-Free Library Prep Kit (Illumina, San Diego, CA).
The prepared library was sequenced on a MiSeq Nano (Illumina, San Diego, CA, USA) with a 500 bp kit at the Cornell University Institute of Biotechnology (Ithaca, NY, USA). Resulting paired end reads were trimmed using Trimmomatic and assembled with SPAdes using k-mer sizes 21, 33, 55, 77, 99, and 127, and an auto coverage cutoff. Assembly quality was checked with QUAST and genome completeness was verified with BUSCO using the proteobacteria_odb9 database for comparison. The resulting 62 contigs were annotated online using RAST (rast.nmpdr.org).
DIAMOND was used to assign annotated protein models with a closest blast hit using the uniref90 database, an E-value threshold of 10−10, and a block size of 10. InterProScan (version 5.50-84.0) was used to assign family and domain information to protein models.
Output from both of these searches was used to assign gene ontologies with BLAST2GO. Gene set enrichment analysis was done with BioConductor topGO package, using the default weight algorithm, the TopGO Fisher test, with a p-value threshold of 0.05.
The transposon insertion plasmid, pMiniHimarFRT was delivered to GoB58 by conjugation with E. coli WM3064. E. coli WM3064 transformed with pMiniHimarFRT was grown overnight to saturation in 50 mL LB (10 g L−1 tryptone, 5 g L−1 yeast extract, and 10 g L−1 NaCl) supplemented with 50 μg mL−1 kanamycin (kan) and 90 μM diaminopimelic acid (DAP), rinsed once with 50 mL LB, then re-suspended in 20 mL YPM.
We used a Monte Carlo numerical simulation (collectionmc) to approximate how many insertions would need to occur before a mutant is found representing a knockout of each gene in the genome. Our calculations demonstrated that we would be able to identify mutants in at least 99% of all G. oxydans B58 genes if we generated and selected at least 55,000 mutants (
GoB58 was grown for approximately 24 hours in YPM, then back-diluted to an optical density (OD) of 0.05 in 750 mL YPM and incubated at 30° C. for two doublings until the OD reached 0.2. GoB58 culture was distributed into 13 50 mL conical tubes, to which rinsed and re-suspended WM3064 was added at a ratio of 1:1 by density (approximately 1 mL WM3064 to 50 mL B58). Bacteria were mixed by inversion then spun down at 1900 g for 5 minutes. Supernatant was poured off, and the mixture was resuspended in the remaining liquid (≈0.5 mL), pipetted onto a YPM plate in 5 spots of 0.1 mL, and allowed to dry on the bench under a flame.
Mating plates were incubated at 30° C. for 24 hours. Mating spots were collected by adding 4 mL YPM to a plate, scraping the spots into the liquid, then suspending by pipetting up and down several times. Suspended cells were collected from each plate, and the suspension was plated onto YPM agar with 100 μg mL−1 kanamycin at 100 μL per plate.
After 3 days of incubation at 30° C., colonies were picked into 96-well microplates using a CP7200 colony picking robot (Norgren Systems, Ronceverte WV, USA). Each well contained 150 μL YPM with 100 μg mL−1 kanamycin. For all experiments, GoB58 was grown in polypropylene microplates sealed with a sterile porous membrane (Aeraseal, Catalog Number BS-25, Excel Scientific) and incubated at 30° C. shaking at 800 rpm. Isolated disruption strains were grown for three days to allow nearly all wells to reach saturation. Wells B2 and E7 of each plate were reserved as no-bacteria controls.
In total 18 matings were required to recover and pick a progenitor collection of 49,256 disruption strains into 525 microplates over the course of about two months. Microplates with saturated wells were maintained at 4° C. for up to 3 weeks and incubated an extra night at 30° C. before pooling.
Combinatorial pooling which was done in three batches. The 525 plates were virtually arranged in a 20 by 27 grid, and combinatorial pooling, cryopreservation, pool amplicon library generation, and sequencing were all done as previously described.
Sequencing data for the progenitor collection was processed into a progenitor collection catalog using the KOSUDOKU suite of algorithms. To create a condensed collection, a disruption strain was chosen for each of the 2,733 disrupted genes available in the progenitor collection, first prioritizing close proximity to the translation start, then the total probability of the proposed progenitor collection address. A second strain was chosen from the remaining strains for each gene that had another available. For 50 genes, both disruption strains selected were ambiguously located, and thus a third strain was selected from the remaining collection.
In total, 5,137 disruption strains were isolated and struck-out for single colonies. Many progenitor wells were predicted to have more than one possible strain per well, so for each strain, the number of colonies isolated was two times the predicted number of strains in the progenitor well, up to ten. The condensed collection, which amounted to 17,706 wells, was pooled, sequenced, and validated as previously described. Unknown disruption strains significantly linked to acidification were identified with Sanger sequencing, also as previously described, with the exception of the transposon-specific primers. For the first and second rounds of nested PCR, the transposon-specific primers were (5′-GTATCGCCGCTCCCG-3′ (SEQ ID NO: 309), and (5′-CATCGCCTTCTATCGCCTTC-3′ (SEQ ID NO: 310)), respectively.
Endpoint acidity was measured using the pH indicator thymol blue (TB, Sigma-Aldrich, St. Louis, MO), which changes from red to yellow below a pH of 2.8 (www.sigmaaldrich.com/US/en/product/sial/114545). The lowest pH of biolixiviant generated by GoB58 was 2.3 (Reed2016a), thus TB allows for distinguishing strains that lower the pH below that of the wild type biolixiviant. To generate biolixiviant, the condensed collection was pin replicated into new growth plates containing 100 μL YPM with 100 μg mL−1 kanamycin per well. After two days of growth, an equal volume of 40% w/v glucose was added to the cultures for a final solution of 20% w/v glucose. The amount of glucose needed to lower the pH below 2.3 via the production of gluconic acid was estimated to be 13% w/v, but the higher concentration was used to account for any use of glucose as a carbon source and still maintain an excess amount. Viability tests demonstrated that the bacteria were still viable after two days of culture in such a solution (data not shown).
Bacteria were incubated with glucose for 48 hours to allow acid production to reach completion. Plates were then centrifuged for 3 minutes at 3200 g (top speed) and 90 μL of the biolixiviant supernatant was removed and add to TB at a final concentration of 40 μg mL−1. After 1 minute of vortexing, absorbance was measured for each well at 435 nm and 545 nm on a Synergy 2 plate reader (Biotek Instruments, Winooski, VT, USA). Because of variation in background absorbance from well to well on each plate, absorbance was measured at these two wavelengths, and their ratio was used as a proxy for pH, which correlates linearly within the range of pH for the majority of biolixiviants produced by the collection.
Acidification rate was measured using the pH indicating dye, Bromophenol Blue (BPB). Knockout collection strains were grown for two days. OD was measured at 590 nm for each well, then 5 μL of culture was transferred to a polystyrene assay plate containing 95 μL of 2% w/v glucose and 20 μg mL−1 BPB in deionized water. The initial pH of the culture is just above 5, and within moments of adding culture to glucose with BPB, the color begins to change rapidly. Assay plates were vortexed for one minute after addition of bacterial culture, then immediately transferred to a plate reader where the change in color was tracked by measuring absorbance at 600 nm every minute for 6 minutes, resulting in 7 reads. Mean rate (V) and R-squared were calculated by the Gen5 microplate reader and imager software (Biotek Instruments). A plot of all V relative to OD demonstrated that the two are correlated, thus V was normalized to OD for each well.
Once every well had its assigned data point (A435/A545 for TB, and V/OD for BPB), hits were determined by first identifying outliers for each plate. The interquartile range and upper and lower bounds were calculated in Microsoft Excel considering all wells with cultured disruption strains. Any data point that was more than 1.5 times over or under the upper or lower bound, respectively, was considered an outlier. A disruption strain was considered a hit if over half of the wells for that strain (or 1 of 2) were outliers.
Acidification End Point and Rate Quantification with Colorimetric Dyes
For each assay, knockout strains identified as hits were isolated from the knockout collection into new microplates, along with several blanks per plate, and proxy wild type strains—GoB58 strains with an intergenic transposon insertion that should not affect the acidification phenotype. OD and acidification phenotypes were measured for each proxy WT strain separately to verify that growth and acidification are unaffected in these strains.
Acidification phenotypes for the disruption strains were compared to that of proxy WT with a Student's t-test in Microsoft Excel, two-tailed with equal variance. A Bonferroni correction was used to determine significance to account for the possibility a comparison is significant by chance alone: a phenotype was considered significant if p>0.05/n, where n is the number of comparisons being made (n=120 or n=242 for endpoint acidity comparisons with pWT set A or set B, respectively; n=60 for rate of acidification comparisons with pWT).
The biolixiviant end point pH and acidification rate of each G. oxydans mutant were compared against a proxy wild-type set of mutants for each phenotype. To account for the presence of a kanamycin cassette in the genome, the proxy wild-type set for each phenotype was constructed of several mutants with the transposon inserted in an intergenic region, that had no growth defect, and no apparent change in phenotype.
As the efficiency of the E. coli WM3064 to G. oxydans mating was low, we constructed the G. oxydans progenitor collection in 18 mating batches. As a result of this, the possibility existed that there might be slight variations in the wild type background from batch to batch.
For the acidification rate, we found that these variations did not affect the wild-type behavior across the collection, and a single set of proxy wild-type strains could be used as a comparison with notable disruption strains in the quantification assays. For the end point pH measurement, we found two distinct proxy wild-type behaviors in the condensed collection. For plates 1 to 76; 110 to 130; and 160 to 185, we used proxy wild-type set A, and for plates 77 to 109 and 130 to 159 we used proxy wild-type set B.
For both wild-type sets, we compared ODs after two days of growth, and endpoint acidity using the TB absorbance ration (A435/A545). For wild-type set A, which we used for the BPB quantification assay, we also compared acidification rate of individual proxy WT strains of set A. Comparisons were all made using a linear model, one-way ANOVA, and post-hoc Tukey HSD in R.
Bacteria were grown for 48 hours in tubes containing 4 mL YPM with 100 μg mL−1 kanamycin. One tube was left uninoculated as a no-bacteria control. OD was normalized to 1.9 and diluted in half with 40% glucose for a final 20% solution in 1.5 mL. Five replicates were created for each strain and controls, and all mixtures were randomly distributed across two deep well plates. 750 μL of mixture was transferred from each well to a second set of deep-well plates for bioleaching experiments. All plates were incubated shaking at 900 rpm at room temperature.
After two days, one set of deep-well plates was centrifuged for 10 minutes at 3200 g (top speed), and the pH of the supernatant was measured by insertion of a micro-probe to the same depth in each well.
Four standards were used for meter calibration—pH 1, 2, 4, and 7—and the meter was re-calibrated after every 12 measurements. pH measurements for each disruption strain were compared with those of proxy WT using a Student's t-test in Microsoft Excel, two-tailed, with equal variance. A biolixiviant pH was considered significantly different if p<0.05/n, with n=27.
The second set of deep-well plates was centrifuged for 10 minutes at 3200 g (top speed), and 500 μL of biolixiviant was transferred from each well to a 1.7 mL Eppendorf tube. 20 mg (4% w/v) of retorted phosphor powder was added to each tube for bioleaching. Tubes were shaken horizontally for 36 hours at room temperature, then centrifuged to pellet remaining solids. Supernatant with leached REE was filtered through a 0.45 μm AcroPrep Advance 96-well Filter Plates (Pall Corporation, Show Low, AZ, USA) by centrifuging at 1500×g for 5 minutes.
All samples were diluted 1/200 in 2% trace metal grade nitric acid (Thermo Fisher Scientific) and analyzed by an Agilent 7800 ICP-MS for all REE concentrations using a rare earth element mix standard (Sigma-Aldrich) and a rhodium in-line internal standard (Sigma-Aldrich). Quality control was performed by periodic measurement of standards, blanks, and repeat samples
An additional 1/20 dilution in 2% nitric acid was analyzed for mgdh and pqqc disruption strains, and the no-bacteria control (glucose).
Bioleaching measurements for each disruption strain were compared with those of proxy WT or glucose using a Student's t-test in Microsoft Excel, two-tailed, with equal variance. Total REE extracted was considered significantly different if p<0.05/n, with n=12 for those compared to pWT, and n=2 for those compared to gluco
This application claims priority to U.S. provisional application No. 63/220,475, filed Jul. 10, 2021, and to U.S. provisional application No. 63/152,798, filed Feb. 23, 2021, the disclosures of each of which are incorporated herein by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/US2022/017101 | 2/18/2022 | WO |
Number | Date | Country | |
---|---|---|---|
63220475 | Jul 2021 | US | |
63152798 | Feb 2021 | US |