Intertillage during Natural Farming Rice Paddy Production Negatively Impacted the Microbial Abundances in Soils but Not Diversities

In natural farming rice paddies, intertillage (tillage between rows, during rice growth period) is often performed mainly to remove weeds without the use of chemicals. Also, the intertillage disturbs soil surfaces, potentially impacting the characteristics of soil microbial communities, such as their diversity and abundance. Natural farming systems aim to maintain biodiversity, but it remains unclear whether the intertillage impacts soil microbes in rice paddies.,us, this study aimed to understand to what extent “five times intertillage” treatment (5T) influences soil bacterial abundance and community structures compared with no tillage (NT), under a natural farming rice paddy system. Soils were sampled at rice proximity, soil surface, and 10 cm depth in a natural farming rice paddy, during the early to late vegetative phase (June to July), in Hokkaido, Japan. ,e 16S rRNA community structures and abundance were analyzed by next generation sequencing (NGS) and quantitative PCR, respectively. We observed that NT had significantly higher bacterial abundances at the soil surface than 5T. However, there were no clear differences between 5T and NT, regarding the bacterial community structures, including their diversity indices. Instead, the sampling timings markedly impacted the bacterial community structures for the rice proximity and soil surface, showing increasing diversity indices at the late vegetative stage, compared to the early vegetative stage, suggesting the interaction between the crop growth and bacterial communities. In this study, we did not observe the significant difference between the rice yield from NT (2.3± 0.7 t·ha) and 5T plots (2.7± 0.9 t·ha); however, the 5Tmight have negatively impacted soil bacterial abundances but not the community structure of the bacteria.


Introduction
Natural farming is an alternative agricultural system to grow food, which is highly supported by microbial functions rather than utilization of various types of fertilizer [1][2][3][4]. e concept of natural farming is similar to conservation agriculture (often involves minimal tillage or no tillage (NT) and rotation) [5][6][7], yet differs from organic farming (still in need of plowing, tilting, spreading organic fertilizer, and weeding).
Regarding natural farming rice paddy systems, Kasubuchi et al. [8] concluded that conducting five times inter tillage (5T) in a natural farming rice paddy during the rice growing period can increase the yield as high as that of the conventional farming about five years later. Multiple intertillage methods had been widely conducted in rice paddies to control weed growth during Edo period which was from the end of 17 th -18 th century. Intertillage is performed during the rice growth period, using a small machine with chainweeder to control weeds amounts (Supplementary Photo 1), enhance soil aeration, plant root growth, and water infiltration [8].
Generally, tillage can soften soil for plant root to spread over a deeper distance underneath the soil surface; shatter soil organic matter to accelerate the release of nutrients for crop growth; remove weed to minimize its competition with crop; incorporate oxygen which benefits aerobic microorganisms active [9][10][11][12]. Contrastingly, tillage decreases microbial abundance and activities in a long-term period due to increasing decomposition rate of soil organic matter which contains dead animal body and plant [13,14]. In contrast to tillage treatment, NT treatment can increase formation of soil organic matter and soil aggregate because there is no disturbance on soil, thus further improving soil microbial activities [15][16][17][18][19].
However, the difference in microbial abundance and diversity under NT and intertillage treated soils in naturally managed rice paddy systems was not well understood. e intertillage is performed during the rice growth period, when the field is submerged in water; thus the impact of the tillage in the rice paddy can be different from the impact of tillage for the upland soils. erefore, better understanding of the impacts of intertillage on soil microbial abundance and diversity, specific for the natural farming rice paddy, is required to elucidate the reasons behind the successfulness of this system. e maintenance of soil microbial activities and functions is critically important for natural farming, in general.
us, the primary goals of this study were to compare the dissimilarity of soil bacterial abundance and diversity under 5T and NT treatment within a paddy farm managed using the natural farming style. We hypothesized the following: (1) soil bacterial abundance and community would sharply change after intertillage compared with NT; (2) 5T and NT treatment would be significantly at variance with soil bacterial abundance and community for the duration of vegetative phase.

Field Sites and Agricultural Management.
e study was conducted on flooded rice paddy research fields of the Field Science Centre of Hokkaido University, Sapporo, Japan (N43°04′39″N/E141°20′03″E) (Supplementary Figure S1). e trial investigates the dynamic of intertillage and NT treatments in a natural farming rice paddy. In this study, we focused on 5T and NT treated plots. e natural managed paddy is 40 × 25 m in size and divided in nine plots with 13.34 × 8.33 m in size each. e natural farming rice paddy alongside with no fertilizer, herbicide, pesticide, and insecticide application was established in 2017. Before 2017, there were three years maintained by flooded conditions during the cropping seasons with no crop growth and herbicide application. Also, the fields were ploughed at one month before seedlings were transplanted to the paddy fields. e treatment of the nine plots is divided in 0, 2, and 5 times intertillage (three plots each) (Supplementary Figure S2).

Soil Sampling Procedure.
In frame of the study, plots from the natural farming rice paddy with 5T and NT treatment were sampled. To investigate short-term effects of 5T treatment compared with NT, samples were taken one day before and two days after management activities for the third, fourth, and fifth intertillage event. e third and fifth intertillage events corresponded to early vegetative phase (June) and late vegetative phase (July), respectively. In total, six plots were sampled corresponding to three replicates per treatment. From each plot, soil samples were taken at three different locations within a square of two rows of five consecutive rice plants. e soil surface samples were taken from soil layer at 0-0.5 cm depth, between plant rows. e rice proximity samples were taken from the rhizosphere zone of the rice plants. A spatula was used to scrape off the soil attached to the rice roots. e rice roots were exposed at the soil surface (in the paddy water) thus the scraping could be performed easily. e 10 cm samples were taken from between the rice plant rows, similar to the soil surface samples, but at 9-11 cm depth. At 10 cm depth, the soils were not disturbed by the intertillage (Supplementary Figure S3).

16S Sequencing Library Preparation.
e DNA of the sampled soils was extracted using a NucleoSpin Soil Kit (Macherey-Nagel GmbH & Co. KG, Duren, Germany) according to the manufacturer's protocol, and buffer SL2 was used. A negative extraction control was performed using empty extraction tubes with the first and last extraction. e DNA extract was quantified using Qubit 2.0 Fluorometer (Invitrogen, Waltham, United States) and a Qubit dsDNA BR Assay Kit. e V4 region of the 16S rRNA gene was amplified using PCR primer F515 (5′-CACGGT CGKCGGCGCCATT-3′) and R806 (5′-GGACTAC HVGGGTWTCTAAT-3′). An AmpliTaq Gold ® 360 Master Mix (Applied Biosystems ™ , Carlsbad, USA) and 5 ng input DNA were used for the PCR. e PCR program included initial denaturation at 95°C for 10 min followed by 20 cycles of denaturation at 95°C for 30 sec, annealing at 57°C for 30 sec, and elongation at 72°C for 1 min and ended with final elongation at 72°C for 7 min. A positive control using E. coli DNA control (Applied Biosystems ™ , Carlsbad, USA) and a negative control were performed with the same conditions. e amplicons were checked on a 1.5% agarose gel. Purification of PCR products was performed using an Agencourt AMPure XP kit (Beckman Coulter Inc., Webster, United States) according to the manufacturer's protocol. Ten nanograms of amplified DNA per sample was barcoded using the IonA-barcode[i]-F515 forward and IonP1-R806 reverse primers. e PCR program included the same conditions like above but with 5 cycles. e barcoded amplicons were purified and quantity-checked as described above. e quality was analyzed on an Agilent 2100 Bioanalyzer (Agilent Technologies, Inc., Santa Clara, USA) using an Agilent DNA 1000 Kit. e library was diluted to 50 pM and loaded into the Ion 314 Chip ( ermo Fisher Scientific K.K., Japan) using the Ion Chef Instruments ( ermo Fisher Scientific K.K., Japan) with the Ion PGM ™ Hi-Q Chef Kit. DNA sequencing was conducted on the Ion PGM Sequencer ( ermo Fisher Scientific K.K., Japan) using the Ion PGM 400 Kit.

Gene Abundance Assay.
Abundance of the 16S rRNA gene was assessed by quantitative polymerase chain reaction (qPCR) using a Stratagene Mx3005P cycler (Agilent Technologies, Inc., Santa Clara, USA). Prior to the experiment, the qPCR assay was optimized using different oligonucleotide concentrations, soil dilutions, and annealing temperatures to reach R 2 > 0.999 and amplification efficiencies (Eff) between 0.8 and 1. e amplification efficiencies were calculated using the following formula: Eff � 10 (−1/slope) − 1.
e qPCR reactions were performed in 20 µl using KAPA SYBR green Master Mix (Takara Bio) and 400 nM of each primer F515 and R806. e cycler program was set with an initial enzyme activation step for 5 min at 95°C and annealing step for 1 min at 58°C. e DNA extracts were 100-fold diluted. e biological replicates were analyzed in technical duplicates. Negative controls and serial dilutions of amplified E. coli DNA (10 2 to 10 6 ng/µl) were included in every qPCR run to calculate standard curves for absolute quantification. Melting curve analyses were performed to check the quality of the generated amplicons.

Measurement of Ammonium and pH.
e determination of the inorganic-N concentrations was performed following the approaches taken by Silva and Hatton and Pickering [20,21]. Firstly, 2 g of the sampled fresh soil was extracted with 2M KCl (10 ml). After shaking the mixture for 60 min, the suspension was filtered through a filter paper (Grade 5C, <5 mm; Advantec, Tokyo, Japan). e extracted solution was stored at minus 20°C until measurement. For the measurement, a colorimetric method was employed with a flow injection analyzer (AQLA-700; Aqualab, Tokyo, Japan) [22,23]. To measure soil pH, 5 g of fresh soil was shaken with 25 ml of 10% KCl for 30 minutes. e pH of the extracts was measured by pH meter (AS800, AS ONE Corporation, Osaka, Japan).

Analysis of the 16S rRNA-Based Bacterial Community
Structures.
e barcoded 16S rRNA gene sequences were demultiplexed, quality-filtered, and assessed using the Quantitative Insights Into Microbial Ecology (QIIME) workflow [24]. e next generation sequencing based method showed the loading density was at 87% and 94%. On average, 58,889 reads were mapped per 16S rRNA sample. Operational Taxonomic Units (OTUs) were prepared by eliminating all the OTUs that matched the GreenGenes 13_5 reference sequence with 97% similarity. To analyze the changes in microbial community and their interactions with other environmental factors, Principal Component Analysis (PCA) was performed; the PCA plots were separated into three parts with same variation. Soil microbial community analysis was performed with permutational multivariate analysis of variance (PERMANOVA) and the Bray-Curtis distance based on 999 permutations of the raw data using the Adonis function in R.

Statistical Analysis.
For the pH and concentration of ammonium, two-way and three-way analysis of variance (ANOVA) were performed to investigate the effect of the location and the duration. One-way ANOVA and Tukey-Kramer multiple comparison tests were performed for the analysis of significant differences among samples on duration or for each treatment. Statistical analyses were performed by R 3.6.1 [25] and we set the significance level at P < 0.05.

Bacterial Abundance.
At the soil surface, the bacterial abundance in the NTsoils was higher than that in the 5Tsoils (P < 0.05); however, the abundances of bacteria in soils were similar between NT and 5T at the rice proximity and 10 cm depth (Figure 1). e bacterial abundance was higher at the rice proximity and soil surface, compared with 10 cm (P < 0.05).
e bacterial abundance was decreasing over time for the soil surface and rice proximity, but it was relatively stable over time at the 10 cm depth soils.
When averaged across the sampling timings, NT soils had a larger amount of DNA, when compared with 5T soils, in the soil surface, but this was not observed for the rice proximity and the 10 cm depth soils ( Figure 2). e raw soil DNA concentrations increased over time at the rice proximity and soil surface under NT treatment. For the 5T at the rice proximity and soil surface, the raw DNA concentrations peaked at day 32 and did not show increase towards the later stage of the experiment. At the 10 cm depth, the raw DNA concentrations under NT and 5T treatment had the similar fluctuation and they were lower than the other locations.

Bacterial Community.
On early vegetative phase, the bacterial community structure was similar at the same location regardless of the tillage treatment ( Figure 3). However, there was a significant difference between 10 cm depth and rice proximity as well as soil surface (Supplementary  Table S1). At the rice proximity location, NT treatment had a higher bacterial relative abundance in Chloroflexi (9.91%), Gemmatimonadetes (6.07%), and Spirochaetes (10.66%) compared to the 5T, while Acidobacteria (8.33%) and Betaproteobacteria (8.63%) were relatively higher in the 5T compared to the NT treatment. At soil surface location, NT treatment had a higher Spirochaetes (10.11%) bacterial abundance, but Chloroflexi (12.26%) had a higher relative abundance under 5T treatment, compared to the NT. At 10 cm depth location, Firmicutes (17.71%) was higher under NT treatment, but Planctomycetes (3.37%) and Spirochaetes (9.55%) were higher under 5T treatment, compared to the NT. e N 2 -fixing genus only appeared at the rice proximity and soil surface, while the relative abundances of Firmicutes and Archaea predominated 10 cm depth location. Nonetheless, there were no specific bacteria predominating rice proximity and soil surface locations.
On late vegetative phase, the bacterial community structure (Figure 3) was alike among three locations. ere were no bacteria particularly predominating rice proximity, soil surface, and 10 cm depth locations. Nevertheless, we observed that several bacterial clusters partly fluctuated from  the early vegetative phase to the late vegetative phase. e relative abundance of Deltaproteobacteria Geobacter increased obviously at the rice proximity and soil surface. Contrastingly, Betaproteobacteria, Firmicutes, and Pseudomonas had decreased at the two locations. e 10 cm depth was a relatively steady location compared to the rice proximity and soil surface. e result of PERMANOVA showed soil bacterial community structure was significantly influenced by vegetative phase (p < 0.001) rather than treatment (Supplementary Table S1). e Shannon diversity index was also significantly influenced by the vegetative phase (P < 0.001) (Supplementary Table S2). Its result showed that the bacterial diversity in the rice proximity and the soil surface increased over time, but it was maintained at the 10 cm (Supplementary Table S3).
Based on the PCA, there were no clear clusters of NT and 5T (Figure 4(a)), suggesting there were no major impacts of the intertillage on the bacterial community structures. For the locations of rice proximity and soil surface, they clustered differently by the two vegetative phases (Figure 4(b)). In contrast to the rice proximity and soil surface locations, the 10 cm depth data made a separate cluster on the PCA, and the NT and 5T showed different clusters for the 10 cm depth. e PCA on genus level (Figure 4(c)) indicated that genus Clostridium predominated the rice proximity and surface under both treatments on early vegetative phase. Stramenopiles (or Heterokonts), a microeukaryotic community whose most were algae [26], predominated rice proximity and surface under both treatments on late vegetative phase. With regard to 10 cm depth location, the uncultured order pGrfC26 of the Miscellaneous Crenarchaeota Group (MCG) was predominant.

Dynamics of pH and Ammonium.
When averaged across the rice vegetative phase, soil pH was the highest on the soil surface (6.8 ± 0.01), followed by the rice proximity  Applied and Environmental Soil Science (6.5 ± 0.02) and the 10 cm depth (6.2 ± 0.02). Soil pH decreased from 21 days to 43 days under NT and 5T treatment at both rice proximity and soil surface location. e pH of rice proximity and soil surface under NT and 5T treatment decreased about 5%. Soil pH at 10 cm depth location under NT and 5T treatment was relatively stable, ranging from 6.25 to 6.08. ere was a significant difference between location and duration (P < 0.001) ( Figure 5). However, soil pH whether disturbed (5T) or not (NT) did not show obvious significant discrepancy. e ammonium declined from 21st day at the early growth after transplanting at the rice proximity and soil surface ( Figure 6). ere was a significant difference between days (P < 0.001), but there was no significant impact on the tillage treatments. e amount of ammonium concentration declined at rice proximity (23 to 9 mg·kg −1 , NT; 19 to 7 mg·kg −1 , 5T) and soil surface (25 to 11 mg·kg −1 , NT; 21 to 9 mg·kg −1 , 5T) from 24 days to 32 days and afterwards slightly increased by 43 days after transplanting. ere was a significant difference at the 10 cm depth location between NT and 5T treatment. e ammonium concentration at the 10 cm depth was relatively stable compared with the rice proximity and soil surface.

Discussion
Our study found that there was a negative impact of the tillage on the bacterial abundances (P < 0.05) at the soil surface ( Figure 1). However, the tillage treatment did not impact the bacterial abundance in soils at the rice proximity and 10 cm depth. We note that the rice grain yields for the 0 and 5 times intertillage treatments were 2.3 ± 0.7 and 2.7 ± 0.9 t·ha −1 , respectively, without the significant difference between them. e tillage disturbs soil surfaces, and the disturbance must have negatively influenced the soil microbial abundance on the soil surfaces in our study. Previous studies also had shown that a tillage treatment results in lower microbial biomass than NT treatment at the soil surface [27,28]. Previous studies stated that the negative impact of the tillage on soil microbial biomass was due to the  interaction between duration (early vegetative phase (gray) and late vegetative phase (black)) and location (rice proximity, soil surface, and 10 cm depth). (c) Arrows show bacteria on genus level. Axes 1 and 2 account for 72 and 16% of the variation, respectively. 6 Applied and Environmental Soil Science  disruption of microbial extracellular enzyme activities. For example, one of the enzymes in soils, Βeta-Glucosidase activity, is positively related to total soil carbon that is closely associated with soil microbial biomass [29][30][31]. In contrary to tillage, NT treatment retains the soil carbon in soils, when compared to the systems with tillage, and this allows microorganisms to grow [32]. Schmidt et al. [33] also indicated that, at a 0-5 cm soil depth, the NT treatment led to higher copy number of bacteria that was approximately 7.5 × 10 8 copy kg −1 soil (this study was from 4.9 × 10 7 to 1.4 × 10 8 copy kg −1 soil at the soil surface) when compared to that under the tillage treatment. However, the same study indicated that this was not the case at deeper soil layers. us, we conclude that the disturbance of soil microbes occurred in rice paddy soils by intertillage occurrence, even when the soils were submerged in water, similar to previous reports performed in upland soils.
Also, we found that the 16S rRNA copy number at the location of rice proximity and soil surface under NT and 5T treatment decreased along with rice growth (Figure 1), but at the 10 cm depth it was more stable compared to other locations. is was consistent with the previous studies that 16S rRNA copy number decreased with increasing plant growth [31,34]. We also observed the highest amount of 16S rRNA copy number appeared on the early vegetative phase.
us, we assumed that one of the reasons for the decrease of bacterial abundances might be caused by rice plants rapidly consuming nutrients resulting in nutrient deficiency for microbial growth [33,35,36]. erefore, the effect of intertillage on the microbial growth may depend on the nutrient availability in the soils; thus the interaction between the soil's nutrient availability and intertillage treatments should be further studied, regarding the soil microbial abundance. In the current study, we aimed to compare tillage and no tillage systems within a natural farming rice paddy. However, future studies should compare the natural farming and conventional systems, regarding the relationships between plant growth and microbial abundance, because we were uncertain whether the use of chemical fertilizers and other practices in the conventional farming systems change the relationship between plant growth and microbial abundance.
As another factor controlling the microbial abundance in the current study, the bacterial abundance peaked at around pH of 7 and decreased with soil pH at the rice proximity and soil surface both under NT and 5T treatments (Figures 1 and  5). is agreed to a previous study reporting that the bacterial growth was correlated with soil pH and tended to be relatively higher at around neutral soil pH environment [37]. However, unexpectedly, our results showed an increase of the raw DNA concentrations with decreasing soil pH and over the experimental period, particularly for the NT (Figure 2). Other biological activities such as fungal, faunal, or plant activities might have contributed to the increase of raw soil DNA amount; for example, previous studies reported that acidic environments promoted fungal activities, while bacterial activities decreased along with reducing soil pH [35,38]. Also, NT systems can improve the soil porosity by increasing 0.5-50 mm macropores, as well as altering water holding capacity (0-10 cm) and the amount of exchangeable ions, compared with tillage systems [35,38,39].
us, further investigation is needed to identify the factors contributing to the increase of the raw DNA concentrations.
For the community structures of bacteria in the soils, there was no clear impact of intertillage treatments, when compared to NT (Figure 4(a)). Rather, the sampling timing influenced the community structures, particularly at soil surface and rice proximity (Figure 4(b)). For example, Clostridium of Firmicutes in rice proximity and surface was present in relatively larger amount during the early vegetative phase compared to the late vegetative phase (Figure 4(c)). Clostridium is a microorganism responsible of the rice straw decomposition in a rice paddy [40,41]. Due to low temperatures in Hokkaido limiting the microbial decomposition ability, we often find a large amount of rice straw left beneath the soil surface in spring.
us, ours results suggested that the activated Clostridium might decompose some of the rice straw in spring. e contribution of Clostridium to the decomposition of rice straw has to be further studied, especially in relation to their low temperature activities.
Contrastingly, Stramenopiles dominated the rice proximity and soil surface soils during the late vegetative phase (Figure 4). Stramenopiles are an assemblage of eukaryotic organism, including unicellular, such as diatoms, to large multicellular forms, such as the brown algae and oomycetes.
ese Microeukaryotes act as grazers which can affect microbial biomass and community as well as releasing ammonium and nitrate [42]. us, these microbes might play an important role in natural farming rice paddy systems, providing available nitrogen to rice plants in the later stages of their growth, although further studies are needed to confirm this.
In the 10 cm depth soils, class MCG of Crenarchaeota showed relatively larger amount, when compared to the soils sampled from other parts (Figure 4(c)). Miscellaneous Crenarchaeota Group (MCG) of Archaea are widely distributed in terrestrial and marine ecosystem [43]. ese microorganisms play an important role in biogeochemical cycles. MCG in anaerobic soil ecosystem are related to degrade carbon from plant, reduce nitrite to ammonium, and produce acetate [44]. Undecomposed rice straw, which was found while soil sampling at the 10 cm depth under both treatments, possibly led MCG to be in the majority as well as lower pH and higher ammonium concentration compared to the rice proximity and soil surface.
Overall diversities of the bacterial communities also did not show the difference between the tillage treatments (Supplementary Table S2). However, the rice proximity and soil surface under NT and 5T treatment showed a significant increase of microbial diversity on late vegetative phase, compared to the 10 cm soils (Supplementary Table S3).
ere were previous studies [45][46][47] that found that NT treatments increased microbial community's diversity when compared to the soils under tillage, but other previous studies concluded that the NT treatment did not influence the microbial community [48,49]. e microbial diversities' fluctuations over the growing period might have been the reasons behind these contrasting results because microbes are influenced by fluctuating soil chemical properties such as pH and nutrients [50].
We expected different soil pH and nutrient characteristics at soil surfaces and near rice plants between the NT and 5T treatments because previous studies stated that, in natural farming rice paddy systems, intertillage can incorporate oxygen to aerate soil, therefore promoting microorganism activities to decompose organic matter and crop residues, facilitating nitrogen mineralization and ammonification [8,9,11,51]. For the ammonium concentrations, we observed relatively higher ammonium concentration in 5T than NT only in a deeper soil layer (P < 0.01). Possibly, the rapid consumption of ammonium ion by rice plants masked the differences between the NT and 5T, regarding the soils' chemical characteristics. In results, we observed the decrease in soil pH and ammonium in both of the tillage treatments and we conclude that the increase in the diversity indices with time might have been caused by increased relative abundances of acid preference bacteria, such as Acidobacteria and Firmicutes [52]. We note that our natural farming rice paddy experiment was only conducted for one year; thus the significant impact of the five times intertillage treatments (compared to NT) on soils' microbial diversity might appear more in a longer term, although it needs to be confirmed.

Conclusion
From this study, we draw the following conclusions: (i) Bacterial abundance was higher at soil surface with no tillage, compared to the soils under intertillage. Soil surface was the most disturbed location due to the tillage; thus we observed the significant decrease in bacterial abundance in this zone, when the tillage was conducted. (ii) Bacterial abundance was reduced over time, regardless of the tillage treatments, at rice proximity and soil surface. Decreasing nutrients and pH because of the plants' nutrient absorbance limited the bacterial growth and therefore led to bacterial abundance reduction. (iii) e bacterial diversity and community structure were affected by soil pH which was decreasing with ammonium levels in soil. However, the impacts of the tillage treatments were not clear.

Data Availability
e raw and analyzed data used to support this study are available from the corresponding author upon request.

Conflicts of Interest
e authors declare that there are no conflicts of interest or personal relationships that could have influenced the work reported in this paper.

Acknowledgments
is study has received grants from the JSPS KAKENHI 443 Grants-in-Aid for Scientific Research (Grant number: JP 18H03963). e authors thank the laboratory members for helping them with sampling and lab works.

Supplementary Materials
Supplementary photo 1: the intertillage machine used to remove weed in five times intertillage plots. Figure S1: the sampling site: e Field Science Centre of Hokkaido University, Sapporo, Japan (N43°04′39″151, E141°20′03″634). Figure S2: the experimental site of natural farming rice paddy. It was divided into nine plots and A, B, and C blocks. Panels were set within block (brown) and between blocks (blue). Each plot was separated with polycarbonate board. Each block includes 0 (no tillage; NT), 2 (2 times intertillage; 2T), and 5 (5 times intertillage; 5T) treatments, where soils were sampled at NT and 5T plots. Each treatment had three replicates. Figure S3: soil samples were taken from ① rice proximity; ② soil surface; ③ 10 cm depth one day before and two days after management activities for the third, fourth, and fifth intertillage event in 5T and NT plots. Table  S1: PERMANOVA (Permutational multivariate analysis of variance) tested effect results of vegetative phase, location, and treatment on bacterial community composition. Table  S2: results of ANOVA for Shannon index at the OTU level.