
==== Front
Microbiol Resour Announc
Microbiol Resour Announc
mra
Microbiology Resource Announcements
2576-098X
American Society for Microbiology 1752 N St., N.W., Washington, DC

39109829
mra00104-24
10.1128/mra.00104-24
mra.00104-24
Omics Data Sets
applied-and-industrial-microbiologyApplied and Industrial MicrobiologyMetagenomes, metagenome-assembled genomes, and metatranscriptomes from a chlorinated ethene-dechlorinating culture amended with biochar pyrolyzed at different temperatures
https://orcid.org/0000-0003-0132-8318
Dang Hongyu 1 Conceptualization Data curation Formal analysis Methodology Software Visualization Writing – original draft
https://orcid.org/0009-0002-9441-6287
Zhao Weilun 1 Conceptualization Data curation Formal analysis Investigation Methodology Writing – original draft
https://orcid.org/0000-0002-4833-6668
Mattes Timothy E. 1 Conceptualization Data curation Funding acquisition Investigation Project administration Resources Supervision Writing – review and editing tim-mattes@uiowa.edu

1 Department of Civil and Environmental Engineering, University of Iowa , Iowa City, Iowa, USA
Editor Champion Matthew University of Notre Dame , Notre Dame, Indiana, USA

Address correspondence to Timothy E. Mattes, tim-mattes@uiowa.edu
The authors declare no conflict of interest

9 2024
07 8 2024
07 8 2024
13 9 e00104-2409 2 2024
05 7 2024
Copyright © 2024 Dang et al.
2024
Dang et al.
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution 4.0 International license.

ABSTRACT

We investigated the effects of biochar and pyrolysis temperature on a chlorinated ethene-dechlorinating anaerobic consortium. Sequencing of nucleic acids from suspended and biochar-attached cells yielded 9 metagenomes, 122 metagenome-assembled genomes, and 18 metatranscriptomes that provide insights into the structure, function, activity, and interactions of the dehalogenating consortium with biochar.

KEYWORDS

Dehalococcoides
chlorinated ethenes
bioaugmentation
biochar
methanogens
HHS | NIH | National Institute of Environmental Health Sciences (NIEHS) NIH R01ES032671 Mattes Timothy E. HHS | NIH | National Institute of Environmental Health Sciences (NIEHS) P42ES032671 Mattes Timothy E. cover-dateSeptember 2024
==== Body
pmcANNOUNCEMENT

Biochar has been used as a supplement to enhance anaerobic biodegradation of chlorinated ethenes, as it provides habitat, facilitates electron transfer, and stimulates microbial activity and microbe-microbe interactions (1, 2). However, specific biochar properties (e.g., pore size, electron donating/accepting, conductivity, etc.) that are beneficial to microbial ecology, growth, and activity are poorly understood. To fill these research gaps, biochars were synthesized at various temperatures (350°C, 500°C, 700°C, and 900°C) to provide a range of material properties. Each biochar type was added to duplicate microcosms constructed with 100 mL RAMM medium in 160 mL serum bottles and inoculated with an anaerobic tetrachloroethene-dechlorinating consortium (SDC-9) (3–5), together with duplicate live controls without biochar (10 bottles in total).

After 51–82 days of incubation, solid biochar was separated from the liquid culture in bottles by passing through 20 µm UV-sterilized filters (Just the Basics, Woonsocket, RI, USA) and washed with RAMM medium. RNA and DNA were extracted from both liquid and biochar samples with the RNeasy PowerSoil total RNA kit plus DNA Elution kit (Qiagen, Germantown, MD, USA). Residual DNA in RNA extracts was removed with the TURBO DNA-free kit (Thermo Fisher Scientific) and the Direct-zol RNA MiniPrep Plus kit (Zymo Research Corp., Irvine, CA, USA). High-quality RNA (i.e., with clear 16S and 23S rRNA peaks and RNA integrity number > 8), as confirmed with a 2100 Bioanalyzer RNA Pico assay (Agilent Technologies, Santa Clara, CA, USA), was submitted for sequencing.

RNA and DNA were sequenced at the Iowa Institute of Human Genetics (Iowa City, IA, USA). Sheared DNA (average size 550 bp) was used to prepare indexed DNA libraries with the KAPA HyperPrep Kit (Roche Sequencing and Life Science, Indianapolis, IN, USA). RNA libraries were indexed and rRNA-depleted with the stranded total RNA preparation with Ribo-Zero Plus Kit (Illumina Inc., San Diego, CA, USA). Supplemental rRNA probes (149) for dominant taxa in the culture, developed from eight 16S rRNA sequences retrieved from Silva 138.1 (6), were included to improve rRNA depletion and mRNA sequencing efficiency. The DNA and RNA libraries were pooled and sequenced with NovaSeq 6000 on SP (150-bp paired-end) and S1 flow cells (100-bp paired-end), respectively.

Raw sequencing reads were quality controlled by trimming adapters and low-quality sequences (average quality score per base < 15 or sequence length < 36 bases) with Trimmomatic (version 0.39) (7). Trimmed DNA and RNA data sets ranged from 38 to 60 million reads (Table 1), which provided sufficient sequencing depth for metagenomic and metatranscriptomic analysis (8). Reads were assembled into contigs using both individual assembly and co-assembly approaches with Megahit (version 1.2.9) (9). Trimmed short reads were mapped to the contigs with bowtie2 (version 2.2.5) (10). Contigs and their corresponding BAM files were used for binning with Metabat2 (version 2.15) (11). Bin completion and contamination were quantified with CheckM (version 1.2.2) (12) and Anvi’o (version 7.1) (13). Bins with a completion ≥ 80% and contamination ≤ 5% were designated as metagenome-assembled genomes (MAGs). Bins with a completion ≥ 80% and contamination > 5% were further refined with Anvi’o. All MAGs were dereplicated with dRep (version 3.4.2) (14). MAG marker genes, determined by checkM, were aligned with Clustal Omega (version 1.2.3) (15) and used to generate a phylogenetic tree with iqtree (version 2.0.3) (16). MAG taxonomy was classified using GTDB-tk (version 2.1.1; database release R214) and visualized in iTOL (version 5) (17) (Fig. 1). MAGs were annotated using the NCBI Prokaryotic Genome Annotation Pipeline (version 6.0) (18).

TABLE 1 Accession numbers and characteristics before and after trim of metagenomes and metatranscriptome short readsa

Sample name	SRA accession no.	BioSample no.	Pyrolysis temperature (°C)	Replicate	Sample type	Number of reads	
Before trim	After trim-paired	After trim-unpaired R1	After trim-unpaired R2	
Char350-A-liquid_DNA	SRR27243496	SAMN38883632	350	–b	Metagenome	42,164,008	40,942,698	841,272	287,877	
Char350-B-attached_DNA	SRR27243495	SAMN38883633	350	–	Metagenome	53,558,356	52,003,807	1,030,067	386,197	
Char500-A-attached_DNA	SRR27243484	SAMN38883634	500	–	Metagenome	44,813,995	43,558,176	816,976	326,457	
Char500-B-liquid_DNA	SRR27243476	SAMN38883635	500	–	Metagenome	47,600,200	46,070,842	1,020,842	364,793	
Char700-A-attached_DNA	SRR27243475	SAMN38883636	700	–	Metagenome	49,407,025	48,054,961	848,396	379,281	
Char700-B-liquid_DNA	SRR27243474	SAMN38883637	700	–	Metagenome	50,301,937	48,829,200	952,342	385,050	
Char900-A-liquid_DNA	SRR27243473	SAMN38883638	900	–	Metagenome	48,444,322	47,201,060	794,642	338,191	
Char900-B-attached_DNA	SRR27243472	SAMN38883639	900	–	Metagenome	48,296,394	46,942,738	849,928	377,311	
NoChar-A-liquid_DNA	SRR27243471	SAMN38883640	NAc	–	Metagenome	46,581,840	45,244,312	833,038	378,801	
Char350-A-attached_RNA	SRR27243470	SAMN38883641	350	A	Metatranscriptome	57,957,767	56,961,929	652,912	284,891	
Char350-A-liquid_RNA	SRR27243494	SAMN38883642	350	A	Metatranscriptome	48,115,301	47,220,372	698,649	145,630	
Char350-B-attached_RNA	SRR27243493	SAMN38883643	350	B	Metatranscriptome	48,588,053	47,579,440	796,578	143,771	
Char350-B-liquid_RNA	SRR27243492	SAMN38883644	350	B	Metatranscriptome	51,770,161	50,826,589	710,920	165,394	
Char500-A-attached_RNA	SRR27243491	SAMN38883645	500	A	Metatranscriptome	57,882,442	56,775,702	831,448	203,673	
Char500-A-liquid_RNA	SRR27243490	SAMN38883646	500	A	Metatranscriptome	57,447,734	55,879,132	1,329,767	158,883	
Char500-B-attached_RNA	SRR27243489	SAMN38883647	500	B	Metatranscriptome	54,685,915	53,406,752	1,054,073	147,996	
Char500-B-liquid_RNA	SRR27243488	SAMN38883648	500	B	Metatranscriptome	53,516,674	52,587,866	702,192	159,253	
Char700-A-attached_RNA	SRR27243487	SAMN38883649	700	A	Metatranscriptome	39,154,854	38,428,073	512,125	169,707	
Char700-A-liquid_RNA	SRR27243486	SAMN38883650	700	A	Metatranscriptome	60,394,555	59,237,767	918,897	169,161	
Char700-B-attached_RNA	SRR27243485	SAMN38883651	700	B	Metatranscriptome	66,681,063	65,239,485	1,150,431	196,916	
Char700-B-liquid_RNA	SRR27243483	SAMN38883652	700	B	Metatranscriptome	58,947,262	57,787,395	921,054	166,568	
Char900-A- attached_RNA	SRR27243482	SAMN38883653	900	A	Metatranscriptome	58,700,889	57,558,788	843,335	211,378	
Char900-A-liquid_RNA	SRR27243481	SAMN38883654	900	A	Metatranscriptome	66,675,677	65,344,965	1,061,992	181,614	
Char900-B-attached_RNA	SRR27243480	SAMN38883655	900	B	Metatranscriptome	62,066,578	60,842,217	969,950	178,548	
Char900-B-liquid_RNA	SRR27243479	SAMN38883656	900	B	Metatranscriptome	54,717,891	53,611,439	873,365	160,818	
NoChar-A-liquid_RNA	SRR27243478	SAMN38883657	NA	A	Metatranscriptome	51,288,938	50,335,534	743,633	146,313	
NoChar-B-liquid_RNA	SRR27243477	SAMN38883658	NA	B	Metatranscriptome	60,643,367	59,652,887	685,617	235,730	
a In total, 9 (of 18) DNA samples and all 18 RNA samples were sequenced.

b –, no replicate for this sample.

c NA, biochar was not available for this sample.

Fig 1 Phylogenetic tree of 122 MAGs after dereplication, quality statistics determined by CheckM and Anvi’o, genome size, N50, GC content, number of contigs, and rRNA gene recovery (i.e., 16S and 23S, “✓” represents a recovery and “-” represents no recovery). Anvi’o was used to refine bins with completion ≥ 80% and contamination > 5%. Redundant sequences were manually removed based on their read coverage according to bin refinement instructions (https://merenlab.org/2015/05/11/anvi-refine/). The tree was generated using marker genes from each MAG as determined by checkM. The contig containing vcrA (vcrA_contig) was manually added next to the Dehalococcoides MAG on the tree as it has a 99% nucleotide identity with Dehalococcoides mccartyi determined by BLAST (19). Because the vcrA_contig and Dehalococcoides MAG had different tetranucleotide frequencies, they were binned separately. The vcrA_contig GenBank accession number is PP061217.1. The 16S rRNA gene sequences of dominant taxa (marked with “∗”) were collected from the Silva database to design probes for rRNA depletion. Silva taxa (and accession numbers) are Eubacterium (QTVG01000004); Dehalococcides mccartyi, (CP011127); Desulfitobacterium hafniense (AP008230); Desulfovibrio biadhensis (LM999902); Methanocorpusculum aggregans (LMVO01000026); Endomicrobium proavitum (CP009498); and Macellibacteroides fermentans, a(HQ020488).

A total of 122 MAGs belonging to 14 phyla were obtained, including one classified as Dehalococcoides mccartyi (Fig. 1). At least 87% of the metagenomic short reads were mapped to these MAGs, which indicated they represent the microbial composition in SDC-9. These data will advance knowledge of chlorinated ethene-dechlorinating microbial communities and their interactions with biochar.

ACKNOWLEDGMENTS

This work was funded by grants from the National Institute of Environmental Health Sciences (NIEHS; NIH R01ES032671 and P42ES032671). Data were generated at the Genomics Division of the IIHG, which is supported, in part, by the University of Iowa Carver College of Medicine.

We thank both Kevin Knudston and Einat Snir (IIHG) for their contributions to this data set. We also thank Paul Hatzinger for providing the SDC-9 culture.

DATA AVAILABILITY

Raw sequence data and MAGs are available under NCBI BioProject number PRJNA1054096, with metagenome and metatranscriptome SRA accession numbers from SRR27243470 to SRR27243496.
==== Refs
REFERENCES

1 Wu Z, Man Q, Niu H, Lyu H, Song H, Li R, Ren G, Zhu F, Peng C, Li B, Ma X. 2022. Recent advances and trends of trichloroethylene biodegradation: a critical review. Front Microbiol 13 :1053169. doi:10.3389/fmicb.2022.1053169 36620007
2 Bolan S, Hou D, Wang L, Hale L, Egamberdieva D, Tammeorg P, Li R, Wang B, Xu J, Wang T, Sun H, Padhye LP, Wang H, Siddique KHM, Rinklebe J, Kirkham MB, Bolan N. 2023. The potential of biochar as a microbial carrier for agricultural and environmental applications. Sci Total Environ 886 :163968. doi:10.1016/j.scitotenv.2023.163968 37164068
3 Steffan RJ, Vainberg S. 2013. Production and handling of dehalococcoides bioaugmentation cultures, p 89–115. In Stroo HF, Leeson A, Ward CH (ed), Bioaugmentation for groundwater remediation. Springer New York, New York, NY.
4 Vainberg S, Condee CW, Steffan RJ. 2009. Large-scale production of bacterial consortia for remediation of chlorinated solvent-contaminated groundwater. J Ind Microbiol Biotechnol 36 :1189–1197. doi:10.1007/s10295-009-0600-5 19521729
5 Shelton DR, Tiedje JM. 1984. General method for determining anaerobic biodegradation potential. Appl Environ Microbiol 47 :850–857. doi:10.1128/aem.47.4.850-857.1984 6721493
6 Quast C, Pruesse E, Yilmaz P, Gerken J, Schweer T, Yarza P, Peplies J, Glöckner FO. 2013. The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res 41 :D590–6. doi:10.1093/nar/gks1219 23193283
7 Bolger AM, Lohse M, Usadel B. 2014. Trimmomatic: a flexible trimmer for Illumina sequence data. Bioinformatics 30 :2114–2120. doi:10.1093/bioinformatics/btu170 24695404
8 Sims D, Sudbery I, Ilott NE, Heger A, Ponting CP. 2014. Sequencing depth and coverage: key considerations in genomic analyses. Nat Rev Genet 15 :121–132. doi:10.1038/nrg3642 24434847
9 Li D, Liu C-M, Luo R, Sadakane K, Lam T-W. 2015. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics 31 :1674–1676. doi:10.1093/bioinformatics/btv033 25609793
10 Langmead B, Salzberg SL. 2012. Fast gapped-read alignment with Bowtie 2. Nat Methods 9 :357–359. doi:10.1038/nmeth.1923 22388286
11 Kang DD, Li F, Kirton E, Thomas A, Egan R, An H, Wang Z. 2019. MetaBAT 2: an adaptive binning algorithm for robust and efficient genome reconstruction from metagenome assemblies. PeerJ 7 :e7359. doi:10.7717/peerj.7359 31388474
12 Parks DH, Imelfort M, Skennerton CT, Hugenholtz P, Tyson GW. 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 25977477
13 Eren AM, Kiefl E, Shaiber A, Veseli I, Miller SE, Schechter MS, Fink I, Pan JN, Yousef M, Fogarty EC, et al. . 2021. Community-led, integrated, reproducible multi-omics with anvi’o. Nat Microbiol 6 :3–6. doi:10.1038/s41564-020-00834-3 33349678
14 Olm MR, Brown CT, Brooks B, Banfield JF. 2017. dRep: a tool for fast and accurate genomic comparisons that enables improved genome recovery from metagenomes through de-replication. ISME J 11 :2864–2868. doi:10.1038/ismej.2017.126 28742071
15 Sievers F, Wilm A, Dineen D, Gibson TJ, Karplus K, Li W, Lopez R, McWilliam H, Remmert M, Söding J, Thompson JD, Higgins DG. 2011. Fast, scalable generation of high-quality protein multiple sequence alignments using Clustal Omega. Mol Syst Biol 7 :539. doi:10.1038/msb.2011.75 21988835
16 Minh BQ, Schmidt HA, Chernomor O, Schrempf D, Woodhams MD, von Haeseler A, Lanfear R. 2020. IQ-TREE 2: new models and efficient methods for phylogenetic inference in the genomic era. Mol Biol Evol 37 :1530–1534. doi:10.1093/molbev/msaa015 32011700
17 Letunic I, Bork P. 2021. InteractiveTree Of Life (iTOL) v5: an online tool for phylogenetic tree display and annotation. Nucleic Acids Res 49 :W293–W296. doi:10.1093/nar/gkab301 33885785
18 Tatusova T, DiCuccio M, Badretdin A, Chetvernin V, Nawrocki EP, Zaslavsky L, Lomsadze A, Pruitt KD, Borodovsky M, Ostell J. 2016. NCBI prokaryotic genome annotation pipeline. Nucleic Acids Res 44 :6614–6624. doi:10.1093/nar/gkw569 27342282
19 Sayers EW, Beck J, Bolton EE, Bourexis D, Brister JR, Canese K, Comeau DC, Funk K, Kim S, Klimke W, Marchler-Bauer A, Landrum M, Lathrop S, Lu Z, Madden TL, O’Leary N, Phan L, Rangwala SH, Schneider VA, Skripchenko Y, Wang J, Ye J, Trawick BW, Pruitt KD, Sherry ST. 2021. Database resources of the National center for biotechnology information. Nucleic Acids Res 49 :D10–D17. doi:10.1093/nar/gkaa892 33095870
