Andressa M. Venturini
Genome-resolved metagenomics was applied to recover archaeal and bacterial metagenome-assembled genomes (MAGs) from forest and pasture soils of the Brazilian Eastern Amazon, which were subjected to different soil moisture levels (Venturini et al., 2022).
The sampling was carried out in July 2015 in a pristine forest of the Tapajós National Forest (3°17'44.4"S 54°57'46.7"W) and an active cattle pasture (3°18'46.7"S 54°54'34.8"W) in the Brazilian Eastern Amazon. Three soil samples from 0 to 10 cm depth were collected at each site. Samples from each land use were combined, sieved, and subjected to a microcosm experiment under increasing soil moisture levels for 30 days at 25 °C. Moisture treatments were established in triplicates for each land use using 1.5 L glass jars filled with 350 g of soil each. These treatments included the original soil moisture of each site (22% for forest and 24% for pasture) and 100% of moisture at field capacity (FC, 50% for forest and 63% for pasture). Forest and pasture soils from these both moisture treatments were DNA-extracted in duplicate using the PowerLyzer PowerSoil DNA Isolation Kit, totaling 12 DNA samples. DNA samples were checked using 1% agarose gel electrophoresis and a Nanodrop 2000c spectrophotometer. Paired-end shotgun metagenomic sequencing (2 x 150 bp) was performed on an Illumina HiSeq platform at Novogene Co. Ltd. using the NEBNext Ultra II DNA Library Prep Kit for Illumina for library construction. Detailed information about the sites, their soil physicochemical properties, the design of the microcosm experiment, DNA extraction and quantification, and the metagenomic sequencing was previously reported (Lemos et al., 2020; Venturini et al., 2022).
Sample table:
ID | Land use | Treatment |
---|---|---|
A1 | Forest | Original |
A2 | Forest | Original |
A3 | Forest | Original |
A10 | Forest | 100% FC |
A11 | Forest | 100% FC |
A12 | Forest | 100% FC |
A13 | Pasture | Original |
A14 | Pasture | Original |
A15 | Pasture | Original |
A22 | Pasture | 100% FC |
A23 | Pasture | 100% FC |
A24 | Pasture | 100% FC |
1.1. Evaluate paired-end reads using FastQC.
1.2. Trim and filter paired-end reads using Trimmomatic.
Changes in Trimmomatic: adapters (TruSeq3-PE-2)/sliding window minimum quality (20)/head crop length (10)/minimum read length (50).
1.3. Re-evaluate paired-end reads using FastQC.
Created Object Name | Type | Description |
---|---|---|
A1_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A1_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A1_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A2_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A2_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A2_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A3_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A3_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A3_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A10_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A10_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A10_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A11_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A11_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A11_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A12_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A12_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A12_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A13_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A13_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A13_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A14_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A14_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A14_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A15_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A15_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A15_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A22_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A22_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A22_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A23_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A23_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A23_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
Created Object Name | Type | Description |
---|---|---|
A24_trimmed_paired | PairedEndLibrary | Trimmed Reads |
A24_trimmed_unpaired_fwd | SingleEndLibrary | Trimmed Unpaired Forward Reads |
A24_trimmed_unpaired_rev | SingleEndLibrary | Trimmed Unpaired Reverse Reads |
2.1. Merge libraries of each land use to serve as inputs for the assemblers using Merge Reads Libraries.
A1_trimmed_paired |
A2_trimmed_paired |
A3_trimmed_paired |
A10_trimmed_paired |
A11_trimmed_paired |
A12_trimmed_paired |
Created Object Name | Type | Description |
---|---|---|
A1-12_trimmed_paired_merged | PairedEndLibrary | Merge trimmed paired-end libraries (from A1 to A12). |
A13_trimmed_paired |
A14_trimmed_paired |
A15_trimmed_paired |
A22_trimmed_paired |
A23_trimmed_paired |
A24_trimmed_paired |
Created Object Name | Type | Description |
---|---|---|
A13-24_trimmed_paired_merged | PairedEndLibrary | Merge trimmed paired-end libraries (from A13 to A24). |
3.1. Assemble paired-end reads using MEGAHIT.
Changes in MEGAHIT: parameter preset (meta-large).
Created Object Name | Type | Description |
---|---|---|
A1-12_MEGAHIT_large.assembly | Assembly | Assembled contigs |
Created Object Name | Type | Description |
---|---|---|
A13-24_MEGAHIT_large.assembly | Assembly | Assembled contigs |
4.1. Cluster the assembled contigs in bins that correspond to putative population genomes using MetaBAT2.
4.2. Cluster the assembled contigs in bins that correspond to putative population genomes using MaxBin2.
Changes in MaxBin2: marker set (107 and 40)/minimum contig length (2500)/plot markers per contig (Y).
4.3. Cluster the assembled contigs in bins that correspond to putative population genomes using CONCOCT.
4.4. Bin optimization approach using DAS Tool.
A1-12_trimmed_paired_merged |
A13-24_trimmed_paired_merged |
A1-12_trimmed_paired_merged |
Created Object Name | Type | Description |
---|---|---|
A1-12_MaxBin2_107m_bins | BinnedContigs | BinnedContigs from MaxBin2 |
A1-12_trimmed_paired_merged |
Created Object Name | Type | Description |
---|---|---|
A1-12_MaxBin2_40m_bins | BinnedContigs | BinnedContigs from MaxBin2 |
A13-24_trimmed_paired_merged |
Created Object Name | Type | Description |
---|---|---|
A13-24_MaxBin2_107m_bins | BinnedContigs | BinnedContigs from MaxBin2 |
A13-24_trimmed_paired_merged |
Created Object Name | Type | Description |
---|---|---|
A13-24_MaxBin2_40m_bins | BinnedContigs | BinnedContigs from MaxBin2 |
A1-12_trimmed_paired_merged |
A13-24_trimmed_paired_merged |
A1-12_MetaBAT2_bins |
A1-12_MaxBin2_107m_bins |
A1-12_MaxBin2_40m_bins |
A1-12_CONCOCT_bins |
A13-24_MetaBAT2_bins |
A13-24_MaxBin2_107m_bins |
A13-24_MaxBin2_40m_bins |
A13-24_CONCOCT_bins |
5.1. Assign bin completeness and contamination scores using CheckM.
Changes in CheckM: reference tree (full tree).
5.2. Filter bins using CheckM.
Changes in CheckM for ≥ 50% of completeness and ≤ 10% of contamination: reference tree (full tree)/completeness (50)/contamination (10).
Created Object Name | Type | Description |
---|---|---|
A1-12_CheckM_bins50.BinnedContigs | BinnedContigs | HQ BinnedContigs A1-12_CheckM_bins50.BinnedContigs |
Created Object Name | Type | Description |
---|---|---|
A13-24_CheckM_bins50.BinnedContigs | BinnedContigs | HQ BinnedContigs A13-24_CheckM_bins50.BinnedContigs |
6.1. Provide taxonomic assignments for medium-quality bacterial and archaeal bins based on the Genome Taxonomy Database using GTDB-Tk.
7.1. Extract selected bacterial and archaeal bins (high-quality bins plus bins of environmental relevance) as individual assemblies using Extract Bins as Assemblies from BinnedContigs.
Forest: Bin.003/Bin.004/Bin.013.
Pasture: Bin.001/Bin.002/Bin.005/Bin.006/Bin.020/Bin.027/Bin.031/Bin.034/Bin.035.
Bin.003.fasta |
Bin.004.fasta |
Bin.013.fasta |
Created Object Name | Type | Description |
---|---|---|
Forest_H-QBs.AssemblySet | AssemblySet | Assembly set of extracted assemblies |
Bin.003.fasta_forest_assembly | Assembly | Assembly object of extracted contigs |
Bin.004.fasta_forest_assembly | Assembly | Assembly object of extracted contigs |
Bin.013.fasta_forest_assembly | Assembly | Assembly object of extracted contigs |
Bin.001.fasta |
Bin.002.fasta |
Bin.005.fasta |
Bin.006.fasta |
Bin.020.fasta |
Bin.027.fasta |
Bin.031.fasta |
Bin.034.fasta |
Bin.035.fasta |
Created Object Name | Type | Description |
---|---|---|
Pasture_H-QBs.AssemblySet | AssemblySet | Assembly set of extracted assemblies |
Bin.001.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.002.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.005.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.006.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.020.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.027.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.031.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.034.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
Bin.035.fasta_pasture_assembly | Assembly | Assembly object of extracted contigs |
8.1. Calculate the abundance of each selected bin in the merged metagenomes using Align Reads using Bowtie2.
Changes in Bowtie2: preset (very-sensitive).
9.1. Create AssemblySet with the selected bins using Build AssemblySet.
Forest: Bin.003/Bin.004/Bin.013.
Pasture: Bin.001/Bin.002/Bin.005/Bin.006/Bin.020/Bin.027/Bin.031/Bin.034/Bin.035.
9.2. Annotate the selected bins using Annotate and Distill Assemblies with DRAM.
Bin.003.fasta_forest_assembly |
Bin.004.fasta_forest_assembly |
Bin.013.fasta_forest_assembly |
Bin.001.fasta_pasture_assembly |
Bin.002.fasta_pasture_assembly |
Bin.005.fasta_pasture_assembly |
Bin.006.fasta_pasture_assembly |
Bin.020.fasta_pasture_assembly |
Bin.027.fasta_pasture_assembly |
Bin.031.fasta_pasture_assembly |
Bin.034.fasta_pasture_assembly |
Bin.035.fasta_pasture_assembly |
Created Object Name | Type | Description |
---|---|---|
H-QBs.AssemblySet | AssemblySet | KButil_Build_AssemblySet |
Created Object Name | Type | Description |
---|---|---|
Bin.001.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.002.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.003.fasta_forest_assembly_DRAM | Genome | Annotated Genome |
Bin.004.fasta_forest_assembly_DRAM | Genome | Annotated Genome |
Bin.005.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.006.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.013.fasta_forest_assembly_DRAM | Genome | Annotated Genome |
Bin.020.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.027.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.031.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.034.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
Bin.035.fasta_pasture_assembly_DRAM | Genome | Annotated Genome |
H-QBs.GenomeSet | GenomeSet | H-QBs from forest and pasture |
10.1. Download genomes from NCBI public links using Upload File to Staging from Web.
The Patescibacterial genomes were previously recovered by Lemos et al. (2020) using the metagenomic data from the high-moisture (100% FC) pasture soils.
WARW00000000: https://sra-download.ncbi.nlm.nih.gov/traces/wgs01/wgs_aux/WA/RW/WARW01/WARW01.1.gbff.gz.
WARV00000000: https://sra-download.ncbi.nlm.nih.gov/traces/wgs01/wgs_aux/WA/RV/WARV01/WARV01.1.gbff.gz.
10.2. Import GenBank files as genomes using Import GenBank File as Genome from Staging Area.
Changes in Import GenBank File as Genome from Staging Area: genome type (metagenome-assembled genome (MAG))/source of the GenBank file (other).
10.3. Extract bins of Patescibacteria as individual assemblies using Extract Bins as Assemblies from BinnedContigs.
Pasture: Bin.013/Bin.036/Bin.038.
10.4. Compare genomes using Compute ANI with FastANI.
Created Object Name | Type | Description |
---|---|---|
WARW01.1.gbff_genome | Genome | Imported Genome |
Created Object Name | Type | Description |
---|---|---|
WARV01.1.gbff_genome | Genome | Imported Genome |
Bin.036.fasta |
Bin.038.fasta |
Bin.013.fasta |
Created Object Name | Type | Description |
---|---|---|
Pasture_Patescibacteria.AssemblySet | AssemblySet | Assembly set of extracted assemblies |
Bin.036.fasta_pasture_Patescibacteria_assembly | Assembly | Assembly object of extracted contigs |
Bin.038.fasta_pasture_Patescibacteria_assembly | Assembly | Assembly object of extracted contigs |
Bin.013.fasta_pasture_Patescibacteria_assembly | Assembly | Assembly object of extracted contigs |
Bin.013.fasta_pasture_Patescibacteria_assembly |
Bin.036.fasta_pasture_Patescibacteria_assembly |
Bin.038.fasta_pasture_Patescibacteria_assembly |
WARW01.1.gbff_genome_assembly |
WARV01.1.gbff_genome_assembly |
Venturini, A.M., Gontijo, J.B., Mandro, J.A., Paula, F.S., Yoshiura, C.A., da França, A.G., Tsai, S.M., 2022. Genome-resolved metagenomics reveals novel archaeal and bacterial genomes from Amazonian forest and pasture soils. Microbial Genomics 8, 7. https://doi.org/10.1099/mgen.0.000853.
Venturini, A.M., Dias, N.M.S., Gontijo, J.B., Yoshiura, C.A., Paula, F.S., Meyer, K.M., Nakamura, F.M., da França, A.G., Borges, C.D., Barlow, J., Berenguer, E., Nüsslein, K., Rodrigues, J.L.M., Bohannan, B.J.M., Tsai, S.M., 2022. Increased soil moisture intensifies the impacts of forest-to-pasture conversion on methane emissions and methane-cycling communities in the Eastern Amazon. Environmental Research 212, 113139. https://doi.org/10.1016/j.envres.2022.113139
Lemos, L.N., Manoharan, L., Mendes, L.W., Venturini, A.M., Pylro, V.S., Tsai, S.M., 2020. Metagenome assembled‐genomes reveal similar functional profiles of CPR /Patescibacteria phyla in soils. Environmental Microbiology Reports 12, 651–655. https://doi.org/10.1111/1758-2229.12880