Single-cell genomics of single soil aggregates: methodological assessment and potential implications with a focus on nitrogen metabolism

Single-cell genomics of single soil aggregates: methodological assessment and potential implications with a focus on nitrogen metabolism | bioRxiv /* */ /* */ <!-- <!-- /*! * yepnope1.5.4 * (c) WTFPL, GPLv2 */ (function(a,b,c){function d(a){return"[object Function]"==o.call(a)}function e(a){return"string"==typeof a}function f(){}function g(a){return!a||"loaded"==a||"complete"==a||"uninitialized"==a}function h(){var a=p.shift();q=1,a?a.t?m(function(){("c"==a.t?B.injectCss:B.injectJs)(a.s,0,a.a,a.x,a.e,1)},0):(a(),h()):q=0}function i(a,c,d,e,f,i,j){function k(b){if(!o&&g(l.readyState)&&(u.r=o=1,!q&&h(),l.onload=l.onreadystatechange=null,b)){"img"!=a&&m(function(){t.removeChild(l)},50);for(var d in y[c])y[c].hasOwnProperty(d)&&y[c][d].onload()}}var j=j||B.errorTimeout,l=b.createElement(a),o=0,r=0,u={t:d,s:c,e:f,a:i,x:j};1===y[c]&&(r=1,y[c]=[]),"object"==a?l.data=c:(l.src=c,l.type=a),l.width=l.height="0",l.onerror=l.onload=l.onreadystatechange=function(){k.call(this,r)},p.splice(e,0,u),"img"!=a&&(r||2===y[c]?(t.insertBefore(l,s?null:n),m(k,j)):y[c].push(l))}function j(a,b,c,d,f){return q=0,b=b||"j",e(a)?i("c"==b?v:u,a,b,this.i++,c,d,f):(p.splice(this.i++,0,a),1==p.length&&h()),this}function k(){var a=B;return a.loader={load:j,i:0},a}var l=b.documentElement,m=a.setTimeout,n=b.getElementsByTagName("script")[0],o={}.toString,p=[],q=0,r="MozAppearance"in l.style,s=r&&!!b.createRange().compareNode,t=s?l:n.parentNode,l=a.opera&&"[object Opera]"==o.call(a.opera),l=!!b.attachEvent&&!l,u=r?"object":l?"script":"img",v=l?"script":u,w=Array.isArray||function(a){return"[object Array]"==o.call(a)},x=[],y={},z={timeout:function(a,b){return b.length&&(a.timeout=b[0]),a}},A,B;B=function(a){function b(a){var a=a.split("!"),b=x.length,c=a.pop(),d=a.length,c={url:c,origUrl:c,prefixes:a},e,f,g;for(f=0;f<d;f++)g=a[f].split("="),(e=z[g.shift()])&&(c=e(c,g));for(f=0;f<b;f++)c=x[f](c);return c}function g(a,e,f,g,h){var i=b(a),j=i.autoCallback;i.url.split(".").pop().split("?").shift(),i.bypass||(e&&(e=d(e)?e:e[a]||e[g]||e[a.split("/").pop().split("?")[0]]),i.instead?i.instead(a,e,f,g,h):(y[i.url]?i.noexec=!0:y[i.url]=1,f.load(i.url,i.forceCSS||!i.forceJS&&"css"==i.url.split(".").pop().split("?").shift()?"c":c,i.noexec,i.attrs,i.timeout),(d(e)||d(j))&&f.load(function(){k(),e&&e(i.origUrl,h,g),j&&j(i.origUrl,h,g),y[i.url]=2})))}function h(a,b){function c(a,c){if(a){if(e(a))c||(j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}),g(a,j,b,0,h);else if(Object(a)===a)for(n in m=function(){var b=0,c;for(c in a)a.hasOwnProperty(c)&&b++;return b}(),a)a.hasOwnProperty(n)&&(!c&&!--m&&(d(j)?j=function(){var a=[].slice.call(arguments);k.apply(this,a),l()}:j[n]=function(a){return function(){var b=[].slice.call(arguments);a&&a.apply(this,b),l()}}(k[n])),g(a[n],j,b,n,h))}else!c&&l()}var h=!!a.test,i=a.load||a.both,j=a.callback||f,k=j,l=a.complete||f,m,n;c(h?a.yep:a.nope,!!i),i&&c(i)}var i,j,l=this.yepnope.loader;if(e(a))g(a,0,l,0);else if(w(a))for(i=0;i (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0];var j=d.createElement(s);var dl=l!='dataLayer'?'&l='+l:'';j.src='//www.googletagmanager.com/gtm.js?id='+i+dl;j.type='text/javascript';j.async=true;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-M677548'); Skip to main content Home About Submit ALERTS / RSS Search for this keyword Advanced Search New Results Single-cell genomics of single soil aggregates: methodological assessment and potential implications with a focus on nitrogen metabolism View ORCID Profile Emi Matsumura , Hiromi Kato , Shintaro Hara , Tsubasa Ohbayashi , Koji Ito , Ryo Shingubara , Tomoya Kawakami , Satoshi Mitsunobu , Tatsuya Saeki , Soichiro Tsuda , Kiwamu Minamisawa , View ORCID Profile Rota Wagai doi: https://doi.org/10.1101/2025.01.09.632267 Emi Matsumura 1 Institute for Agro-Environmental Sciences (NIAES), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Emi Matsumura For correspondence: matsumurae509{at}affrc.go.jp rota{at}affrc.go.jp Hiromi Kato 2 Graduate School of Life Science, Tohoku University , Miyagi, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Shintaro Hara 1 Institute for Agro-Environmental Sciences (NIAES), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Tsubasa Ohbayashi 1 Institute for Agro-Environmental Sciences (NIAES), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Koji Ito 1 Institute for Agro-Environmental Sciences (NIAES), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Ryo Shingubara 3 Research Center for Advanced Analysis (NAAC), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Tomoya Kawakami 1 Institute for Agro-Environmental Sciences (NIAES), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Satoshi Mitsunobu 4 Graduate School of Agriculture, Ehime University , Ehime, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Tatsuya Saeki 5 bitBiome Inc. , Tokyo, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Soichiro Tsuda 5 bitBiome Inc. , Tokyo, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Kiwamu Minamisawa 2 Graduate School of Life Science, Tohoku University , Miyagi, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site Rota Wagai 1 Institute for Agro-Environmental Sciences (NIAES), National Agriculture and Food Research Organization (NARO) , Ibaraki, Japan Find this author on Google Scholar Find this author on PubMed Search for this author on this site ORCID record for Rota Wagai For correspondence: matsumurae509{at}affrc.go.jp rota{at}affrc.go.jp Abstract Full Text Info/History Metrics Supplementary material Preview PDF Abstract Soil particles in plant rooting zone are largely clustered to form complex porous structural unit called aggregates where highly diverse microbes coexist and drive biogeochemical cycling. The complete extraction of microbial cells and DNA from soil is a substantial task as certain microbes exhibit strong adhesion to soil surfaces and/or inhabit deep within aggregates. Yet, the degree of aggregate dispersion and the efficacy of extraction have rarely been examined, and thus adequate cell extraction method from soil remain unclear. We aimed to develop an optimal method of cell extraction for single-cell genomics (SCG) analysis of single soil aggregates by focusing on water-stable macroaggregates (diameter: 5.6-8.2 mm) isolated from a topsoil of cultivated Acrisol. Using the same six individual aggregates, we performed both SCG sequencing and amplicon analysis. While both bead-vortexing and sonication dispersion methods improved the extractability of bacterial cells compared to previous studies, the latter yielded higher number and more diverse microbes compared to the former. The analyses of nitrogen-cycling and exopolysaccharides-related genes suggested that the sonication-assisted extraction led to greater recovery of microbes strongly attached to soil particles and/or inhabited the aggregate subunits that were more physically stable (e.g., aggregate core). Further SCG analysis revealed that all six aggregates held intact microbes having the genes (i.e., potentials) to convert nitrate into all possible nitrogen forms while some low-abundance genes showed inter-aggregate heterogeneity. In addition, all six aggregates studied showed overall similarity in pore characteristics, phylum-level composition, and the microbial functional redundancy. Together, these results suggest that water-stable macroaggregates may act as a functional unit in soil and show potential as a useful experimental unit in soil microbial ecology. Our study also suggest that conventional methods employed for the extraction of cell and DNA may not be optimal. The current findings underscore the necessity to advance extraction methodologies, thereby facilitating a more comprehensive understanding of the microbial diversity and functioning within soil environments. 1 Introduction Soil microorganisms inhabit a highly heterogeneous environment that varies both spatially and temporally, driving biogeochemical cycles essential to terrestrial ecosystems. Although microbial communities are highly diverse ( Gans et al., 2005 ), most biogeochemical transformations are mediated by a limited set of metabolic pathways shared across various taxonomic groups ( Louca et al., 2018 ). Among these cycles, the nitrogen (N)-cycling is tightly controlled by specific groups of soil microbes ( Hayatsu et al., 2008 , Levy-Booth et al., 2014 , Kuypers et al., 2018 ). Nitrous oxide (N₂O) is a potent greenhouse gas and ozone-depleting substance ( Ravishankara et al., 2009 ). The production and transformation of N₂O in the soil is controlled by microbial communities, soil properties, and environmental conditions ( Wan et al., 2024 , Wrage-Mönnig et al., 2018 , Zhuang et al., 2024 ). Denitrification and N₂O reduction basically occur under anaerobic conditions ( Schlüter et al., 2024 ). In the final step of the denitrification process (N₂O → N₂), the primary taxa of nosZ -harboring bacteria responsible for nitrous oxide reduction in soil include Alpha-proteobacteria (e.g., Rhizobiaceae), Beta-proteobacteria (e.g., Burkholderiales), Bacteroidota, Gemmatimonadota, and Delta-proteobacteria ( Hallin et al., 2018 ). However, much remains to be understood. Gaining insights into the ecology of nitrogen-cycle-associated microorganisms in soil is crucial for a comprehensive understanding of N₂O production mechanisms in agricultural soils. Soil structure, in particular aggregates, play a critical role in the creation of heterogeneous microenvironments that are physically stable habitat for microorganisms compared to their surrounding soil matrix ( Oades and Waters, 1991 ). Some aggregates are distinct in surface soils and can be separated by water-sieving as water-stable aggregates ( Oades and Waters, 1991 ). A porous soil-aggregate are formed as a consequence of development of the surface soil layer in pedogenes ( Totsche et al., 2018 ). In a conceptual model of aggregate hierarchy, it has been proposed that the small particles including microaggregate ( 250 µm diameter ( Tisdall and Oades, 1982 ). Conversely, the micro-aggregate exhibits greater physical stability due to the strong binding of iron (hydr) oxide and short-range order minerals via microbial debris and extracellular polymeric substances ( Chenu and Cosentino, 2011 , Six et al., 2004 , Totsche et al., 2018 ). The diverse physicochemical properties of the soil solid phase within the hierarchic structure contribute to form the heterogeneous microenvironments ( Gupta and Germida, 2015 ), and they can be isolated by dispersion at progressively increasing energy level ( Kaiser and Asefaw Berhe, 2014 ). Good soil structure, which relies on effective aggregation, is an integral part of ecosystem functioning ( Rabot et al., 2018 , Costa et al., 2018 , Bronick and Lal, 2005 ). The heterogeneity of soil matrix contributes to the distribution of distinct microbial community compositions, which may vary according to aggregate size ( Sessitsch et al., 2001 , Hemkemeyer et al., 2015 ) or the interior and exterior of aggregates ( Ranjard et al., 1998 , Hattori, 1967 ). Recent research by ( Mitsunobu et al., 2025 ) has shown that the dynamics of microbial N₂O functional groups and O₂/N₂O concentrations are regulated by soil microstructure (e.g., pore networks) using a detailed anatomical approach to single soil aggregates. This study suggests that the physical structure at the scale of individual aggregates influences the habitats of N₂O-reducing bacteria. Specifically, it was observed that nosZ -harboring bacteria were localized deep within the interiors of macroaggregates. However, fundamental questions remain regarding the interactions between soil structural heterogeneity and microbial function. For instance, although aggregates can be considered relatively discrete physical units that maintain natural soil structure, the degree to which microbial community composition varies across individual aggregates is still unknown. Key questions include: What determines the microbial community structure within individual aggregates? What determines the microbial community structure within individual aggregates? What are the taxonomic identities and genetic profiles of bacteria within these aggregates, and where are they specifically located within the aggregate matrix? Addressing these questions will require innovative methods beyond conventional bulk soil genome analysis, which typically assumes soil is a homogeneous environment. Advances in molecular biology have made metagenomics a powerful tool for studying the relationship between functional and taxonomic diversity in microbial communities ( Tringe et al., 2005 , Pan et al., 2014 , Leff et al., 2015 , Fierer et al., 2012b , Fierer et al., 2012a , Fierer et al., 2013 ). The two primary metagenomic methods are amplicon analysis and shotgun sequencing. Shotgun metagenomics, which provides both functional and taxonomic insights, has several advantages and limitations. One major advantage of shotgun metagenomics is that it allows the abundance of each gene to be associated with specific ecological processes, enabling the simultaneous examination of multiple ecosystem functions within a single soil sample ( Allison and Martiny, 2008 ). This multifaceted approach recognizes the significance of ecosystem multifunctionality ( Hector and Bagchi, 2007 ). However, shotgun metagenomics also has some limitations. For example, ribosomal protein genes are often absent from metagenome-assembled genomes (MAGs), complicating functional profiling ( Mise and Iwasaki, 2022 ). Moreover, metagenomic sequencing struggles to achieve strain-resolved genomes ( Arikawa and Hosokawa, 2023 ), despite the substantial functional diversity that can exist between strains ( Lin et al., 2016 , Krause et al., 2006 , Hwangbo et al., 2016 ). High-resolution genomic analyses capable of distinguishing between strains are therefore necessary. Additionally, because shotgun sequencing data often includes extracellular DNA, the presence of exogenous DNA can artificially inflate estimates of microbial diversity and genomic potential ( Carini et al., 2016 , Alteio et al., 2020 ). In recent years, single-cell genomics (SCG) has emerged as a complementary approach to shotgun metagenomics. This method bypasses the binning step inherent in shotgun analysis and enables the analysis of individual, intact cells ( Gawad et al., 2016 ). Most of the single amplified genomes (SAGs) are derived from intact cells ( Ide et al., 2022 ), allowing them to be treated as complete genomic representations of individual organisms. By applying single-cell genomics to specific soil components, we may be able to determine the spatial distribution of bacteria involved in N-cycling, providing insights into both their habitats and taxonomic identities. For the application of single-cell genomics analysis to soil samples, it is essential to first disperse the soil sample and separate microbial cells from soil particles. The soil’s characteristics, especially clay and organic matter content, significantly influence the efficiency of cell extraction (El Mujtar et al., 2022 ). The choice of dispersion method would be crucial in single-cell analysis, as harsher dispersion conditions can increase cell recovery but may compromise cell viability ( Ouyang et al., 2021 , Lee et al., 2021 ). To detach cells from soil particles, both physical dispersion methods (e.g., blending, sonication, vortexing, shaking) and chemical dispersion methods (using ionic or non-ionic buffers) are commonly employed in combination ( Williamson et al., 2011 , Lindahl and Bakken, 1995 , Courtois et al., 2001 , Bakken, 1985 ). For small sample volumes or high-throughput analyses, shaking, vortexing (with or without beads), and sonication are frequently used ( Khalili et al., 2019 , Frossard et al., 2016 ). Sonication (also called ultrasonic treatment), applied at varying energy levels, has been used to partially or fully disperse soil aggregates. This method can reveal factors affecting soil aggregate stability through stepwise disaggregation ( Asano and Wagai, 2014 ). Sonication has demonstrated a higher recovery rate for cell extraction compared to other methods and has proven effective in isolating cells (Uhlířová and Šantrůčková, 2003, Neumann et al., 2013 , Frostegård et al., 1999 ). However, some studies have reported negative results for cell viability following sonication ( Lindahl and Bakken, 1995 , Khalili et al., 2019 , Ellery and Schleyer, 1984 ). Currently, it remains uncertain which extraction technique is most suitable for SCG analysis of soil samples. Many technical challenges remain to understand the linkage between biophysical complexity and genomic diversity within bulk soils. To this end, the current study aimed to characterize soil bacterial community using both metagenomic and SCG approaches by focusing on water-stable macroaggregates, the soil physical subunits that (i) retain the intact micro-structure and microbial habitat, and (ii) can be isolated reproducibly. We intended to fully disperse soil aggregates to maximize cell extraction because some soil microbes are localized in the interior of aggregates (e.g., N 2 O reducers, ( Mitsunobu et al., 2025 )) and/or difficult to detach from soil due presumably to persistent biofilm formation and inhabitation to physically-stable subunits (El Mujtar et al., 2022 , Costa et al., 2018 , Almås et al., 2005 ). The specific objectives of this study were: 1) to establish an extraction method that allows SCG analysis on water-resistant macroaggregates, and 2) to clarify the similarity and differences in taxonomic and functional diversity of bacterial community among the individual aggregates. Using the water-stable macroaggregates from a cropland topsoil (light clay, Acrisol) that were previously used to examine physical control on the distribution of N 2 O-reducing bacteria ( Mitsunobu et al., 2025 ), we optimized the soil dispersion method by comparing cells obtained using two soil dispersion methods (bead-vortexing method and ultrasonic method) using microscopic observation, taking into account the trade-off between soil dispersion and cell damage. Next, we compared two dispersion techniques (bead-vortexing and sonication) and evaluated the recovery of the bacterial community based on quantitative PCR and amplicon-based bacterial composition. Finally, we characterized the N-cycle and other functional aspects of single aggregates based on high-resolution genome analysis and taxonomic and functional diversity evaluation using SCG analysis data. 2 Materials and Methods 2.1 Soil sample and incubation Soil aggregates were collected and incubated in the same manner by Mitsunobu et al. (2025) . The soil aggregates were collected from the 0- to 20cm layer of an Acrisol containing bark compost from the long-term field trial plot at a research farm of the Shizuoka Prefectural Research Institute of Agriculture and Forestry in Iwata, Shizuoka (34°45’12’’N, 137°50’’33E). The average annual temperatures and annual precipitation rates were 16.8 ℃ and 1843 mm/year, respectively (Hamamatsu, 1991-2020). In the basic soil properties, bulk density was 1.2 g/cm 3 , clay was 37.7 %, silt was 25.8 %, sand was 36.5 %, total organic carbon was 3.01 %, C/N ratio was 15, and pH (H 2 O) was 7.0. We transported the soil samples to the laboratory and stored at 4 ℃ until wet sieving. Water-stable soil macroaggregates were isolated by immersing the field-moist soil sample in autoclaved Milli-Q water on top of a 2 mm autoclaved sieve for 10 min, followed by manual vertical movement at the rate of 60 times over a distance of 1.5 cm for 2 min. We selected the aggregates of 5 - 8 mm in diameter (Table S1) and air-dried them for three days under a fume hood covered with a paper towel to prevent the contamination of the dust. We incubated the soil aggregates with artificial soil water which contained low level of nitric acid and glucose to induce N-cycling in a chamber simulating an atmospheric environment. Using the same methods as in Mitsunobu et al. (2025) , we incubated the soil aggregates with the medium simulating soil water based on Yamazaki (2004) with a slight modification. The medium contains 210 μM glucose, 350 μM KNO 3 , and a solution of 410 μM NaH 2 PO 4 ・2H 2 O, 160 μM Na 2 HPO 4 ・ 12H 2 O, 84 μM MgSO 4 ・7H 2 O, 200 μM CaCl 2 ・2H 2 O which adjusted to pH 6.4. All media were autoclaved. We placed the soil aggregate on the autoclaved glass filter paper in a Petri dish and then dropped 1-1.5 mL of the medium per an aggregate onto the paper, which allowed it to soak into the soil aggregate a little at a time without breaking the aggregate. The incubation was carried out in the dark, at 25 °C and 48 h. We changed the air four times during incubation to simulate atmospheric conditions. The weight of incubated aggregates was measured before and after the incubation. Aggregate diameter and volume were measured by the aggregate images picture taken from directly above using Image J tool. The outline of the aggregate was extracted from the image taken from directly above, and the Feret diameter was used as the diameter. We then calculated the aggregate volume from the diameter. Using non-incubated water-stable macroaggregates (n = 5), aggregate pore characteristics were assessed using X-ray micro-computed tomography (µCT). The µCT images of three of the five aggregates were obtained from previous study ( Mitsunobu et al., 2025 ) which gives detail method description. Briefly, the aggregates were scanned at beamline 20B2 of SPring-8 with a resolution of 2.71–2.72 µm per voxel. The µCT images were segmented using grey scale values to distinguish pores from the solid phase. Porosity was calculated as the ratio of pore volume to total aggregate volume. The pore depth distribution was analyzed by dividing each aggregate into three regions: 0– 600 µm from the surface, 600–1200 µm, and 1200 µm to the center. 2.2 Cell extraction method development 2.2.1 Pilot1. Conventional method First, we tested extraction method used in previous studies that applied SCG to soils, with a slight modification ( Hosokawa et al., 2022 , Nishikawa et al., 2022 ). We first mixed 1 g bulk soil samples and a buffer (soil: buffer ratio = 1 g: 3 mL) in a 15 mL tube. The buffer we used in this study was 20 mM potassium phosphate (KPi, pH 7.5) with 0.05% Tween 80. The mixture was dispersed using a mechanical shaking technique at 120 rpm for 10 min, then allowed to stand at room temperature for 5 min. 2 mL of the supernatant was collected in 1.5 mL tube and then centrifuged at 8000×g for 5 min using a benchtop centrifuge (75004251, Thermo Fisher Scientific). These steps (collection and centrifugation) were repeated twice for purification. After the last mixture was centrifuged at 300×g for 5 min, the supernatant was collected. The supernatant was stained with fluorescent dyes, LIVE/DEAD TM BacLightTM Bacterial Viability Kit (Thermo Fisher SCIENTIFIC, USA) to estimate cell numbe. Eight μL of the supernatant was mixed with 1 μL of SYTO9 (50 μM) and 1μL of propidium iodide (PI, 0.5 mg/ μL), incubated for 10 min in the dark at room temperature, and observed under a fluorescence microscopy (Eclipse Ni-U, Nikon, Tokyo, Japan). 2.2.2 Pilot2. Comparison of two aggregate dispersion techniques (bead-vortexing and ultrasonication) We compared the following two dispersion treatments using 0.5 and 0.25 g of bulk soil (dry weight). In the bead-vortexing treatment, the incubated soil was dispersed in KPi buffer (soil: buffer ratio = 1 g :3 mL) with sterile zirconia beads (3 mm in diameter, soil: beads ratio = 1 g: approximately 3 g, equivalent to 16 beads). The dispersion was performed by vortexing for 1 min using a small vortex mixer (N-81, Nissin, Japan). For the 0.25 g treatment, vortexing was conducted for 0.5 min, according to Whiteley et al. (2003) and Griffiths et al. (2003) . For the sonication treatment, the soil plus buffer mixture prepared in the same way as for the bead-vortexing treatment was dispersed using an ultrasonic homogenizer VP-300 (Taitec, Japan) at a low energy intensity of 30 J/ mL (15% intensity cycle of 20 s on, 20 s off), according to Neumann et al. (2013) . The tip of the probe (VP-MT03) was immersed 10 mm into the soil suspension and the sonication was done in an ice bath to reduce cell damage. The mixture after each dispersion treatment was centrifuged at 300×g for 5 min, and the supernatant was collected. The degree of aggregate dispersion was assessed by comparing the particle size distribution of single aggregates after the three dispersion treatments (mechanical shaking, bead-vortexing and sonication) by laser scattering particle size distribution analyzer ((LA-920, HORIBA Corporation, Japan). Three replicates (aggregates) of each treatment were performed. In other words, the volume of particles included in the distribution model on the small size side (< 1.98 μm) was compared. This threshold was identified as the volume-based local minimum between the two modal peaks, located in the 1.73 -1.98 μm particle size class range in the sonication and bead-vortexing dispersion samples. We compared the cell number on LIVE/ DEAD staining in each dispersion treatment as in Pilot1. We observed large amounts of fine soil particles especially after the sonication dispersion (Fig. S1). We thus explored purification techniques to remove the soil particles from the suspension. 2.2.3 Pilot3 : Purification technique after dispersion We firstly tested the density gradient method using Nycodenz ( Kallmeyer et al., 2008 ). The supernatant was collected from 10 g and 0.5 g incubated soil after the sonication treatment. The Nycodenz (Serumwerk Bernburg AG, Bernburg, Germany) with 1.3 g/ml density (60% w/v) was carefully layered under the supernatant in a ratio of supernatant: Nycodenz = 2:7. The sample tube were centrifuged at 10000× g for 20 min, resulting in a sharp band of bacterial cells. Then we carefully collected that to a new tube, then moved to microscopic observation. For 10 g soil, more than 10 8 – 9 cells / soil sample was observed. For 0.5 g soil, however, the sharp band of bacterial cells was not observed after density separation. We thus concluded that Nycodenz purification was not suitable for single aggregates. Next, we tested the sequential washing technique with 0.5 g and 0.25g of bulk soil, following Bakken (1985) . The ca 500 μL and 250 μL of supernatant (S0) was transferred to new tubes after each dispersion treatments described above and centrifuged at 600×g for 15 min, resulting in the first supernatant (S1). After S1 was removed, 500 μL KPi buffer was added to the residue (R1) and subjected to repeated mixing-centrifugation steps, resulting in a series of supernatants (S2 and S3) and RS3. We observed the degree of dispersion and the cell number of each supernatant (S1, S2 and S3 in Figs. 1 and S1). Download figure Open in new tab Fig. 1 Live/dead staining of thesoilsupernatants obtained from repeated washing steps with sonication and dispersion. (A, C, E) 1st supernatant (S1) and (B, D, F) 3rd supernatant (S3) were stained with SYTO9 and PI. SYTO9 (A, B), PI (C, D) and merged (E, F) images are shown. The bar indicates 50μm. S3 has fewer PI-stained particles, indicating that it was purified to cells only. 2.3 Final cell extraction method (Fig. 2) For the bead-vortexing treatment, sterile zirconia beads (3 mm-diameter) in each 1.5 mL tube and vortexed for 1 min. For the sonication treatment, an incubated aggregate was filled up with KPi buffer to 4 mL in 15 mL centrifuge tube. The mixture was dispersed using an ultrasonic homogenizer (30 J/ mL, 4W, cycle of 20 s on, 20 s off) cooling on ice. In case that a soil aggregate was above 0.1 g (dry weight), filled up to 8 mL. The mixture treated with each dispersion technique were centrifuged at 300×g for 5 min. The ca 200 μL in the bead-vortexing treatment or 3.8 mL in the sonication treatment of supernatant (S0) then was collected to 1.5 mL tube. The first supernatant (S1) was removed, and the residue (R1) and new 100μL KPi buffer were subjected to repeated mixing-centrifugation steps with same manner as Pilot3. We proceeded the third supernatant (S3) to the following steps. 2.4 Single-cell genomic analysis We used total of six aggregates (three per dispersion treatment) for single-cell genomics. SAGs of soil microbes were obtained using the SAG-gel method ( Chijiiwa et al., 2020 , Nishikawa et al., 2022 ) from the third supernatant (S3) resulting from two dispersal techniques and purification (2.3). The concentration of cells was determined using LIVE/DEAD BacLight bacterial viability assay (Thermo Fisher Scientific), and the cells were then suspended in DPBS with 1.5% low-gelling-temperature agarose (Sigma-Aldrich, MO, United States) at 1 cell/capsule. After microfluidic single-cell encapsulation in the capsules, the single-cell-encapsulating gel capsules were recovered in the aqueous phase. Then, gel capsules were immersed in Buffer D2 to denature DNA, and multiple displacement amplification was performed for 3 h using the REPLI-g Single Cell Kit (QIAGEN, Hilden, Germany). After multiple displacement amplification, gel capsules were stained with SYBR Green (Thermo Fisher Scientific). FACSMelody cell sorter (Becton, Dickinson and Company, NJ, United States) equipped with a 488 nm excitation laser was used to sort the gel capsules with confirmed DNA amplification into 384-well plates at 1 bead/ well. Following capsule sorting, the 384-well plates were stored at −20°C. For the sequencing analysis, SAG libraries were prepared from the capsule-sorted plates using the xGen DNA Lib prep EZ UNI (Integrated DNA Technologies, Inc., IA, United States). Ligation adaptors were modified to TruSeqCompatible Full-length Adapters UDI (Integrated DNA Technologies, Inc., IA, United States). Each SAG library was sequenced using the DNBSEQ-G400 2 × 150 bp configuration (MGITech CO., Ltd., Beijing, China) with the MGIEasy Universal Library Conversion Kit. Adapter sequences and low-quality reads were eliminated from raw sequence reads of single-cell genome sequences using bbduk.sh (version 38.90) with following options (ktrim=r ref=adapters k=23 mink=11 hdist=1 tpe tbo qtrim=r trimq=10 minlength=40 maxns=1 minavgquality=15).The reads mapped to the masked human genome were eliminated using bbmap.sh (version 38.90) with following options (quickmatch fast untrim minid=0.95 maxindel=3 bwr=0.16 bw=12 minhits=2 path=human_masked_index qtrim=rl trimq=10). The data of human masked index, hg19_main_mask_ribo_animal_allplant_allfungus.fa.gz, was obtained from https://zenodo.org/record/1208052#.X1hBFWf7SdY . These quality-controlled reads of single-cell genomes were assembled de novo into contigs using SPAdes (v3.15.2) ( Bankevich et al., 2012 ) with the following options (--sc --careful --disable-rr --disable-gzip-output –k 21,33,55,77,99,127). Contigs shorter than 200 bp were excluded from the SAG assemblies. CDSs, rRNAs, and tRNAs were predicted from the SAGs using Prokka (v1.14.5) ( Seemann, 2014 ). with the following options (--rawproduct --mincontiglen 200). The quality of the contigs was evaluated using QUAST (v5.0.2) with the default options. The completeness and contamination of SAGs were evaluated using CheckM v1.1.3 ( Parks et al., 2015 ) in lineage workflow with the options (-r --nt) or in taxonomy workflow with the options (--nt domain Bacteria). Taxonomy identification was performed using GTDB-Tk v2.1.0 ( Chaumeil et al., 2022 ) with default options, and GTDB release 207. The quality of each SAG was determined based on Bowers et al. (2017) . SAGs that were < 50% estimated completeness were considered low quality SAGs. SAGs that had ≥ 50% estimated completeness and 90% estimated completeness and <5% estimated contamination, SAGs need to have 23 S, 16 S, and 5 S rRNA genes and at least 18 tRNAs present in the final assembly. The coding sequence (CDS) and amino acid sequence of each SAG was deduced with Bakta v1.4.2 ( Schwengers et al., 2021 ) annotation pipeline, and annotations of each CDS, were determined using eggNOG-mapper v2.1.9 ( Cantalapiedra et al., 2021 ) with eggNOG DB version: 5.0.2 ( Huerta-Cepas et al., 2019 ). For the SAGs of Psudoarthrobacter 87 strains that matched each other in almost full length of 16S rRNA, the gene structures of N-cycling related genes were compared and visualized using GenomeMatcher. The presence or absence of the following genes related to extracellular polysaccharide production ( Cania et al., 2019 ) was then checked based on the KEGG IDs The percentage of SAGs with at least one detected for each exopolysaccharides (EPS) gene was calculated, and the significance difference was tested using Student’s t-test with a threshold of 0.05. Genes encoding enzymes relevant for N-cycling were detected by diamond ( Buchfink et al., 2015 ) blastp search of all CDSs against NCycDB ( Tu et al., 2019 ), a curated database of N-cycling genes, with the threshold of both identity and query coverage greater than or equal to 70%. To evaluate the pathway coverage in the N cycle ( Tu et al., 2017 ), the abundance of functional genes annotated by NCycDB were compared between reaction paths and between aggregates. The sample number of SAGs for each aggregate was 226, 241, 269 and 206, 227, 203 SAGs, for Beads # 1-3 and Sonic # 1-3, respectively. With respect to SAGs with nosZ -like genes, we examined whether nos accessary genes locate adjacent to the nosZ gene by using interProScan 5.65-97.0 ( Jones et al., 2014 ) in EMBL-EBI with default parameters to detect protein signatures of nosDFLY . The structures of nos gene cluster of SAGs were compared with reference genomes of related isolates or MAGs by using GenomeMatcher v3.04 ( Ohtsubo et al., 2008 ). 2.5 Metagenomic analysis We compared the supernatants and the soil residues to assess the difference in cell number and bacteria community composition to assess whether the extracted cells were representative of the soil aggregate. DNA was extracted from a total of 11 aggregate-derived R0 (0.085-0.302 g) and S2 (0.062 – 0.232 g) using an Extrap Soil DNA Kit Plus ver. 2 (BioDynamics Laboratory Inc., Tokyo, Japan) according to the manufacturer’s instructions. These 11 aggregates contained the 4 aggregates analyzed for SCG. 16S rRNA, nosZ -I and - nosZ -II amplicon sequencing and subsequent bioinformatics analysis are the same as the previous studies ( Bamba et al., 2024 , Hara et al., 2024 ). For 16S rRNA amplicon, taxonomy was assigned to ASVs using the SINTAX algorithm ( Edgar, 2016 ) implemented in USEARCH (v11.0.667) against the SILVA database v123 ( Quast et al., 2013 ). The sequence reads of each sample were rarefied to 27,326 reads per sample. For nosZ amplicon, to according to previous study ( Mitsunobu et al., 2025 ), the generated OTU sequences were subjected to Diamond blastx against the database of NosZ amino acid sequences retrieved from the Fungene nosZ repository ( Fish et al., 2013 ). Class level analysis was done based on NCBI taxonomy. The sequence reads of each sample were rarefied to 282 reads in nosZ -I and 856 reads per sample in nosZ -II. 2.6 Quantitative PCR analysis Quantitative PCR (qPCR) assays were conducted using the fluorescent dye SYBR Green (THUNDERBIRD Next SYBR qPCR mix, TOYOBO) by a QuantStudio 3.0 real-time PCR System (Applied Biosystems/Thermo Fisher Scientific). 16S rRNA and genes coding for the two known clades of N 2 O reductase, nosZ -I and nosZ -II, were quantified using the primer pairs Bact1369F/ProK1492R ( Suzuki et al., 2000 ) for the V3-V4 region of the 16S rRNA, nosZ-F/nosZ-R ( Rich et al., 2003 ) for nosZ -I and nosZ-II-F/nosZ-II-R ( Jones et al., 2013 ) for nosZ -II. The PCR reactions for 16S rRNA started with an initial denaturing step of 95 °C for 30 s, followed 40 cycles at 95 °C for 5 s and 60 °C for 30 s. Thermal cycling conditions for nosZ -I were initial denaturation at 95 °C for 30 s, 40 cycles at 95 °C for 5 s, 58 °C for 30 s, 72 °C for 30 s. For nosZ -II, annealing temperature was 54 °C. Melting curve analyses involved a denaturing step at 95 °C for 15 s, annealing at 65 °C for 1 min and melting in 0.1 °C steps up to 95 °C. Standard curves for each assay were generated by serial dilutions of linearized plasmids with cloned fragments of environmental DNA. Amplification efficiencies were 101.7% for 16S rRNA, 110.8% for nosZ -I and 84.2% for nosZ -II. 2.7 Data analysis For alpha diversity analysis, observed ASVs or OTUs, Chao1, Shannon, InvSimpson and Faith’s phylogenetic diversity index were calculated for amplicon and SCG analysis sample. Rarefaction interpolation and extrapolation analysis of taxonomic richness (the observed ASVs) and functional richness (number of KEGG-ID and number of different kinds of N-cycling related genes annotated by NCycDB) was performed using iNEXT v 3.0.0.packages ( Hsieh et al., 2016 ). Furthermore, we investigated the shape of the relationship between species (OTUs based on 16S rRNA) and functional (KEGG-ID) richness to compare redundancy patterns within individual aggregate ( Gamfeldt et al., 2008 , Teichert et al., 2017 ). This relationship is expected to be linear when all species of an assemblage support singular functions, meaning that the loss of any species will produce an important and equivalent decline in functional richness. On the contrary, functionally redundant assemblages will display curvilinear relationships, i.e. saturation trends, as some functional traits are shared by multiple species. For beta diversity analysis, a weighted UniFrac distance matrix was employed in a principal coordinate analysis (PCoA). Subsequently, a permutational multivariate analysis of variance (PERMANOVA) was conducted. The linear discriminant analysis effect size (LEfSe) method determines the features (clades, operational taxonomic units) most likely to explain differences between classes by combining standard tests of statistical significance with additional tests that encode biological consistency and effect relevance ( Segata et al., 2011 ). Soil unit bioindicators at the phylum level were determined using LEfSe analysis with a p-value of 0.05. For molecular phylogenetic analysis focused on Gemmatimonadota, as the most abundant in soil environments and potentially important role in reducing the N 2 O ( Chee-Sanford et al., 2019 , Jones et al., 2013 , Oshiki et al., 2022 ), the alignment of nucleic acid was performed using MAFFT ( Katoh et al., 2019 ). A phylogenetic tree was constructed in MEGA 7 ( Kumar et al., 2016 ) using the maximum likelihood method by best fit model, (Tamura-Nei model for 16S rRNA, General Time Reversible model for nosZ ) selected best-fit model ( Tamura et al., 2011 ) in bootstrap analyses based on 500 replicates. All statistical analyses in community ecological analysis were performed in R software v4.2.1 using phyloseq v1.42.0 ( McMurdie and Holmes, 2013 ), ggplot2 v3.5.0 ( Wickham, 2011 ), microeco v1.4.0 ( Liu et al., 2021 ), vegan v2.6-4 ( Oksanen et al., 2013 ), cowplot v1.1.3 ( Wilke, 2015 ). and ape v5.7-1 ( Paradis and Schliep, 2019 ). 3 Results In this study, we first showed that the combination of dispersion (bead-vortexing or ultrasonic) and sequential washing as the feasible method to extract intact cells from the soil aggregates for SCG analysis. Second, we performed SCG analysis using this method and showed the difference between the two dispersion techniques in their extractability of bacterial cells and DNA with respect to genome quality as well as their taxonomic and functional diversity. Third, by analyzing selected functional genes related to exopolysaccharides (EPS) and N-cycling, we also showed the two dispersion techniques resulted in important differences in the recovery of these genes. 3.1 Development and evaluation of soil bacterial extraction method 3.1.1 Comparison of protocols for soil bacterial extraction We compared previous and new methods of extracting bacterial cell from the soil aggregates for subsequent SCG analysis. Pilot1, using the previous extraction method ( Hosokawa et al., 2022 , Nishikawa et al., 2022 ), resulted in insufficient cell extraction for the studied soil (Fig. S1A). Next, we compared the two dispersion techniques (bead-vortexing and sonication treatment, Pilot2). The post-dispersion soil suspension had a bimodal particle size distribution with greater liberation of smaller-sized particles (ca. 0.5 µm) at the expense of larger-sized particles (ca. 15-90 µm) in the sonication treatment compared to the bead-vortexing treatment. In contrast, the suspension of mechanical shaking showed clearly incomplete dispersion – much greater volume of >50 µm particles and much less volume of < 1 µm particles compared to the other two dispersion techniques (Fig. S2). Specifically, the total volume of particles smaller than the bimodal distribution boundary value of 1.98 µm was 54.5 ± 1.6% after the sonication and 45.4 ± 9.0% after the bead-vortexing treatments (t-test, p = 0.1), which confirmed a greater degree of the aggregate dispersion by the sonication treatment in agreement with our visual inspection. Aggregate pore characterization by X-ray µCT showed that the majority of pore was present in outer zones and the interior zone (deeper than 1.2 mm from the aggregate surface to the center) accounted for only 0-13 % of total pore volumes (Fig. S3) among the five randomly-picked aggregates (diameter range: 4.9-7.1 mm). The pore distribution pattern suggest that outer zones are less physically stable and thus more susceptible for dispersion than the interior zones of these aggregates. Cell counts after the two dispersion treatments was >10 6 (cell/ sample), which was above the minimum counts (10 6 cell/ sample) for SCG analysis (Fig. S1B, C). In the extracted solution, we, however, observed many suspended soil particles (particularly after the sonication treatment) that would interfere with SCG analysis. We thus conducted a sequential washing after the dispersion treatments to remove the suspended particles that would interfere SCG analysis (Pilot3). The supernatants recovered after the first ( Fig. 1A, C, E ) and second washings (data not shown) still contained some soil particles (microaggregates) as well as bacterial cells even after the sonication treatment. In the supernatant after the third washing, we hardly observed any aggregated particles and, based on nucleic acid staining, the number of non-bacterial particles was reduced ( Fig. 1E, F ) while achieving the bacteria count of 10 6-7 per sample (Fig. S1E, F). Based on these results, we adopted the soil extraction method which combines the aggregate dispersion with the washing technique for the aggregate-scale SCG analysis ( Fig. 2 ). Download figure Open in new tab Fig. 2 Illustration of a protocol for extraction of soil bacteria. We used R0 and S2 for qPCR and amplicon analysis and S3 for single cell genomics analysis. 3.1.2 Comparison of bacteria community extracted from soil residues and supernatants We compared the difference in copy number and bacteria community composition between the second supernatants (S2 in Fig. 2 ) and the soil residues after the initial centrifugation (R0 in Fig. 2 ) to assess the extent to which the extracted cells represent the microbiome of the single soil aggregates. The copy number in the supernatant was one to two orders of magnitude lower than that in the soil residues (Table S1). When comparing the two dispersion techniques, the copy number of 16S rRNA, nosZ -I, and nosZ -Ⅱ tended to be higher in the sonication compared to the bead-vortexing treatment while the significant difference was detected only for nosZ -Ⅱ ( Fig. 3C , t-test, p= 0.02). Download figure Open in new tab Fig 3 Copy number and ASV/OTU richness of the second supernatant (S2). The copy number for 16S rRNA(a), nosZ-I(b) and -Ⅱ (c) and the richness (observedASV/OTU)for 16S rRNA(d)and nosZ –I(e). We didn’t perform amplicon analysis for nosZ –Ⅱ regiondueto not enough sample volume. The significant difference were determined by t-test after outlier exclusion. The amplicon analysis of individual aggregates showed that S2/R0 ratio of ASV for 16S rRNA was 0.43 on average (range: 0.31 – 0.57), that for nosZ -I OTU was 0.24 (0.10 – 0.41) (Table S1). On average in each aggregate, 53% (range to 30-70%) and 34% (14 -57%) of ASVs for 16S rRNA and OTUs for nosZ -I in S2 were shared between R0 and S2 of total ASVs or OTUs. The Shannon diversity index of the S2 community (16S rRNA) was higher by in the sonication (5.68 ± 0.29) compared to the breads treatment (4.70 ± 0.17) (Table S1, t-test, p=0.01). The difference in bacterial community composition between the supernatant and the residue as well as between the two dispersion treatments was shown by PCoA of the weighted UniFrac distance based on16S rRNA region (PERMANOVA, Fraction: r 2 = 0.40, F =14.8, p = 0.001, Dispersion: r 2 = 0.065, F =2.4, p = 0.047, Fig. 4A ). The distance between S2 and R0 in the sonication treatment (0.31 ± 0.02) was significantly smaller (t-test, p < 0.001) than that in the bead-vortexing treatment (0.36 ± 0.03). In contrast to 16S rRNA, no significant difference between the two soil fractions and two dispersion treatments was found for nosZ -I and -Ⅱ regions (Fig. S4A, C). We noted that Bradyrhizobium in nosZ -I and Cloacibacterium and Runella (both Bacteroidota) in nosZ -Ⅱ appeared to be higher in the sonication compared to the bead-vortexing treatment (Fig. S4B, D). The characteristic phylum in S2 shown by the LEfSe analysis was Actinobacteriota, Parcubacteria for the bead-vortexing treatment and Armatimonadota, Bacteroidota, Candidatus_Saccharibacteria and Gemmatimonadota for the sonication treatment (Table S2). Download figure Open in new tab Fig 4 Community differences are shown by a PCoAplot of the weighted UniFrac distance matrix (a) and a bar chart at the phylum level (b) across the different fractions and dispersion. R0, residue; S2, second supernatant; Bd,after beads treatment; Sn, after the sonication treatment. The numbers indicate aggregate sample number. Sn2-3 and Bd2-3 Sn3 highlighted in bold with black border, indicate the Sonic and Beads treatment sample in the SCG analysis. 3.2 Profiles of single amplified genomes extracted from single aggregates 3.2.1 SAGs data quality We evaluated basic genomic information including its quality based on the 2,304 SAGs analyzed. Of all SAGs, 2 SAGs were classified as high-quality (HQ), 241 SAGs as medium-quality (MQ), and 2047 SAGs as low-quality (LQ) genomes ( Table 1 ). The remaining 14 SAGs had a contamination rate of 10-34 %, and most of them were Actinobacteriota. The average total length was 3.72, 2.52, and 0.91 Mb, the average N50 was 101, 60, and 25 Kb, and the average number of tRNAs were 39.5, 27.9, and 10.6, for HQ, MQ, LQ, respectively. When comparing the two dispersion techniques, SAGs of HQ and MQ tended to be higher in the bead-vortexing treatment than in the sonication treatment with an average of 15.5 % and 5.4 %, respectively. Similar trends were also observed for the completeness, the number of tRNA and total CDS while no significant difference was shown for N50 ( Table 1 , Fig. S5). When analyzing after excluding the most dominant Actinobacteriota SAGs, all indexes tended to diverge to a greater extent between the aggregates compared to between the dispersion treatments (Fig. S5). View this table: View inline View popup Download powerpoint Table 1. Basic genomic information and quality, which are the evaluation indices for single-cell analysis 3.2.2 Taxonomic and functional diversity Based on phylogenetic annotation based on the GTDB, Actinobacteriota (66.0%, 1034/ 1566 SAGs identified to the phylum level) was the most frequent phyla ( Fig. 5A ). When comparing the two dispersion techniques among SAGs identified by phylum level, the average number of phyla detected was 12.6 in the sonication and 9.3 in the bead-vortexing treatment. The first dominant Actinobacteriota was significantly more frequent after the bead-vortexing treatment (209.0 ± 11.9 SAGs) than after the sonication (135.7 ± 16.8 SAGs) (t-test, p= 0.007). In contrast, the second dominant Proteobacteria and the third dominant Acidobacteriota tended to be higher after the sonication ( Fig. 5A , Fig. S6). The results indicated that the cell extraction with the sonication dispersion allowed the detection of a wider range of phylum compared to that with the bead-vortexing dispersion technique. We calculated Faith’s phylogenetic diversity using the full-length 16S rRNA data. The phylogenetic diversity appeared to be higher after the sonication treatment (175.1 ± 15.4) compared to the bead-vortexing treatment (155.6 ± 14.7) but without statistically significant difference (t test, p =0.2). Similar trend was also observed in the rarefaction curve of ASV ( Fig. 5B ). Download figure Open in new tab Fig. 5 Taxonomic and functional diversity of SAGs. The taxonomic composition of SAGs, annotated based on the GTDB at the phylum level, is presented for all SAGs (A). Rarefaction curves illustrate alpha diversity, including species richness based on 16S rRNA gene sequences (B), functional diversity (C), and functional redundancy (D). The distribution of SAG counts at the genus level for the three most abundant genera within each phylum is detailed in Supplementary Figure S6. We also assessed the functional diversity of the SAGs. When comparing the two dispersion techniques among SAGs, the number of different CDS based on KEGG-ID tended to be higher in the sonication compared to the bead-vortexing treatment ( Fig. 5C ). The comparison of the number of functional genes with species richness showed that total number of functional genes was equally high (4000-5000) with similar plateau shapes for all the aggregates, suggesting high bacterial functional redundancy across the six single aggregates ( Fig. 5D ). However, the inter-aggregate variation became noticeable when focusing on N-cycle genes (Fig. S7A), especially on specific processes (Fig. S7B-D) (see Section 3.3.2 ). 3.2.3 Genetic variation in comparative arrangement of gene cluster within a species One of the advantages of single-cell analysis over MAGs is its capability of detecting the differences among the genes from closely related strains We analyzed the SAGs of Psudoarthrobacter 87 strains that matched each other in almost full length of 16S rRNA, and compared the contigs in which the assimilatory nitrate reductase and nitrite reductase clusters were located at loci of Y7B10_sc-00145, Y7B8_sc-00268, Y7B8_sc-00316, and Y7B8_sc-00152 ( Fig. 6A ). The results showed overall very high homology, except for a loss of about 3500 bp in Y7B10_sc-00145. In Y7B10_sc-00145, a region of about 3500 bp encoding the Radical SAM domain protein gene and the prolipoprotein diacylglyceryl transferase gene is missing, while it was conserved in the other three SAGs. Download figure Open in new tab Fig.6 Comparative arrangement of gene cluster among SAGs. (a) Cluster for nos between SAGs (blue letter) and reference strains (black letter) . (b) Cluster of the assimilatory nitrate reductase and nitrite reductase cluster among four SAGs belong to Pseudoarthrobacter. Colour coding represents % homology based on nucleotide sequences of genes calculated by GenomeMatcher ( Ohtsubo et al., 2008 ). 3.3 Functional gene analysis of single amplified genomes The results above (3.2.2) showed that the dominant phyla were the same between the two dispersion treatments while the sonication treatment obtained higher functional diversity. Next, we carried out more detailed profiling on functional genes that are related to soil aggregation and N-cycling. 3.3.1 Comparison of SAGs harboring EPS related genes We investigated a total of 242 SAGs of middle and high quality for the presence of exopolysaccharides (EPS) related genes as EPS play a major role in cell attachment to soil surface and in aggregate formation ( Wagai et al., 2023 ). A total of 242 SAGs of MQ and HQ were investigated for the presence of EPS related genes (Fig. S8). The proportion of SAGs harboring the polysaccharide export outer membrane protein gene ( wza , KO number K0991) was significantly lower in the bead-vortexing treatment compared with the sonication treatment, 6.1% ± 3.1% and 33.7% ± 8.5%, respectively. In addition, 36 out of 61 SAGs in the sonication treatment and 16 out of 181 SAGs in the bead-vortexing treatment belonged to actinomycetes, but none possessed the wza gene. In contrast, many SAGs (86.2% ± 12.2% in the sonication and 72.2% ± 25.5% in the bead-vortexing treatment) in bacteria other than actinomycetes harbored the wza gene. The same trend was observed for LptBFGC lipopolysaccharide export complex permeasegene ( lptF , KO no. K07091), LptBFGC lipopolysaccharide export complex permease gene ( lptG , KO no. K11720), LptBFGC lipopolysaccharide export complex inner membrane protein gene ( lptC , KO no. K11719), capsular polysaccharide export system permease gene ( kpsE , KO no. K10107) (Fig. S8). 3.3.2 Nitrogen cycling genes We then narrowed down the functional analysis to N cycling genes and assessed how the frequency of major N-cycling genes CDS (annotated by NCycDB) differed among the reaction paths (arrows, Fig. 7 ) and among the six soil aggregates (six-cell boxes with heatmap, Fig. 7 ). Relatively abundant genes (> 20%, bold arrows) were nirBD related to NO 2 - →NH 4 + and narB , nasAB , narGHIJ related to NO 3 - →NO 2 - . In addition, nmo and gdh genes that are related to organic nitrogen degradation, and a GS gene involved in glutamine synthesis were abundant. On the other hand, functional genes related to nitrification, nitrogen fixation and anamox were little detected. Low abundance genes (< 5%, thin arrow) were nosZ , nrfA , norBC and napAB . When comparing among the aggregates, the high abundance genes were typically found in the aggregates isolated after the bead-vortexing dispersion treatment (the upper 3 cells in each box, Fig. 7 ), while the low abundance genes tended to be found in the aggregates under the sonication dispersion treatment. A rarefaction estimates CDS richness tended to indicate the detection of more diverse functional genes by the sonication dispersion except for assimilatory nitrate reduction (ANR) -related genes (Fig. S7). Download figure Open in new tab Fig.7 The potential coverage of pathways in the nitrogen cycle in thesingle aggregate.The sample number of SAGs for each aggregate was 226,241,269 and 206,227,203 SAGs, for Beads# 1-3 and Sonic# 1-3, respectively. We also assessed whether each taxon (phylum) differed in the diversity of functional gene arrangements (gene arrangement/synteny) involved in the N-cycling. Actinobacteriota, the most dominant phylum, harbored multiple genes related to dissimilatory nitrate reduction to ammonium (DNRA), denitrification (nitrate and nitrite reduction), and ANR within single cells except for Thermoleophilia (class) that harbored few N-cycling genes. Proteobacteria, the second dominant phylum also harbored the genes for multiple N-cycling pathways. In contrast with Actinobacteriota, many classes of Proteobacteria harbored genes for denitrification-specific processes, such as nirKS and nosZ -I. The only SAGs with genes involved in nitrification belonged to archaea (Nitrososphaeraceae) Homology search of CDSs of all SAGs by blastp against NcycDB identified NosZ-like CDSs from 15 SAGs ( Table 2 ). Eleven out of these 15 SAGs were taxonomically assigned by Gtdbtk: five SAGs from Acidobacteriota, three from Proteobacteria, two from Gemmatimonadota, one from Verrucomicrobiota. Although the rest of four SAGs were unclassified by Gtdbtk, their NosZ-like CDSs showed homology to inferred NosZ of Bacteroidota, Proteobacteria, and Gemmatimonadota. View this table: View inline View popup Download powerpoint Table 2. Taxonomy of nosZ harboring SAGs The 16S rRNA and nosZ classifications did not match for two SAGs from the bead-vortexing treatment that were not unclassified by Gtdbtk. The comparing the results of the single-cell and amplicon analyses by focusing on Gemmatimonadota, revealed that the sequences of nosZ and of 16S rRNA from SCG were closer to the dominant ASV of Gemmatimonadota identified by the amplicon analyses of both nosZ and 16S rRNA than any of the known isolated strains in Gemmatimonadota (Fig. S9). Furthermore, we examined weather nos accessory genes located in nosZ -flanking region of SAGs from Acidobacteriota, Gemmatimonadota and Verrucomicrobiota, the phyla known to be difficult to culture. Blastp against NCBI and InterProScan analysis detected CDSs with homology to nosDFY , important for formation of the active site of nosZ , and other nos accessory genes in their nosZ -flanking regions. Comparative analysis of nosZ -flanking regions revealed that the gene arrangement of the nos gene cluster was highly conserved among the SAGs and reference genomes ( Fig. 6B ). 4 Discussion 4.1 Optimal extraction methods for SCG and general characteristics of extracted bacteria We were able to obtain genomic information based on SCG analysis from single soil aggregates for the first time by optimizing the extraction method after a series of pilot tests assessing aggregate dispersion techniques (zirconia bead-vortexing and sonication) as well as the purification of dispersed soil suspension. In theory, a high recovery of bacterial cells from soil is achievable by sufficient, if not complete, dispersion of the aggregates and detachment of cells from soil surfaces while minimizing the physical damage to the cells. Our results suggested that the ultrasonic dispersion was more effective than the bead-vortexing for the recovery of bacterial cells and that of particular bacterial groups (especially nosZ harboring bacteria in Fig. 3C ). Below, we first discuss the methodological aspects in terms of the representativeness of the number and composition of bacteria recovered. A fundamental question when applying cell or DNA extraction techniques to soil is the extent to which the bacteria extracted represent the whole microbiome in the soil sample ( Pathan et al., 2021 , Wydro, 2022 ). To a limited extent, we assessed this by comparing the extracted bacterial genes between the soil residue obtained after the soil dispersion/centrifugation (R0) and the second supernatant of the soil suspension (S2) because the third supernatant (S3) was used for the SCG analysis and thus unavailable ( Fig. 2 ). We found that the 40% of ASVs in S2, based on 16S rRNA amplicon analysis, was common between S2 and R0 despite that the number of bacteria in the supernatant was two orders of magnitude lower than the residue (Table S1). These results suggest that the microbial community in the third supernatant used for SCG analysis is likely to be not drastically different from the community present in the bulk sample (i.e., individual macroaggregates). While the assessment of the bacteria extraction efficiency in soil is rather rare in the literature, our results are generally consistent with the previous study which examined the bacteria recovery from pasture soils in New Zealand. Highton et al. (2023) used an extraction method comparable to ours and compared between bulk soils and their water extracts. After the soil dispersion in water using a blender followed by low-speed centrifugation and washing, they used an epifluorescence microscopy based (EFM) quantification method and showed that the number of extracted cells ranged from 3.3 to 9.4% of the bulk soil. In comparison, our current study showed that the copy numbers of S2 relative to that of R0 for 16S rRNA, nosZ -I and -II were 5.2 %, 4.1%, and 9.4%, respectively, using qPCR method (Table S1). To the extent that the copy number of the 16S rRNA gene correlates with the number of cells detected by EFM quantification methods ( Deng et al., 2019 ), the bacteria recoveries in our study were comparable to those of Highton et al. (2023) . Using 16S rRNA amplicon analysis, they further showed that the bacterial community in the extracts and bulk soils shared roughly half of the ASVs and the only major difference at the family level was the greater recovery of Bacillaecae (Firmicutes) in the bulk soil. In our study, the S2/ R0 ratio for species richness of 16S rRNA and nosZ -I was 38.9% and 17.9%, respectively (Table S1), suggesting that roughly a quarter of the bulk soil bacteria community was extracted. A trade-off is likely to be present between dispersion-assisted cell extraction and physical damage to bacterial cells during the dispersion especially for strongly aggregated soils. In fact, we found that greater dispersion by the sonication treatment led to the recovery of more diverse bacteria while their SAGs were, on average, lower in quality ( Table 1 and Fig. S5). Genome quality is strongly affected by DNA damage caused by physical and chemical processes taking place during the sonication treatment. The sonication treatment was able to disrupt the aggregates of ca. 20-100 µm diameter sizes compared to the bead-vortexing treatment (Fig. S2). As soil bacterial community composition can be significantly different among particle/aggregate size fraction ( Biesgen et al., 2020 , Ranjard et al., 1998 ), the greater abundance and diversity of bacteria released by the sonication ( Figs. 3 and 5 ) may be attributable to those associated with the 20-100 µm aggregates. The proportion of the total SAGs higher than the medium quality accounted for on average 16% (bead-vortexing treatment) and 5% (sonication treatment). The sonication-assisted extraction allowed the detection of nos and its surrounding genes even from low-quality category of SAGs ( Fig. 6B ). Furthermore, the sonication treatment led to higher recovery and diversity of DNA ( Figs. 3 and 4A ) and greater diversity of SAGs including their functional genes ( Fig. 5A-C ) compared to the bead-vortexing treatment. Thus, our sonication-assisted extraction method was more effective in characterizing the bacterial community present in the soil aggregates compared to the extraction with the bead-vortexing dispersion. While cavitation effect of ultrasound by sonication has been demonstrated to disrupt cell structure and associated chemical effects including the generation of free radicals can damage DNA ( Li et al., 2018 , Lv et al., 2019 , Miller et al., 1991 , Oyane et al., 2009 ), these negative effect may not be significant. The current extraction method possibly led to higher recovery of soil bacterial community compared to previous soil SCG studies because no other studies used sonication to disperse soil aggregates. At the same time, the sonication-assisted extraction used in the current study likely led to greater degrees of cell damages – medium and high quality SAGs accounted for 5% (sonication treatment) and 16% (bead-vortexing treatment) of total SAGs. In comparison, Nishikawa et al. (2022) performed a mixing (details not reported) plus washing method with high-speed centrifugation for the extraction and showed that 20% of the total SAGs was in medium and high quality SAGs. Aoki et al. (2022) , using a mixing pretreatment (details unknown) followed by density gradient centrifugation (Nycodenz), yielded on average 41% medium- and high-quality SAGs. These previous studies showed higher proportions of medium- and high-quality SAGs compared to the current study, which likely reflect the fact that their soil samples are typically poorly aggregated (beach, desert, mangrove, paddy soils). In our study, the dispersion using the conventional method was insufficient (Figs. S1A and S2). For a further study targeting at the bacteria in aggregates and well-aggregated soils (e.g., high in soil C, clay, and/or metal oxides), it is important to compare different methods such as the sequential washing/sonication method ( Hattori, 1967 , Kohno and Hashimoto, 2008 ) and the sequential dispersion/density gradient centrifugation method ( Nadeem et al., 2013 ). The physical stability of soil aggregate may exert a primary control on cell extraction efficiency given the levels of dispersive energy required for effective aggregate disruption. This hypothesis appears reasonable as cell detachment from soil surfaces is likely to be severely impeded when their habitats (e.g., pores) are present in aggregate interiors and not fully exposed to extracting solution. The aggregates studied here are relatively stable as our Acrisol is relatively high in clay and iron oxides contents ( Mitsunobu et al., 2025 ). In these agricultural topsoils that contain significant amounts of water-stable aggregates, it would be important to sufficiently disperse aggregates for the extraction of cells that reasonably represent the soil bacterial community. 4.2 Greater recovery of “strongly-attached bacteria” by sonication-assisted extraction The two aggregate dispersion techniques used in this study clearly affected the recovery of soil bacteria. The supernatant after the sonication treatment showed the bacterial community composition more similar to the residual soil microbiome compared to those after the bead-vortexing treatment ( Fig. 4 ). The 16S rRNA amplicon sequencing results at the phylum level also showed the superiority of the sonication dispersion treatment compared to the bead-vortexing treatment but with a few exceptions. The Actinobacteriota detected after the bead-vortexing treatment was higher in relative abundance compared to the sonication treatment based on both the amplicon analysis and SAGs ( Fig. 4B and Table S2). This result may be explained by the tendency that Actinobacteriota are easily dispersible ( Ouyang et al., 2021 ) and tolerant to physical stress ( Schimel et al., 2007 ). In contrast, the Bacteroidota were most abundant after the sonication treatment (LEfSe analysis in Table S2). Several genera of Bacteroidota are considered "hard-to-extract" taxa due presumably to their strong binding to the surfaces of soil organic and mineral particles ( Ouyang et al., 2021 ). The following three lines of results suggest that the sonication treatment preferentially released microbial cells that are more strongly attached to soil particle surface and/or present in the interior of aggregates compared to the bead-vortexing treatment. First, the sonication treatment liberated more nosZ -harboring bacteria cells and N 2 O-related functional genes compared to the bead-vortexing treatment. Specifically, a larger number of the cells containing nosZ -I and -II were detected after the sonication treatment ( Fig. 3B, C ) and the S2/ R0 ratio of nosZ -II (12%) was particularly higher after the sonication treatment compared to the bead treatment (1%). In addition, the sonication treatment led to a significantly higher abundance of Bacteroidota and Gemmatimonadota (Table S2), the main phylum possessing nosZ -II ( Hallin et al., 2018 ). Our SCG analysis further showed that Proteobacteria and Acidobacteriota, which detected nosZ in this study, were more diverse with the sonication treatment ( Fig. 5 ). Furthermore, the relative abundance of nosZ key denitrifying functional genes, was higher in the sonication treatment ( Fig. 7 ). Second, the possible localization of the extracted cells in the original aggregates can be inferred by comparing the current results with the previous study where the single aggregates (approximately 6 mm in diameter) isolated from the same bulk soil and incubated in the same way were sliced from the top to the center at ca 300 μm intervals followed by the amplicon sequencing of each slice ( Mitsunobu et al., 2025 ). We found that two nos Z sequences, detected from 2 SAGs (YK46S3_sc-0263 belonging to Gemmatimonadota and YK46S3_sc-0319 belonging to Acidobacteriota), had 100% homology to the nosZ ASVs isolated from the aggregate interior in the previous study ( Mitsunobu et al., 2025 ). Further analysis showed that the 2 SAGs corresponded to 2 ASVs (Gemmatimonadota) and 2 ASVs (Acidobacteriota) in the previous study (Fig. S10). These four ASVs were detected only in the interior region of the aggregates This comparison therefore supports the idea that the sonication dispersion liberated the microbial cells located in the interior of aggregate as well as its exterior. Third, the analysis of EPS-related genes implies the preferential extraction of the cells strongly attached to soil surface and/or reside in physically-stable subunits of the aggregates. The strongly-attached bacteria are difficult to extract for two reasons. First, they are strongly attached to the surface of the soil solid phase using sticky compounds such as extracellular polymeric substances ( Wagai et al., 2023 ), or they selectively reside in subunits of the studied macroaggregate that are difficult to break up (e.g., physically stable microaggregates present within water-stable macroaggregates, ( Six et al., 2004 , Yudina et al., 2022 ). Given that extracellular polymeric substances acts as a major binding agent for soil aggregation ( Chenu and Cosentino, 2011 , Costa et al., 2018 , Oades and Waters, 1991 ), the strongly-attached bacteria are likely associated with physically-stable subunits of the incubated macroaggregates. We therefore hypothesized if the sonication treatment causes more dispersion of physically-stable aggregate subunits, the bacteria possessing EPS-related genes are released more after the sonication compared to the bead-vortexing treatment. The comparison of the functional genes related to the production of EPS ( Cania et al., 2019 ) between the two dispersion treatments revealed that, among the 9 EPS genes detected, 3 of them were significantly higher in the sonication treatment (p < 0.05) while none was higher in the bead-vortexing treatment (Fig. S8). These differences likely resulted from non-Actinobacteriota because (i) 66 % of the bacteria extracted by SCG ( Fig. 5A ) and 83% of the high-quality bacteria belong to Actinobacteriota (Fig. S5), and (ii) the Actinobacteriota was characterized by a lower number of EPS genes (see Result 3.3.1) whereas these numbers were higher for the non-Actinobacteriota that were more extracted from the sonication treatment (see Result. 3.3.1.). There are two possible explanations for the higher abundance of Actinobacteriota under the bead-vortexing treatment. First, Actinobacteriota extracted from the sonication treatment were more easily damaged than those from the bead-vortexing treatment. Second, the relative proportion of Actinobacteriota in the extracted cells was lower in the sonication treatment because the stronger aggregate dispersion by sonication released greater number of the bacteria that were strongly attached to soil particles and/or resided in more stable microaggregates, thereby diluting the Actinobacteriota population. As a further result of supporting the second, gram-positive bacteria having physical tough membrane ( Schimel et al., 2007 ) such as Actinobacteriota, Firmicutes, and Saccharibacteria were more abundant in the bead-vortexing treatment. In contrast, gram-negative bacteria such as Proteobacteria, Acidobacteriota, Dependentiae, Verrucomicrobiota, and Bacteroidota were more abundant in the sonication treatment ( Fig. 5A ). The observed differences in recovered bacteria between the two aggregate dispersion treatments may give important insights into the habitat and ecology of N 2 O reducing bacteria. Nadeem et al. (2013) distinguished two types of habitats for denitrifying microbes based on their attachment strength to soil particles using sequential dispersion/density gradient centrifugation method. For a sandy loam grassland soil, the authors showed that strongly-attached cells produced less N 2 O than loosely-attached cells and the reduced N 2 O production by the strongly-attached cells was at least partially attributable to greater N 2 O reduction (higher activity of N 2 O-reductase, nosZ ). Thus, the strongly-attached cells were presumed to be localized at “inner” habitats such as crevices and cavities of the soil particles and/or present as persistent biofilms that are more resistant to physical dispersion ( Nadeem et al., 2013 ). The sonication treatment in our study dispersed single aggregates more effectively into finer subunits than the bead-vortexing treatment did (Fig. S2) and released more nosZ -harboring bacteria ( Fig. 3 ), implying that nosZ -harboring bacteria may be localized more in the aggregate subunits that are physically stable (e.g., the aggregate interiors where pore was much less abundant, Fig. S3) against the bead-vortexing and/or present in the cell attachment mode that was more susceptible to the fine air bubbles released by the sonication treatment compared to the bead-vortexing. 4.3 Characteristics of functional genes isolated from extracted intact cells We were able to draw three additional inferences that were not obtainable by the metagenomics approach. For the most dominant species among the most extractable phylum, Actinobacteriota, we found inter-individual variation in the arrangement of N-cycling genes (syntenic block) within the same species. Comparison of genes around nasC in 4 SAGs with 100% 16S rRNA sequence match and higher (medium) genomic quality ( Fig. 6A ) showed that one SAG lacked the 3500 bp encoding five genes, while others did not. These results indicate that, even among very closely related SAGs, where the full 16S rRNA length is perfectly matched, genetic variations due to sequence deletions (or insertions) were present. Such deletions would be difficult to detect in MAGs constructed from shotgun metagenomic sequences as the sequences before and after the deletion are almost perfectly matched. In other words, the current study showed evidence of the variation in functional gene arrangement within one species, Actinobacteriota (i.e. even at the population level) for the first time in soil environment. Analysis of nosZ -flanking regions suggested that the intact cell derived SAGs we obtained were most likely from a potentially functional nosZ -harboring bacteria. The genetic arrangement including nos cluster was similar to that of known closely related species that are shown to have nosZ gene expression ( Fig. 6B ). In particular, Gemmatimonas aurantiaca T-27 ( Chee-Sanford et al., 2019 , Park et al., 2017 ) has been studied in detail for N 2 O reduction. Three SAGs belonging to Gemmatimonas (esp. SAG of Y7B10_sc-00065) showed quite similar arrangement in the nos operon, indicating that the obtained SAGs are from a potentially functional nosZ -harboring bacteria. Therefore, not only did they retain the nosZ , but they were a potentially functional nosZ -harboring bacteria because the nosZ was detected as part of the nos operon. The comparison of the results from the single-cell and amplicon analyses by focusing on Gemmatimonadota revealed that the sequences of nosZ and 16S rRNA from SCG were closer to the dominant ASV of Gemmatimonadota identified by the amplicon analyses of both nosZ and 16S rRNA than any of the known isolated strains in Gemmatimonadota (Fig. S9 for nosZ and 16S rRNA). The current study thus showed that single-cell genomics of single soil aggregate has the potential to provide genomic information on microorganisms that are closely related to those dominant in the soil aggregates based on the amplicon analysis. 4.4 Insights from single aggregate analyses The rationale underlying our focus on individual soil aggregates was to achieve a more refined understanding of the taxonomic diversity and functional roles of microbial communities within intact soil microhabitats, without assuming homogeneity of bulk soil samples. We compared species and functional diversity among 6 individual aggregates based on SAG data and 11 aggregates based on amplicon data. Our results showed both uniqueness and similarity in microbial diversity and selected functional gene profiles among the studied individual aggregates. We detected some ASVs unique to individual aggregates. Mean species accumulation plot showed the increase in observed ASVs with the increase in aggregates at least up to eighth aggregates despite that all aggregates were incubated under the same condition prior to the DNA extraction (Fig. S11A). By comparing individual aggregates and the homogenized bulk soil from the same set of soil cores, Simon et al. (2024) showed that five aggregates captured higher diversity than the bulk and the ASVs unique to these aggregates accounted for 20% of total ASVs obtained. Their results indicate that these rare species were heterogeneously distributed at the small scale ( Simon et al., 2024 ). Similarly, we found that 27% of bacterial ASVs was unique to one of the eleven aggregates (Fig. S11B). However, on average, the unique ASVs made up for only 0.5% of the total bacterial abundances (Fig. S11B). The comparison of N-cycling genes also showed some variations among the six aggregates. The gene families involved to anammox, nitrification, nitrogen fixation, denitrification, DNRA, were low abundant. Two of the functional genes were not detected in any aggregates ( hao , hzsA ). These genes showed high inter-aggregate variations. Three factors may be considered to account for the high variability of the undetected and low frequency genes. First, nitrification may be limited due to the low oxygen condition shown in our incubated aggregates ( Mitsunobu et al., 2025 ) and the low concentration of ammonium relative to nitrate in our incubation. Second, nifH , hao , nosZ and hzo are often regarded as low abundant gene families despite their important roles in N cycling ( Kuypers et al., 2018 , Tu et al., 2017 ). Detection of low abundance genes are more difficult and subject to more errors. Third, only a small portion of the entire microbial community in each aggregate was extractable and the strongly-attached cells were likely to be less represented as discussed in 4.1 and 4.2. For example, nosZ -harboring bacteria, which showed increasing trend from the outer surface towards aggregate core based on qPCR ( Mitsunobu et al., 2025 ), were undetectable in one of the three replicates in the bead-vortexing treatment ( Fig. 7 ). In contrast, they were detected from all three replicates of the sonication treatment which led to more effective dispersion of the aggregates (Fig. S2). Furthermore, more diverse N-cycling genes (especially denitrification and DNRA) tended to be extracted after the sonication than the bead-vortexing treatment (Fig. S7), confirming the importance of the dispersion method (Discussion 4.2.). The heatmap of the relative gene abundance for each dispersion method (three boxes aligned side by side in Fig. 7 box) showed some uniqueness among the individual aggregates. The variation among the triplicates was large for nxrAB , nrfA , nosZ (to a limited extent, napAB , nirA , nirKS ) in both dispersion methods, for amoA (to a limited extent, norBC ) in the bead-vortexing treatment only, and for nifH (to a limited extetent, narB , nasAB , nao , nmo ) in the sonication treatment. These inter-aggregate variation in the N-cycling genes may suggest the presence of unique microbial community structure in different aggregates ( Szoboszlay and Tebbe, 2021 ). The current study, nevertheless, revealed overall similarity in the bacterial community structure and their potential functions among the six aggregates, implying that single water-stable macroaggregates may represent microbially and biophysically stable units in soil. While the current study is limited in the number of cells analyzed by SCG (ca 400 per aggregate) and the extraction was conducted after the lab incubation (as opposed to freshly isolated from the soil), the three lines of evidence support this view. First, the inter-aggregate variations in the number of bacteria were relatively small both in the residues (R0) and the supernatant (S2) (quantitative PCR, Fig. 3 ). The variation of species richness in each fraction was also small (R0 in Table S1, S2 in Fig. 3 , S3 in Fig.5B ) although the effect of the two dispersion treatments was present as discussed above ( section 4.2 ). Second, when we sum up the ASVs between the residue and supernatant (R0 and S2), more than 60% of ASVs were common across the eleven aggregates. In addition, the three most dominant phyla were also similar among the ten aggregates in R0 fraction and six aggregates in S3 fraction ( Figs. 4B and 5A ). Third, the bacterial functional diversity and redundancy were high in a similar degree among the six single aggregates as indicated by the comparison between the number of functional genes and SAGs ( Fig. 5C ) and OTUs ( Fig. 5D ). While such high functional redundancy of soil microbial community has been shown at bulk soil scales ( Chen et al., 2021 , Louca et al., 2018 ), our study further showed that the high redundancy is maintained even at the individual aggregate scale ( Fig. 5 ). Fourth, we found rather similar relative abundance of the N-cycling genes across the six aggregates ( Fig. 7 ). Of the seven major pathways we analyzed, ANR, DNRA and organic N metabolism had the highest gene abundance as depicted in the arrow thickness ( Fig. 7 ). Our results at the aggregate level are consistent with the FACE grassland study ( Tu et al., 2017 ) which also showed Actinobacteriota as the main phylum and organic N metabolism and nitrate reduction as the two major N transformation pathways at the bulk soil level, presumably because microbes gain energy and nutrients by these two processes ( Condron et al., 2010 , Moreno-Vivián et al., 1999). Furthermore, we found that the genes associated with denitrification, ANR, DNRA, and N metabolism showed high relative abundance in all six aggregates (see 3x2 box next to each arrow, Fig. 7 ). Notably, all six aggregates contained microbes with genes enabling the conversion of nitrate into all possible nitrogen forms, suggesting that aggregates act as a minimal functional unit for N-cycling in soil. Similarly, Simon et al. (2024) suggested that soil aggregates may act as functional units of soil organic matter turnover based on the correlations among microbial community composition, organic matter content, and its recycling status among the intact aggregates (ca 2 mm in diameter). Similarly high functional redundancy among the aggregates shown in our study further highlights the potential value of soil aggregate as a basic experimental unit to examine microbial diversity-function relationship as aggregates largely maintain the physical, chemical, and microbial condition of in-situ soils. Our results pose a question regarding the efficacy of contemporary cell extraction techniques, such as bead-beating or vortexing, in soils that exhibit strong aggregation due to their suboptimal recovery of cells and DNA. It is, thus, imperative to advance extraction methodologies that will facilitate a more comprehensive understanding of the microbial diversity and functioning within soil environments. Data Availability Statement The datasets, sequence raw data of SAGs and amplicon sequence data for this study can be found in the the DDBJ. All data were registered using BioProject PRJDB19546. Acknowledgement The authors are grateful to S. Sato and M. Bamba (Tohoku University) for amplicon analysis, M. Hayatsu for valuable comments on extraction method development. We also thank Shizuoka Prefectural Research Institute of Agriculture and Forestry for kindly sharing soils from the long-term experimental field. This paper is based on results obtained from a project, JPNP18016, commissioned by the New Energy and Industrial Technology Development Organization (NEDO). References ↵ Allison , S. D. and Martiny , J. B . ( 2008 ) Resistance, resilience, and redundancy in microbial communities . Proceedings of the National Academy of Sciences . 105 , 11512 – 11519 . OpenUrl Abstract / FREE Full Text ↵ Almås , Å. R. , Mulder , J. and Bakken , L. R . ( 2005 ) Trace metal exposure of soil bacteria depends on their position in the soil matrix . Environmental science & technology . 39 , 5927 – 5932 . OpenUrl CrossRef PubMed ↵ Alteio , L. V. , Schulz , F. , Seshadri , R. , Varghese , N. , Rodriguez-Reillo , W. , Ryan , E. , et al. ( 2020 ) Complementary Metagenomic Approaches Improve Reconstruction of Microbial Diversity in a Forest Soil . mSystems . 5 , 10.1128/msystems.00768-19. doi: 10.1128/mSystems.00768-19 . OpenUrl Abstract / FREE Full Text ↵ Aoki , W. , Kogawa , M. , Matsuda , S. , Matsubara , K. , Hirata , S. , Nishikawa , Y. , et al. ( 2022 ) Massively parallel single-cell genomics of microbiomes in rice paddies . Front Microbiol . 13 , 1024640 . doi: 10.3389/fmicb.2022.1024640 . OpenUrl CrossRef PubMed ↵ Arikawa , K. and Hosokawa , M . ( 2023 ) Uncultured prokaryotic genomes in the spotlight: An examination of publicly available data from metagenomics and single-cell genomics . Computational and Structural Biotechnology Journal . ↵ Asano , M. and Wagai , R . ( 2014 ) Evidence of aggregate hierarchy at micro- to submicron scales in an allophanic andisol . Geoderma . 216 , 62 – 74 . doi: 10.1016/j.geoderma.2013.10.005 . OpenUrl CrossRef GeoRef ↵ Bakken , L. R . ( 1985 ) Separation and purification of bacteria from soil . Appl Environ Microbiol . 49 , 1482 – 1487 . doi: 10.1128/aem.49.6.1482-1487.1985 . OpenUrl Abstract / FREE Full Text ↵ Bamba , M. , Akyol , T. Y. , Azuma , Y. , Quilbe , J. , Andersen , S. U. and Sato , S . ( 2024 ) Synergistic effects of plant genotype and soil microbiome on growth in Lotus japonicus . FEMS Microbiol Ecol . 100 , fiae056. doi: 10.1093/femsec/fiae056 . OpenUrl CrossRef ↵ Bankevich , A. , Nurk , S. , Antipov , D. , Gurevich , A. A. , Dvorkin , M. , Kulikov , A. S. , et al. ( 2012 ) SPAdes: a new genome assembly algorithm and its applications to single-cell sequencing . Journal of computational biology . 19 , 455 – 477 . OpenUrl CrossRef PubMed ↵ Biesgen , D. , Frindte , K. , Maarastawi , S. and Knief , C . ( 2020 ) Clay content modulates differences in bacterial community structure in soil aggregates of different size . Geoderma . 376 , 114544 . doi: 10.1016/j.geoderma.2020.114544 . OpenUrl CrossRef ↵ Bowers , R. M. , Kyrpides , N. C. , Stepanauskas , R. , Harmon-Smith , M. , Doud , D. , Reddy , T. B. K. , et al. ( 2017 ) Minimum information about a single amplified genome (MISAG) and a metagenome-assembled genome (MIMAG) of bacteria and archaea . Nature Biotechnology . 35 , 725 – 731 . doi: 10.1038/nbt.3893 . OpenUrl CrossRef PubMed ↵ Bronick , C. J. and Lal , R . ( 2005 ) Soil structure and management: a review . Geoderma . 124 , 3 – 22 . OpenUrl CrossRef GeoRef Web of Science ↵ Buchfink , B. , Xie , C. and Huson , D. H. ( 2015 ) Fast and sensitive protein alignment using DIAMOND . Nature methods . 12 , 59 – 60 . OpenUrl CrossRef PubMed ↵ Cania , B. , Vestergaard , G. , Krauss , M. , Fliessbach , A. , Schloter , M. and Schulz , S . ( 2019 ) A long-term field experiment demonstrates the influence of tillage on the bacterial potential to produce soil structure-stabilizing agents such as exopolysaccharides and lipopolysaccharides . Environmental Microbiome . 14 , 1 – 14 . OpenUrl CrossRef PubMed ↵ Cantalapiedra , C. P. , Hernández-Plaza , A. , Letunic , I. , Bork , P. and Huerta-Cepas , J . ( 2021 ) eggNOG-mapper v2: functional annotation, orthology assignments, and domain prediction at the metagenomic scale . Molecular biology and evolution . 38 , 5825 – 5829 . OpenUrl CrossRef PubMed ↵ Carini , P. , Marsden , P. J. , Leff , J. W. , Morgan , E. E. , Strickland , M. S. and Fierer , N . ( 2016 ) Relic DNA is abundant in soil and obscures estimates of soil microbial diversity . Nature microbiology . 2 , 1 – 6 . OpenUrl ↵ Chaumeil , P.-A. , Mussig , A. J. , Hugenholtz , P. and Parks , D. H . ( 2022 ) GTDB-Tk v2: memory friendly classification with the genome taxonomy database . Bioinformatics . 38 , 5315 – 5316 . OpenUrl CrossRef PubMed ↵ Chee-Sanford , J. , Tian , D. and Sanford , R . ( 2019 ) Consumption of N2O and other N-cycle intermediates by Gemmatimonas aurantiaca strain T-27 . Microbiology . 165 , 1345 – 1354 . OpenUrl CrossRef PubMed ↵ Chen , H. , Ma , K. , Huang , Y. , Yao , Z. and Chu , C . ( 2021 ) Stable Soil Microbial Functional Structure Responding to Biodiversity Loss Based on Metagenomic Evidences . Front Microbiol . 12 , 716764 . doi: 10.3389/fmicb.2021.716764 . OpenUrl CrossRef PubMed ↵ Chenu , C. and Cosentino , D. ( 2011 ). "Microbial Regulation of Soil Structural Dynamics" . In: The architecture and Biology of Soils: Life in Inner Space . Ritz , K. & Young , I. (eds.). ( Wallingford: CABI ), 37 – 70 . ↵ Chijiiwa , R. , Hosokawa , M. , Kogawa , M. , Nishikawa , Y. , Ide , K. , Sakanashi , C. , et al. ( 2020 ) Single-cell genomics of uncultured bacteria reveals dietary fiber responders in the mouse gut microbiota . Microbiome . 8 , 5 . doi: 10.1186/s40168-019-0779-2 . OpenUrl CrossRef PubMed ↵ Condron , L. , Stark , C. , O’Callaghan , M. , Clinton , P. and Huang , Z . ( 2010 ) The role of microbial communities in the formation and decomposition of soil organic matter . Soil microbiology and sustainable crop production . 81 – 118 . ↵ Costa , O. Y. A. , Raaijmakers , J. M. and Kuramae , E. E . ( 2018 ) Microbial Extracellular Polymeric Substances: Ecological Function and Impact on Soil Aggregation . Front Microbiol . 9 , 1636 . doi: 10.3389/fmicb.2018.01636 . OpenUrl CrossRef ↵ Courtois , S. , Frostegård , Å. , Göransson , P. , Depret , G. , Jeannin , P. and Simonet , P . ( 2001 ) Quantification of bacterial subgroups in soil: comparison of DNA extracted directly from soil or from cells previously released by density gradient centrifugation . Environmental microbiology . 3 , 431 – 439 . OpenUrl CrossRef PubMed Web of Science ↵ Deng , L. , Fiskal , A. , Han , X. , Dubois , N. , Bernasconi , S. M. and Lever , M. A . ( 2019 ) Improving the Accuracy of Flow Cytometric Quantification of Microbial Populations in Sediments: Importance of Cell Staining Procedures . Front Microbiol . 10 , 720 . doi: 10.3389/fmicb.2019.00720 . OpenUrl CrossRef PubMed ↵ Edgar , R. C . ( 2016 ) SINTAX: a simple non-Bayesian taxonomy classifier for 16S and ITS sequences . bioRxiv [Preprint]. Available at : doi: 10.1101/074161 (Accessed January 08, 2025). OpenUrl Abstract / FREE Full Text ↵ El Mujtar , V. A. , Chirdo , F. , Lagares , A. , Wall , L. and Tittonell , P. ( 2022 ) Soil bacterial biodiversity characterization by flow cytometry: The bottleneck of cell extraction from soil . Methods in Ecology and Evolution . 13 , 1388 – 1401 . doi: 10.1111/2041-210x.13876 . OpenUrl CrossRef ↵ Ellery , W. and Schleyer , M . ( 1984 ) Comparison of homogenization and ultrasonication as techniques in extracting attached sedimentary bacteria . Marine ecology progress series. Oldendorf . 15 , 247 – 250 . OpenUrl CrossRef ↵ Fierer , N. , Ladau , J. , Clemente , J. C. , Leff , J. W. , Owens , S. M. , Pollard , K. S. , et al. ( 2013 ) Reconstructing the microbial diversity and function of pre-agricultural tallgrass prairie soils in the United States . Science . 342 , 621 – 624 . OpenUrl Abstract / FREE Full Text ↵ Fierer , N. , Lauber , C. L. , Ramirez , K. S. , Zaneveld , J. , Bradford , M. A. and Knight , R . ( 2012a ) Comparative metagenomic, phylogenetic and physiological analyses of soil microbial communities across nitrogen gradients . The ISME journal . 6 , 1007 – 1017 . OpenUrl CrossRef PubMed ↵ Fierer , N. , Leff , J. W. , Adams , B. J. , Nielsen , U. N. , Bates , S. T. , Lauber , C. L. , et al. ( 2012b ) Cross-biome metagenomic analyses of soil microbial communities and their functional attributes . Proceedings of the National Academy of Sciences . 109 , 21390 – 21395 . OpenUrl Abstract / FREE Full Text ↵ Fish , J. A. , Chai , B. , Wang , Q. , Sun , Y. , Brown , C. T. , Tiedje , J. M. , et al. ( 2013 ) FunGene: the functional gene pipeline and repository . Front Microbiol . 4 , 291 . doi: 10.3389/fmicb.2013.00291 . OpenUrl CrossRef PubMed ↵ Frossard , A. , Hammes , F. and Gessner , M. O . ( 2016 ) Flow cytometric assessment of bacterial abundance in soils, sediments and sludge . Frontiers in microbiology . 7 , 903 . OpenUrl PubMed ↵ Frostegård , Å. , Courtois , S. , Ramisse , V. , Clerc , S. , Bernillon , D. , Le Gall , F. , et al. ( 1999 ) Quantification of bias related to the extraction of DNA directly from soils . Applied and environmental microbiology . 65 , 5409 – 5420 . OpenUrl Abstract / FREE Full Text ↵ Gamfeldt , L. , Hillebrand , H. and Jonsson , P. R . ( 2008 ) Multiple functions increase the importance of biodiversity for overall ecosystem functioning . Ecology . 89 , 1223 – 1231 . OpenUrl CrossRef PubMed Web of Science ↵ Gans , J. , Wolinsky , M. and Dunbar , J . ( 2005 ) Computational improvements reveal great bacterial diversity and high metal toxicity in soil . Science . 309 , 1387 – 1390 . OpenUrl Abstract / FREE Full Text ↵ Gawad , C. , Koh , W. and Quake , S. R . ( 2016 ) Single-cell genome sequencing: current state of the science . Nature Reviews Genetics . 17 , 175 – 188 . OpenUrl CrossRef PubMed ↵ Griffiths , R. I. , Whiteley , A. S. , O’Donnell , A. G. and Bailey , M. J . ( 2003 ) Physiological and community responses of established grassland bacterial populations to water stress . Appl Environ Microbiol . 69 , 6961 – 8 . doi: 10.1128/AEM.69.12.6961-6968.2003 . OpenUrl Abstract / FREE Full Text ↵ Gupta , V. V. S. R. and Germida , J. J . ( 2015 ) Soil aggregation: Influence on microbial biomass and implications for biological processes . Soil Biology and Biochemistry . 80 , A3 – A9 . doi: 10.1016/j.soilbio.2014.09.002 . OpenUrl CrossRef ↵ Hallin , S. , Philippot , L. , Löffler , F. E. , Sanford , R. A. and Jones , C. M . ( 2018 ) Genomics and Ecology of Novel N2O-Reducing Microorganisms . Trends in Microbiology . 26 , 43 – 55 . doi: 10.1016/j.tim.2017.07.003 . OpenUrl CrossRef PubMed ↵ Hara , S. , Kakizaki , K. , Bamba , M. , Itakura , M. , Sugawara , M. , Suzuki , A. , et al. ( 2024 ) Does Rhizobial Inoculation Change the Microbial Community in Field Soils? A Comparison with Agricultural Land-use Changes . Microbes and environments . 39 , ME24006. ↵ Hattori , T . ( 1967 ) Microorganisms and soil aggregates as their micro-habitat (in Japanese with English summary) . Bull. Inst. Agr. Res. Tohoku Univ . 18 , 159 – 193 . OpenUrl ↵ Hayatsu , M. , Tago , K. and Saito , M . ( 2008 ) Various players in the nitrogen cycle: diversity and functions of the microorganisms involved in nitrification and denitrification . Soil Science and Plant Nutrition . 54 , 33 – 45 . OpenUrl CrossRef ↵ Hector , A. and Bagchi , R . ( 2007 ) Biodiversity and ecosystem multifunctionality . Nature . 448 , 188 – 190 . OpenUrl CrossRef PubMed Web of Science ↵ Hemkemeyer , M. , Christensen , B. T. , Martens , R. and Tebbe , C. C. ( 2015 ) Soil particle size fractions harbour distinct microbial communities and differ in potential for microbial mineralisation of organic pollutants . Soil Biology and Biochemistry . 90 , 255 – 265 . doi: 10.1016/j.soilbio.2015.08.018 . OpenUrl CrossRef ↵ Highton , M. P. , Bakken , L. R. , Dörsch , P. , Tobias-Hunefeldt , S. , Molstad , L. and Morales , S. E . ( 2023 ) Soil water extract and bacteriome determine N2O emission potential in soils . Biology and Fertility of Soils . 59 , 217 – 232 . doi: 10.1007/s00374-022-01690-5 . OpenUrl CrossRef ↵ Hosokawa , M. , Endoh , T. , Kamata , K. , Arikawa , K. , Nishikawa , Y. , Kogawa , M. , et al. ( 2022 ) Strain-level profiling of viable microbial community by selective single-cell genome sequencing . Sci Rep . 12 , 4443 . doi: 10.1038/s41598-022-08401-y . OpenUrl CrossRef PubMed ↵ Hsieh , T. , Ma , K. and Chao , A . ( 2016 ) iNEXT: an R package for rarefaction and extrapolation of species diversity (H ill numbers) . Methods in ecology and evolution . 7 , 1451 – 1456 . OpenUrl CrossRef ↵ Huerta-Cepas , J. , Szklarczyk , D. , Heller , D. , Hernández-Plaza , A. , Forslund , S. K. , Cook , H. , et al. ( 2019 ) eggNOG 5.0: a hierarchical, functionally and phylogenetically annotated orthology resource based on 5090 organisms and 2502 viruses . Nucleic acids research . 47 , D309 – D314 . OpenUrl CrossRef PubMed ↵ Hwangbo , K. , Um , Y. , Kim , K. Y. , Madhaiyan , M. , Sa , T. M. and Lee , Y . ( 2016 ) Complete genome sequence of Bacillus velezensis CBMB205, a phosphate-solubilizing bacterium isolated from the rhizoplane of rice in the Republic of Korea . Genome Announcements . 4 , doi: 10.1128/genomea.00654-16 . OpenUrl CrossRef ↵ Ide , K. , Nishikawa , Y. , Maruyama , T. , Tsukada , Y. , Kogawa , M. , Takeda , H. , et al. ( 2022 ) Targeted single-cell genomics reveals novel host adaptation strategies of the symbiotic bacteria Endozoicomonas in Acropora tenuis coral . Microbiome . 10 , 220 . doi: 10.1186/s40168-022-01395-9 . OpenUrl CrossRef PubMed ↵ Jones , C. M. , Graf , D. R. H. , Bru , D. , Philippot , L. and Hallin , S . ( 2013 ) The unaccounted yet abundant nitrous oxide-reducing microbial community: A potential nitrous oxide sink . ISME Journal . 7 , 417 – 426 . doi: 10.1038/ismej.2012.125 . OpenUrl CrossRef PubMed ↵ 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 . OpenUrl CrossRef PubMed Web of Science ↵ Kaiser , M. and Asefaw Berhe , A . ( 2014 ) How does sonication affect the mineral and organic constituents of soil aggregates? —A review. Journal of Plant Nutrition and Soil Science . 177 , 479 – 495 . doi: 10.1002/jpln.201300339 . OpenUrl CrossRef ↵ Kallmeyer , J. , Smith , D. C. , Spivack , A. J. and D’Hondt , S . ( 2008 ) New cell extraction procedure applied to deep subsurface sediments . Limnology and Oceanography: Methods . 6 , 236 – 245 . doi: 10.4319/lom.2008.6.236 . OpenUrl CrossRef ↵ Katoh , K. , Rozewicki , J. and Yamada , K. D . ( 2019 ) MAFFT online service: multiple sequence alignment, interactive sequence choice and visualization . Briefings in bioinformatics . 20 , 1160 – 1166 %@ 1467-5463. OpenUrl CrossRef PubMed ↵ Khalili , B. , Weihe , C. , Kimball , S. , Schmidt , K. T. and Martiny , J. B . ( 2019 ) Optimization of a method to quantify soil bacterial abundance by flow cytometry . MSphere . 4 , doi: 10.1128/msphere.00435-19 . OpenUrl CrossRef ↵ Kohno , N. and Hashimoto , T . ( 2008 ) Distribution of bacterial and fungal populations between water-dispersed soil components and water-stable soil aggregates determined using the partial-sterilization and wash-sonic methods . Soil Science and Plant Nutrition . 54 , 725 – 731 . doi: 10.1111/j.1747-0765.2008.00301.x . OpenUrl CrossRef ↵ Krause , A. , Ramakumar , A. , Bartels , D. , Battistoni , F. , Bekel , T. , Boch , J. , et al. ( 2006 ) Complete Biotechnology genome 24: 1385-1391. of the mutualistic, N2-fixing grass endophyte Azoarcus sp. strain BH72 . Nature . ↵ Kumar , S. , Stecher , G. and Tamura , K . ( 2016 ) MEGA7: molecular evolutionary genetics analysis version 7.0 for bigger datasets . Molecular biology and evolution . 33 , 1870 – 1874 . OpenUrl CrossRef PubMed ↵ Kuypers , M. M. , Marchant , H. K. and Kartal , B . ( 2018 ) The microbial nitrogen-cycling network . Nature Reviews Microbiology . 16 , 263 – 276 . OpenUrl CrossRef PubMed ↵ Lee , J. , Kim , H. S. , Jo , H. Y. and Kwon , M. J . ( 2021 ) Revisiting soil bacterial counting methods: Optimal soil storage and pretreatment methods and comparison of culturedependent and -independent methods . PLoS ONE . 16 . doi: 10.1371/journal.pone.0246142 . OpenUrl CrossRef ↵ Leff , J. W. , Jones , S. E. , Prober , S. M. , Barberán , A. , Borer , E. T. , Firn , J. L. , et al. ( 2015 ) Consistent responses of soil microbial communities to elevated nutrient inputs in grasslands across the globe . Proceedings of the National Academy of Sciences . 112 , 10967 – 10972 . OpenUrl Abstract / FREE Full Text ↵ Levy-Booth , D. J. , Prescott , C. E. and Grayston , S. J . ( 2014 ) Microbial functional genes involved in nitrogen fixation, nitrification and denitrification in forest ecosystems . Soil Biology and Biochemistry . 75 , 11 – 25 . OpenUrl CrossRef ↵ Li , J. , Ma , L. , Liao , X. , Liu , D. , Lu , X. , Chen , S. , et al. ( 2018 ) Ultrasound-induced Escherichia coli O157: H7 cell death exhibits physical disruption and biochemical apoptosis . Frontiers in Microbiology . 9 , 2486 . OpenUrl CrossRef PubMed ↵ Lin , H. , Hu , S. , Liu , R. , Chen , P. , Ge , C. , Zhu , B. , et al. ( 2016 ) Genome sequence of Pseudomonas koreensis CRS05-R5, an antagonistic bacterium isolated from rice paddy field . Frontiers in Microbiology . 7 , 1756 . OpenUrl PubMed ↵ Lindahl , V. and Bakken , L. R . ( 1995 ) Evaluation of methods for extraction of bacteria from soil . FEMS Microbiology Ecology . 16 , 135 – 142 . OpenUrl CrossRef Web of Science ↵ Liu , C. , Cui , Y. , Li , X. and Yao , M . ( 2021 ) microeco: an R package for data mining in microbial community ecology . FEMS microbiology ecology . 97 , fiaa255. ↵ Louca , S. , Polz , M. F. , Mazel , F. , Albright , M. B. , Huber , J. A. , O’Connor , M. I. , et al. ( 2018 ) Function and functional redundancy in microbial systems . Nature ecology & evolution . 2 , 936 – 943 . OpenUrl CrossRef PubMed ↵ Lv , R. , Mingming , Z. , Chen , W. , Wang , D. , Zhou , J. , Ding , T. , et al. ( 2019 ) A hurdle approach of acidic electrolyzed water simultaneous with ultrasound to inactivate Bacillus cereus . Journal of Food Processing and Preservation . 43 , e14148 . OpenUrl ↵ McMurdie , P. J. and Holmes , S . ( 2013 ) phyloseq: an R package for reproducible interactive analysis and graphics of microbiome census data . PloS one . 8 , e61217 . OpenUrl CrossRef PubMed ↵ Miller , D. , Thomas , R. and Frazier , M . ( 1991 ) Ultrasonic cavitation indirectly induces single strand breaks in DNA of viable cells in vitro by the action of residual hydrogen peroxide . Ultrasound in medicine & biology . 17 , 729 – 735 . OpenUrl CrossRef PubMed ↵ Mise , K. and Iwasaki , W . ( 2022 ) Unexpected absence of ribosomal protein genes from metagenome-assembled genomes . ISME communications . 2 , 118 . OpenUrl CrossRef PubMed ↵ Mitsunobu , S. , Wagai , R. , Shimada , H. , Kato , H. , Ito , K. , Sato , S. , et al. ( 2025 ) First microscale data on depth profiles of microbial N₂O reduction, O2 availability, and pore networks inside contrasting single soil aggregates . Soil Biology and Biochemistry . 202 , 109684 . doi: 10.1016/j.soilbio.2024.109684 . OpenUrl CrossRef Moreno-Vivián , C. , Cabello, P. n., Martínez-Luque, M., Blasco, R. and Castillo, F. ( 1999 ) Prokaryotic nitrate reduction: molecular properties and functional distinction among bacterial nitrate reductases . Journal of bacteriology . 181 , 6573 – 6584 . OpenUrl FREE Full Text ↵ Nadeem , S. , Almås , Å. R. , Dörsch , P. and Bakken , L. R . ( 2013 ) Sequential extraction of denitrifying organisms from soils; strongly attached cells produce less N2O than loosely attached cells . Soil Biology and Biochemistry . 67 , 62 – 69 . doi: 10.1016/j.soilbio.2013.08.007 . OpenUrl CrossRef Web of Science ↵ Neumann , D. , Heuer , A. , Hemkemeyer , M. , Martens , R. and Tebbe , C. C . ( 2013 ) Response of microbial communities to long-term fertilization depends on their microhabitat . FEMS Microbiology Ecology . 86 , 71 – 84 . OpenUrl CrossRef PubMed Web of Science ↵ Nishikawa , Y. , Kogawa , M. , Hosokawa , M. , Wagatsuma , R. , Mineta , K. , Takahashi , K. , et al. ( 2022 ) Validation of the application of gel beads-based single-cell genome sequencing platform to soil and seawater . ISME Communications . 2 . doi: 10.1038/s43705-022-00179-4 . OpenUrl CrossRef ↵ Oades , J. and Waters , A . ( 1991 ) Aggregate hierarchy in soils . Soil Research . 29 , 815 – 828 . OpenUrl CrossRef ↵ Ohtsubo , Y. , Ikeda-Ohtsubo , W. , Nagata , Y. and Tsuda , M . ( 2008 ) GenomeMatcher: a graphical user interface for DNA sequence comparison . BMC Bioinformatics . 9 , 376 . doi: 10.1186/1471-2105-9-376 . OpenUrl CrossRef PubMed ↵ Oksanen , J. , Blanchet , F. G. , Kindt , R. , Legendre , P. , Minchin , P. R. , O’hara , R. , et al. ( 2013 ) Package ‘vegan’ . Community ecology package, version . 2 , 1 – 295 . OpenUrl ↵ Oshiki , M. , Toyama , Y. , Suenaga , T. , Terada , A. , Kasahara , Y. , Yamaguchi , T. , et al. ( 2022 ) N(2)O Reduction by Gemmatimonas aurantiaca and Potential Involvement of Gemmatimonadetes Bacteria in N(2)O Reduction in Agricultural Soils . Microbes Environ . 37 . doi: 10.1264/jsme2.ME21090 . OpenUrl CrossRef ↵ Ouyang , Y. , Chen , D. , Fu , Y. , Shi , W. , Provin , T. , Han , A. , et al. ( 2021 ) Direct cell extraction from fresh and stored soil samples: Impact on microbial viability and community compositions . Soil Biology and Biochemistry . 155 . doi: 10.1016/j.soilbio.2021.108178 . OpenUrl CrossRef ↵ Oyane , I. , Takeda , T. , Oda , Y. , Sakata , T. , Furuta , M. , Okitsu , K. , et al. ( 2009 ) Comparison between the effects of ultrasound and γ-rays on the inactivation of Saccharomyces cerevisiae: Analyses of cell membrane permeability and DNA or RNA synthesis by flow cytometry . Ultrasonics sonochemistry . 16 , 532 – 536 . OpenUrl CrossRef PubMed ↵ Pan , Y. , Cassman , N. , de Hollander , M. , Mendes , L. W. , Korevaar , H. , Geerts , R. H. , et al. ( 2014 ) Impact of long-term N, P, K, and NPK fertilization on the composition and potential functions of the bacterial community in grassland soil . FEMS microbiology ecology . 90 , 195 – 205 . OpenUrl CrossRef PubMed ↵ Paradis , E. and Schliep , K . ( 2019 ) ape 5.0: an environment for modern phylogenetics and evolutionary analyses in R . Bioinformatics . 35 , 526 – 528 . OpenUrl CrossRef PubMed ↵ Park , D. , Kim , H. and Yoon , S . ( 2017 ) Nitrous oxide reduction by an obligate aerobic bacterium, Gemmatimonas aurantiaca strain T-27 . Applied and Environmental Microbiology . 83 , e00502 – 17 . OpenUrl CrossRef PubMed ↵ 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 research . 25 , 1043 – 1055 . OpenUrl Abstract / FREE Full Text ↵ Pathan , S. I. , Arfaioli , P. , Taskin , E. , Ceccherini , M. T. , Puglisi , E. and Pietramellara , G . ( 2021 ) The extracellular DNA can baffle the assessment of soil bacterial community, but the effect varies with microscale spatial distribution . FEMS Microbiology Letters . 368 . doi: 10.1093/femsle/fnab074 . OpenUrl CrossRef PubMed ↵ Quast , C. , Pruesse , E. , Yilmaz , P. , Gerken , J. , Schweer , T. , Yarza , P. , et al. ( 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 . OpenUrl CrossRef PubMed Web of Science ↵ Rabot , E. , Wiesmeier , M. , Schlüter , S. and Vogel , H.-J . ( 2018 ) Soil structure as an indicator of soil functions: A review . Geoderma . 314 , 122 – 137 . OpenUrl CrossRef ↵ Ranjard , L. , Poly , F. , Combrisson , J. , Richaume , A. and Nazaret , S . ( 1998 ) A single procedure to recover DNA from the surface or inside aggregates and in various size fractions of soil suitable for PCR-based assays of bacterial communities . European journal of soil biology . 34 , 89 – 97 . OpenUrl CrossRef ↵ Ravishankara , A. , Daniel , J. S. and Portmann , R. W . ( 2009 ) Nitrous oxide (N2O): the dominant ozone-depleting substance emitted in the 21st century . science . 326 , 123 – 125 . OpenUrl Abstract / FREE Full Text ↵ Rich , J. , Heichen , R. , Bottomley , P. , Cromack Jr , K. and Myrold , D . ( 2003 ) Community composition and functioning of denitrifying bacteria from adjacent meadow and forest soils . Applied and environmental microbiology . 69 , 5974 – 5982 . OpenUrl Abstract / FREE Full Text ↵ Schimel , J. , Balser , T. C. and Wallenstein , M . ( 2007 ) Microbial stress-response physiology and its implications for ecosystem function . Ecology . 88 , 1386 – 94 . doi: 10.1890/06-0219 . OpenUrl CrossRef PubMed Web of Science ↵ Schlüter , S. , Lucas , M. , Grosz , B. , Ippisch , O. , Zawallich , J. , He , H. , et al. ( 2024 ) The anaerobic soil volume as a controlling factor of denitrification: a review . Biology and Fertility of Soils . doi: 10.1007/s00374-024-01819-8 . OpenUrl CrossRef ↵ Schwengers , O. , Jelonek , L. , Dieckmann , M. A. , Beyvers , S. , Blom , J. and Goesmann , A . ( 2021 ) Bakta: rapid and standardized annotation of bacterial genomes via alignment-free sequence identification . Microbial genomics . 7 , 000685 . OpenUrl PubMed ↵ Seemann , T . ( 2014 ) Prokka: rapid prokaryotic genome annotation . Bioinformatics . 30 , 2068 – 2069 . OpenUrl CrossRef PubMed Web of Science ↵ Segata , N. , Izard , J. , Waldron , L. , Gevers , D. , Miropolsky , L. , Garrett , W. S. , et al. ( 2011 ) Metagenomic biomarker discovery and explanation . Genome biology . 12 , 1 – 18 . OpenUrl CrossRef PubMed ↵ Sessitsch , A. , Weilharter , A. , Gerzabek , M. H. , Kirchmann , H. and Kandeler , E . ( 2001 ) Microbial population structures in soil particle size fractions of a long-term fertilizer field experiment . Applied and environmental microbiology . 67 , 4215 – 4224 . OpenUrl Abstract / FREE Full Text ↵ Simon , E. , Guseva , K. , Darcy , S. , Alteio , L. , Pjevac , P. , Schmidt , H. , et al. ( 2024 ) Distinct microbial communities are linked to organic matter properties in millimetre-sized soil aggregates . ISME J . 18 . doi: 10.1093/ismejo/wrae156 . OpenUrl CrossRef ↵ Six , J. , Bossuyt , H. , Degryze , S. and Denef , K . ( 2004 ) A history of research on the link between (micro)aggregates, soil biota, and soil organic matter dynamics . Soil and Tillage Research . 79 , 7 – 31 . doi: 10.1016/j.still.2004.03.008 . OpenUrl CrossRef Web of Science ↵ Suzuki , M. T. , Taylor , L. T. and DeLong , E. F . ( 2000 ) Quantitative analysis of small-subunit rRNA genes in mixed microbial populations via 5′-nuclease assays . Applied and environmental microbiology . 66 , 4605 – 4614 . OpenUrl Abstract / FREE Full Text ↵ Szoboszlay , M. and Tebbe , C. C . ( 2021 ) Hidden heterogeneity and co-occurrence networks of soil prokaryotic communities revealed at the scale of individual soil aggregates . Microbiologyopen . 10 , e1144 . doi: 10.1002/mbo3.1144 . OpenUrl CrossRef PubMed ↵ Tamura , K. , Peterson , D. , Peterson , N. , Stecher , G. , Nei , M. and Kumar , S . ( 2011 ) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods . Molecular biology and evolution . 28 , 2731 – 2739 . OpenUrl CrossRef PubMed Web of Science ↵ Teichert , N. , Lepage , M. , Sagouis , A. , Borja , A. , Chust , G. , Ferreira , M. T. , et al. ( 2017 ) Functional redundancy and sensitivity of fish assemblages in European rivers, lakes and estuarine ecosystems . Sci Rep . 7 , 17611 . doi: 10.1038/s41598-017-17975-x . OpenUrl CrossRef PubMed ↵ Tisdall , J. M. and Oades , J. M . ( 1982 ) Organic matter and water-stable aggregates in soils . Journal of soil science . 33 , 141 – 163 . OpenUrl CrossRef GeoRef Web of Science ↵ Totsche , K. U. , Amelung , W. , Gerzabek , M. H. , Guggenberger , G. , Klumpp , E. , Knief , C. , et al. ( 2018 ) Microaggregates in soils . Journal of Plant Nutrition and Soil Science . 181 , 104 – 136 . OpenUrl CrossRef ↵ Tringe , S. G. , Von Mering , C. , Kobayashi , A. , Salamov , A. A. , Chen , K. , Chang , H. W. , et al. ( 2005 ) Comparative metagenomics of microbial communities . Science . 308 , 554 – 557 . OpenUrl Abstract / FREE Full Text ↵ Tu , Q. , He , Z. , Wu , L. , Xue , K. , Xie , G. , Chain , P. , et al. ( 2017 ) Metagenomic reconstruction of nitrogen cycling pathways in a CO2-enriched grassland ecosystem . Soil Biology and Biochemistry . 106 , 99 – 108 . doi: 10.1016/j.soilbio.2016.12.017 . OpenUrl CrossRef ↵ Tu , Q. , Lin , L. , Cheng , L. , Deng , Y. and He , Z . ( 2019 ) NCycDB: a curated integrative database for fast and accurate metagenomic profiling of nitrogen cycling genes . Bioinformatics . 35 , 1040 – 1048 . doi: 10.1093/bioinformatics/bty741 . OpenUrl CrossRef PubMed Uhlířová, E. and Šantrůčková, H . ( 2003 ) Growth rate of bacteria is affected by soil texture and extraction procedure . Soil Biology and Biochemistry . 35 , 217 – 224 . OpenUrl CrossRef ↵ Wagai , R. , Yang , P.-T. and Kaiser , K. ( 2023 ). "Interfacial reactions of microorganisms with minerals and organic matter" . In: Encyclopedia of Soils in the Environment ( Second Edition ). Goss , M. J. & Oliver , M. (eds.). ( Oxford: Academic Press ), 458 – 469 . ↵ Wan , X. , Yang , X. , Zhang , Y. , Hu , P. , Duan , P. , Li , D. , et al. ( 2024 ) Lithology-driven soil properties control of N2O production by ammonia oxidizers in subtropical forest soils . Geoderma . 446 , 116910 . OpenUrl CrossRef ↵ Whiteley , A. S. , Griffiths , R. I. and Bailey , M. J . ( 2003 ) Analysis of the microbial functional diversity within water-stressed soil communities by flow cytometric analysis and CTC+ cell sorting . Journal of Microbiological Methods . 54 , 257 – 267 . doi: 10.1016/S0167-7012(03)00066-6 . OpenUrl CrossRef PubMed ↵ Wickham , H . ( 2011 ) ggplot2 . Wiley interdisciplinary reviews: computational statistics . 3 , 180 – 185 . OpenUrl CrossRef ↵ Wilke , C. O. ( 2015 ) cowplot: Streamlined plot theme and plot annotations for’ggplot2’ . CRAN: Contributed Packages . ↵ Williamson , K. E. , Kan , J. , Polson , S. W. and Williamson , S. J . ( 2011 ) Optimizing the indirect extraction of prokaryotic DNA from soils . Soil Biology and Biochemistry . 43 , 736 – 748 . OpenUrl CrossRef ↵ Wrage-Mönnig , N. , Horn , M. A. , Well , R. , Müller , C. , Velthof , G. and Oenema , O . ( 2018 ) The role of nitrifier denitrification in the production of nitrous oxide revisited . Soil Biology and Biochemistry . 123 , A3 – A16 . OpenUrl CrossRef ↵ Wydro , U . ( 2022 ) Soil microbiome study based on DNA extraction: a review . Water . 14 , 3999 . OpenUrl CrossRef Yamazaki , S. , Yamaguchi , T. , Araki , N. , Sumino , H. , Ohashi , A. & Harada , H . ( 2004 ) Measurements of oxygen and nitrogen profiles in attached biofilm on river rocks by the use of microelectrodes (in Japanese with English abstract) . Doboku-Gakkai Ronbun Shu . 769 , 89 – 94 . OpenUrl ↵ Yudina , A. V. , Klyueva , V. V. , Romanenko , K. A. and Fomin , D. S . ( 2022 ) Micro-within macro: How micro-aggregation shapes the soil pore space and water-stability . Geoderma . 415 . doi: 10.1016/j.geoderma.2022.115771 . OpenUrl CrossRef ↵ Zhuang , S. , Ding , J. , Lin , W. , Zheng , Q. , Kou , X. , Li , Q. , et al. ( 2024 ) Transient anoxic conditions boost N2O emissions by stimulating denitrification capacity and decreasing N2O reduction ratio in soils with different carbon substrates . Soil Biology and Biochemistry . 192 , 109351 . OpenUrl CrossRef View the discussion thread. Back to top Previous Next Posted January 11, 2025. Download PDF Supplementary Material Email Thank you for your interest in spreading the word about bioRxiv. NOTE: Your email address is requested solely to identify you as the sender of this article. Your Email * Your Name * Send To * Enter multiple addresses on separate lines or separate them with commas. You are going to email the following Single-cell genomics of single soil aggregates: methodological assessment and potential implications with a focus on nitrogen metabolism Message Subject (Your Name) has forwarded a page to you from bioRxiv Message Body (Your Name) thought you would like to see this page from the bioRxiv website. Your Personal Message CAPTCHA This question is for testing whether or not you are a human visitor and to prevent automated spam submissions. Share Single-cell genomics of single soil aggregates: methodological assessment and potential implications with a focus on nitrogen metabolism Emi Matsumura , Hiromi Kato , Shintaro Hara , Tsubasa Ohbayashi , Koji Ito , Ryo Shingubara , Tomoya Kawakami , Satoshi Mitsunobu , Tatsuya Saeki , Soichiro Tsuda , Kiwamu Minamisawa , Rota Wagai bioRxiv 2025.01.09.632267; doi: https://doi.org/10.1101/2025.01.09.632267 Share This Article: Copy Citation Tools Single-cell genomics of single soil aggregates: methodological assessment and potential implications with a focus on nitrogen metabolism Emi Matsumura , Hiromi Kato , Shintaro Hara , Tsubasa Ohbayashi , Koji Ito , Ryo Shingubara , Tomoya Kawakami , Satoshi Mitsunobu , Tatsuya Saeki , Soichiro Tsuda , Kiwamu Minamisawa , Rota Wagai bioRxiv 2025.01.09.632267; doi: https://doi.org/10.1101/2025.01.09.632267 Citation Manager Formats BibTeX Bookends EasyBib EndNote (tagged) EndNote 8 (xml) Medlars Mendeley Papers RefWorks Tagged Ref Manager RIS Zotero Tweet Widget Facebook Like Google Plus One Subject Area Microbiology Subject Areas All Articles Animal Behavior and Cognition (7629) Biochemistry (17660) Bioengineering (13881) Bioinformatics (41913) Biophysics (21436) Cancer Biology (18578) Cell Biology (25482) Clinical Trials (138) Developmental Biology (13372) Ecology (19889) Epidemiology (2067) Evolutionary Biology (24302) Genetics (15599) Genomics (22483) Immunology (17728) Microbiology (40365) Molecular Biology (17163) Neuroscience (88540) Paleontology (666) Pathology (2830) Pharmacology and Toxicology (4821) Physiology (7637) Plant Biology (15136) Scientific Communication and Education (2045) Synthetic Biology (4290) Systems Biology (9818) Zoology (2269)