Responses of Ammonia-Oxidizing Microorganisms to Intercropping Systems in Different Seasons

: Intercropping plays an essential role in agricultural production, impacting the soil’s physical and chemical properties and microbial communities. However, the responses of ammonia-oxidizing microorganisms in the continuous-cropping soil to different intercropping systems in different growing seasons are still insufﬁciently studied. Here, we investigated the effects of seven intercropping systems (alfalfa ( Medicago sativa L.)/cucumber, trifolium ( Trifolium repens L.)/cucumber, wheat ( Triticum aestivum L.)/cucumber, rye ( Secale cereale L.)/cucumber, chrysanthemum ( Chrysanthemum coronrium L.)/cucumber, rape ( Brassica campestris L.)/cucumber, mustard ( Brassica juncea L.)/cucumber) on soil physical and chemical properties, potential nitriﬁcation rate (PNR), soil ammonia-oxidizing archaea (AOA), and ammonia-oxidizing bacteria (AOB) communities in the greenhouse in spring and autumn. The results showed that, compared with cucumber monoculture, intercropping increased the soil NH + 4 -N and NO − 3 -N. The chrysanthemum–cucumber, rape–cucumber, and mustard–cucumber treatments increased soil PNR. Intercropping increased the AOA and AOB abundances in two seasons, especially in rape–cucumber, wheat–cucumber, chrysanthemum–cucumber, and trifolium–cucumber treatments. The ratio of AOA and AOB decreased with seasonal variation. The wheat–cucumber and rape–cucumber treatments increased soil AOA community diversity. Seasonal variation had a signiﬁcant effect on the relative abundance of the AOB community. Nonmetric multidimensional scaling analysis showed that the AOA and AOB community structures were obviously different from spring to autumn. Redundancy analysis showed that the AOA community was signiﬁcantly regulated by moisture, NO − 3 -N, and available potassium (AK), while the AOB community was signiﬁcantly regulated by moisture, available phosphorus (AP), AK, NO − 3 -N, and pH. Network analysis showed that the co-occurrence relationship and complexity of AOA and AOB communities were different in two growing seasons. The AOB community may play a critical role in ammonia oxidation in autumn. Taken together, intercropping improved soil physicochemical state, increased soil PNR and signiﬁcantly altered soil AOA and AOB communities. Seasonal variation signiﬁcantly altered the AOA and AOB communities’ structure and interaction between them. The effect of seasonal variation on AOA and AOB communities was greater than intercropping.


Introduction
Intercropping is a commonly used agricultural practice, which has been proven to have significant advantages in yield increase, disease control, and resource utilization [1,2]. Simultaneously, intercropping changes the plant-microbe or microbe-microbe interaction, thereby affecting the nitrogen transformation in the soil [3]. Nitrification is a crucial step in the nitrogen cycle and plant nutrition. It is a biology-mediated process and plays a crucial role in nitrogen form transformation and the nitrogen cycle. It directly affects environmental quality and nitrogen utilization efficiency [4,5]. Ammonia oxidation and nitrite oxidation are two critical steps of nitrification. Ammonia oxidation is the rate-limiting step of nitrification, which is mainly completed with the participation of soil ammonia-oxidizing archaea (AOA) and ammonia-oxidizing bacteria (AOB) [6,7]. Ammonia-oxidizing microorganisms can respond to changes in the soil micro-habitat of the intercropping system. For example, maize-faba bean intercropping increased the abundances of total archaea and the AOA and AOB community diversities and structures were affected by plant growth period and competition among plant species, which changed the soil nitrification rate [8]. Pepper-alfalfa intercropping and pepper-maize intercropping stimulated the activity of AOB in the topsoil, while AOA decreased with planting along with the soil profile and this deep and shallow root intercropping system could absorb nitrogen more effectively [9]. Previous studies have also confirmed that seasonal variation can also change the abundance and structure of ammonia-oxidizing microorganisms [10,11]. However, the effects of intercropping on AOA and AOB communities in the continuous-cropping soil are poorly understood in different seasons.
Due to the increase in plant diversity, intercropping often affects soil properties. Many studies have shown that ammonia-oxidizing microorganisms were affected by many soil environmental factors [12,13]. For example, the AOA and AOB community structures were regulated by the soil moisture, pH, and TN in wheat-maize rotation soil [14]. In low-pH and ammonia nitrogen soil environments, AOA were most represented and functionally dominant, while in neutral, alkaline or nitrogen-rich soil, AOB were dominant in nitrification [15,16]. The soil pH directly affects the effectiveness of NO − 3 -N as a substrate for ammonia oxidation. The soil nitrification rate represents the ability of soil nitrification microorganisms to convert NH +4 into NO − 3 , so it is also an essential indicator of the abundances of nitrification microorganisms [17,18]. The availability of NH +4 is also an important factor determining the ecological niche of AOA and AOB [6,19]. AOA are more effective in utilizing low concentrations of NH +4 produced by organic nitrogen mineralization [20]. In the range of high NH +4 concentrations, the AOB abundance and contribution to ammonia oxidation are positively correlated with the NH +4 concentration [21,22]. AOA are more widely distributed than AOB in soil, which means that AOA may contribute more to ammonia oxidation [23,24]. However, some studies have also confirmed that AOA were dominant in abundance, while AOB were dominant in function in intercropping systems [14]. Therefore, the effects of different intercropping systems on the continuous-cropping soil ammonia-oxidizing microorganisms need to be studied in depth.
Cucumber (Cucumis sativus L.) is a nitrogen demanding crop. The availability of nitrogen in the soil is critical for the growth of cucumber. Because continuous cropping obstacle is common in facilities, the availability of soil nutrients will be reduced, which is not conducive to the growth of crops. However, intercropping can alleviate the obstacle of continuous cropping, increase the abundances of microorganisms in the soil, improve the community structure, accelerate the transformation of nutrients in the soil, and improve the quality of crop nutrient absorption [25][26][27]. However, the interaction of ammonia-oxidizing microorganisms in the continuous-cropping soil of intercropping system in the greenhouse is still unclear.
In this study, the physical and chemical properties, PNR, and ammonia-oxidizing microbial communities of the continuous-cropping soil were studied in cucumber monoculture and cucumber intercropping with seven crops from the Gramineae, Leguminosae, Cruciferae, and Compositae families in spring and autumn. The abundances, diversities, and composition of soil AOA and AOB communities in monoculture and intercropping systems of cucumber were analyzed by quantitative PCR (qPCR) and Illumina MiSeq sequencing to assess the potential effects of different intercropping systems and planting seasons on soil AOA and AOB communities. Our aim was to determine the effects of seven intercropping systems on soil AOA and AOB in two cropping seasons. We hypothesized Agriculture 2021, 11,195 3 of 17 that (1) the structure and composition of soil AOA and AOB communities could be changed by intercropping seven crops with cucumber, (2) the difference of physical and chemical factors caused by seasonal variation changed the community structure and coexistence mode of ammonia-oxidizing microorganisms.

Experimental Location and Design
The experiment was carried out in a greenhouse at the experimental base of Horticulture of Northeast Agricultural University (45 •  A single-factor random block design was adopted, and each treatment was repeated with three replicates in the greenhouse in spring (21 April-30 June) and autumn (29 July-6 October) in 2015. The size of each plot was 6 × 0.5 m 2 . Seven crops (alfalfa (A), trifolium (T), wheat (W), rye (Ry), chrysanthemum (C), rape (Ra), and mustard (M)) were intercropped with cucumbers as the experimental group, and a cucumber monoculture (CM) was the control. Each plot consisted of one row of 12 cucumber plants with two protection lines on each side. When cucumbers were planted for 10 days, intercropped seeds were sown approximately 10 cm away from cucumber seedlings ( Figure 1). The number of intercrop seedlings for wheat and rye was 30 seedlings, that for trifolium and alfalfa was 40 seedlings, that for rape and mustard was 3 seedlings, and that for chrysanthemum was 5 seedlings. When the intercropping crops reached 20 cm, to avoid an impact on cucumber growth, the intercropping crops were trimmed by 10 cm. Compound NPK fertilizer (16:16:8) was applied to each plot at a rate of 300 kg ha −1 .
systems of cucumber were analyzed by quantitative PCR (qPCR) and Illumina MiSe quencing to assess the potential effects of different intercropping systems and plan seasons on soil AOA and AOB communities. Our aim was to determine the effec seven intercropping systems on soil AOA and AOB in two cropping seasons. We hyp esized that (1) the structure and composition of soil AOA and AOB communities coul changed by intercropping seven crops with cucumber, (2) the difference of physical chemical factors caused by seasonal variation changed the community structure and existence mode of ammonia-oxidizing microorganisms.

Experimental Location and Design
The experiment was carried out in a greenhouse at the experimental base of Hort ture of Northeast Agricultural University (45°41′ N, 126°37′ E), Harbin, Heilongjiang P ince, China. The soil was a black soil (Mollisol) under continuous cucumber croppin 3 years with a sandy loam texture, and contained sand, 326 g kg −1 ; silt, 384 g kg −1 ; clay g kg −1 ; organic matter, 3.51%; pH, 6.64; NH A single-factor random block design was adopted, and each treatment was repe with three replicates in the greenhouse in spring (April 21-June 30) and autumn (July October 6) in 2015. The size of each plot was 6 × 0.5 m 2 . Seven crops (alfalfa (A), trifo (T), wheat (W), rye (Ry), chrysanthemum (C), rape (Ra), and mustard (M)) were i cropped with cucumbers as the experimental group, and a cucumber monoculture ( was the control. Each plot consisted of one row of 12 cucumber plants with two protec lines on each side. When cucumbers were planted for 10 days, intercropped seeds w sown approximately 10 cm away from cucumber seedlings ( Figure 1). The number o tercrop seedlings for wheat and rye was 30 seedlings, that for trifolium and alfalfa w seedlings, that for rape and mustard was 3 seedlings, and that for chrysanthemum w seedlings. When the intercropping crops reached 20 cm, to avoid an impact on cucum growth, the intercropping crops were trimmed by 10 cm. Compound NPK ferti (16:16:8) was applied to each plot at a rate of 300 kg ha −1 .

Soil Sampling and Analysis
The soil samples were obtained from 24 test districts in the greenhouse on 30 and 6 October, the collection depth was 0-20 cm, and then they were mixed and si

Soil Sampling and Analysis
The soil samples were obtained from 24 test districts in the greenhouse on 30 June and 6 October, the collection depth was 0-20 cm, and then they were mixed and sieved through a 2 mm mesh. Afterward, the soil sample was divided into two parts in the laboratory: one was kept at −80 • C for DNA extraction, and the other at −4 • C for chemical analysis. For soil inorganic nitrogen concentration (NH +4 -N, NO − 3 -N) analysis; each sample was extracted by 2 M KCl, and the soil available phosphorus (AP) and available potassium (AK) were extracted by 0.5 M NaHCO 3 (pH = 8.5) and 1 M NH 4 + OAc (pH = 7.0), respectively. The soil filtrates were measured with a continuous flow analyzer (SAN++, Skalar, Breda, Netherlands). Soil pH was determined by inserting electrodes into a soil-water suspension (1:2.5). The soil moisture content was calculated by drying at 105 • C for 8 h [28].

Potential Nitrification Rate (PNR)
PNR was analyzed as previously described [29]. Briefly, 5 g soil was put into per conical flask respectively, and 20 mL of phosphate buffer solution (1 mM, pH 7.4) and 0.1 mL of NaClO 3 solution (1.5 M) were added. They were gently mixed and sealed with sealing film. The two conical flasks were put onto a shaking table for 5 h (samples). The other sample was stored at −20 • C for 5 h (control). After culturing, the control was dissolved at room temperature, and then 5 mL of 2M KCl solution was added to the sample and control and filtered immediately. Next, 5 mL of filtrate, 3 mL of NH 4 Cl buffer solution, and 2 mL of color-developing solution were added to the test tube and mixed evenly. After being placed at room temperature for 15 min, the color was determined spectrophotometrically at 520 nm.

DNA Extraction and Quantitative PCR (qPCR)
According to the manufacturer's instructions, the total DNA of the soil was extracted from 0.5 g soil using the MoBioPowerSoil TM DNA Isolation Kit (Mobio Laboratories Inc. Carlsbad, CA, USA). The DNA concentration and purity were determined with a NanoDrop 2000 Spectrophotometer (Thermo Scientific, Waltham, MA, USA).

Illumina MiSeq Sequencing and Data Processing
The PCR amplification primers were the same as qPCR, and the Barcode sequence was added at the front of the primer to distinguish samples during amplification. The amplification was carried out in 25 uL volume, including 2 µL DNA extract, 1 µL each of the front and back primers, 5 µL 5 × PCR buffer, 2 µL (2.5 mmol·L −1 ) dNTP, 5 µL 5 × Q5High GC Enhancer buffer, 0.25 µL (0.02U·µL −1 ) Q5High-Fidelity DNA polymerase (NEB) and 8.75 µL ddH 2 O. The reaction conditions were as follows: 98 • C for 5 min, 35 amplification cycles of 98 • C for 30 s, 55 • C for 30 s, 72 • C for 45 s. The amplified PCR products were detected by 2% agarose gel electrophoresis, and AxyPrep DNA Gel Extraction Kit (Axygen Biosciences, Union City, CA, USA) were used to recover the products. The PCR products were analyzed by QuantiFluor™-ST (Promega, Madison, WI, USA) blue fluorescence quantitative system, and the samples were mixed in proportion. Then, the high-throughput amplicon sequencing was carried out on Illumina MiSeq platform (2 × 300) at Majorbio Bio-Pharm Technology Co., Ltd. (Shanghai, China), each treatment was done in triplicate.
The raw sequence files were quality filtered and processed using QIIME (Version 1.8.0). The operational taxonomic units (OTUs) were clustered at 97% sequence similarity using UPARSE (Version 7.1) [32], and the chimeric sequences were detected and removed using UCHIME [33]. The archaeal and bacterial OTUs were classified by using ribosomal database project (RDP) classifier against the functional genomics resource (FGR) functional gene database [34]. To preclude bias due to several sequencing depths, all samples were subsequently subsampled based on the minimum number of soil microbial sequencing Agriculture 2021, 11,195 5 of 17 depths of this study. All sequences have been deposited in the NCBI SRA database (AOA: SRP 160125; AOB: SRP160403).

Statistical Analysis
The alpha diversity indices of the samples were generated by QIIME [35]. For beta diversity, the AOA and AOB community structures were analyzed by nonmetric multidimensional scaling (NMDS) based on the Bray-Curtis distance dissimilarity, and ANOSIM and Adonis were used to compare the microbial community differences of two growing seasons with the Bray-Curtis distance and 999 permutations, and the soil physical and chemical factors less than 10 were screened by variance inflation factor (VIF) for redundancy analysis (RDA), and they were performed using the "vegan" package in "R" (Version 3.3.1, R Foundation for Statistical Computing, Vienna, Austria). All data were checked for normality and homogeneity of variance. Network analysis based on spearman's rank analysis was performed using the 50 most abundant OTUs of AOA and AOB communities. Based on strong (ρ > 0.8) and significant correlations (p < 0.01), the symbiotic patterns of ammonia oxidizing microorganisms in different seasons were visualized by Gephi (Version 0.9.2). The correlation and analysis of variance were performed at 0.05 probability level in SPSS (Version 17.0).

Soil Physical and Chemical Properties and PNR
The physical and chemical properties of the different treatments were illustrated in Table 1. Compared with cucumber monoculture, all intercropping treatments significantly reduced the soil AP in spring (p < 0.05) but not in autumn. The rye-cucumber treatment significantly increased the soil NH +4 -N in spring (p < 0.05), and the alfalfa-cucumber and chrysanthemum-cucumber treatments significantly increased the soil NO − 3 -N in autumn (p < 0.05). There was no significant difference in the soil moisture, pH, or AK between intercropping and monoculture. The PNR of different treatments was shown in Figure 2. Except for the PNR of the alfalfa-cucumber treatment, other intercropping treatments were significantly lower than the cucumber monoculture in spring (p < 0.05). The PNR of the mustard-cucumber, rapecucumber, and chrysanthemum-cucumber treatments were significantly higher than the cucumber monoculture in autumn (p < 0.05). T-trifolium; W-wheat; Ry-rye; C-chrysanthemum; Ra-rape; M-mustard; CM-cucumber monocul-, AK indicate soil available phosphorus and available potassium, respectively. Columns with different letters are ntly different at 0.05 probability level.

Soil AOA and AOB Abundances
The AOA abundance of the wheat-cucumber, chrysanthemum-cucumber, rape-cucumber, and mustard-cucumber treatments were significantly higher than the cucumber monoculture (p < 0.05), but that of the rye-cucumber treatment was significantly lower compared to the cucumber monoculture in spring (p < 0.05). The AOA abundance of the alfalfa-cucumber, trifolium-cucumber, and wheat-cucumber treatments were significantly higher than the cucumber monoculture in autumn (p < 0.05). The AOA abundance of the wheat-cucumber treatment was the highest in both growing seasons (Figure 3a).
The AOB abundance of the chrysanthemum-cucumber, rye-cucumber, and alfalfacucumber treatments were significantly higher than the cucumber monoculture in spring (p < 0.05). The AOB abundance of the trifolium-cucumber treatment was significantly higher than the cucumber monoculture (p < 0.05), while that of the rye-cucumber, chrysanthemum-cucumber, rape-cucumber, and mustard-cucumber treatments were significantly lower than the cucumber monoculture in autumn (p < 0.05) (Figure 3b).

Soil AOA and AOB Abundances
The AOA abundance of the wheat-cucumber, chrysanthemum-cucumber, rapecucumber, and mustard-cucumber treatments were significantly higher than the cucumber monoculture (p < 0.05), but that of the rye-cucumber treatment was significantly lower compared to the cucumber monoculture in spring (p < 0.05). The AOA abundance of the alfalfa-cucumber, trifolium-cucumber, and wheat-cucumber treatments were significantly higher than the cucumber monoculture in autumn (p < 0.05). The AOA abundance of the wheat-cucumber treatment was the highest in both growing seasons (Figure 3a).
The AOB abundance of the chrysanthemum-cucumber, rye-cucumber, and alfalfacucumber treatments were significantly higher than the cucumber monoculture in spring (p < 0.05). The AOB abundance of the trifolium-cucumber treatment was significantly higher than the cucumber monoculture (p < 0.05), while that of the rye-cucumber, chrysanthemumcucumber, rape-cucumber, and mustard-cucumber treatments were significantly lower than the cucumber monoculture in autumn (p < 0.05) (Figure 3b).
In addition, the soil AOA abundance of all treatments was significantly higher than the AOB in both growing seasons (p < 0.05), and the ratio AOA/AOB was 9.16 to 58.35 in spring and 9.63 to 31.50 in autumn (Figure 3c,d).
The correlation analysis showed that the AOA abundance was negatively correlated to the moisture and positively correlated to the NO − 3 -N (p < 0.05), and the AOB abundance was negatively correlated to the pH and positively correlated to the NO − 3 -N in spring (p < 0.01). The AOA abundance was negatively correlated to the NO − 3 -N, AP, and AK (p < 0.05), and the AOB abundance was positively correlated to the NO − 3 -N in autumn (p < 0.05) ( Table 2). In addition, the soil AOA abundance of all treatments was significantly higher than the AOB in both growing seasons (p < 0.05), and the ratio AOA/AOB was 9.16 to 58.35 in spring and 9.63 to 31.50 in autumn (Figure 3c,d). The correlation analysis showed that the AOA abundance was negatively correlated to the moisture and positively correlated to the NO

Soil AOA and AOB Community Alpha and Beta Diversities
The alpha-diversity indices of AOA and AOB communities were observed for Table  3. The sequence numbers of amoA genes in AOA and AOB ranged from 11,467 to 24,485, and the coverage ranged from 99.90% to 99.96%. At the 97% similarity level, AOA and AOB sequences can be divided into 50 to 57 and 51 to 57 OTUs. The average fragment

Soil AOA and AOB Community Alpha and Beta Diversities
The alpha-diversity indices of AOA and AOB communities were observed for Table 3. The sequence numbers of amoA genes in AOA and AOB ranged from 11,467 to 24,485, and the coverage ranged from 99.90% to 99.96%. At the 97% similarity level, AOA and AOB sequences can be divided into 50 to 57 and 51 to 57 OTUs. The average fragment lengths of AOA and AOB were 272 bp and 473 bp, respectively. The AOA Simpson index values of the wheat-cucumber, rye-cucumber, and rape-cucumber treatments were significantly higher than the cucumber monoculture, and the AOB Simpson index value of the rape-cucumber treatment was significantly higher than the cucumber monoculture in spring (p < 0.05). Moreover, the AOA Simpson index values of the wheat-cucumber, chrysanthemum-cucumber, and mustard-cucumber treatments were significantly higher than the cucumber monoculture in autumn (p < 0.05). The AOA Chao 1 index values of the alfalfa-cucumber, trifolium-cucumber, wheat-cucumber, and chrysanthemum-cucumber treatments were significantly higher than the cucumber monoculture in spring (p < 0.05). The AOA Chao 1 index values of all treatments were significantly higher than the cucumber monoculture except for the wheat-cucumber treatment. The AOB Chao 1 index value of the alfalfa-cucumber treatment was significantly higher than the cucumber monoculture in autumn (p < 0.05). The AOA and AOB community structures from spring to autumn were analyzed by NMDS (Figure 4). Non-parametric multivariate statistical analysis showed that the community structures of different treatments were similar (date not shown). However, the soil AOA (ANOSIM, R = 0.1778, p = 0.001; Adonis, R 2 = 0.1418, p = 0.001) and AOB (ANOSIM, R = 0.6588, p = 0.001; Adonis, R 2 = 0.5151, p = 0.001) community structures of all samples were different from spring to autumn.

Soil AOA and AOB Community Composition
The predominant phyla of the AOA community were Crenarchaeota and Thaumarchaeota, occupying more than 75% of the total sequences (Figure 5a). The genus Nitrososphaera (AOA) was only represented 0.31% and 0.61% in two growing seasons. (Figure  5c). The predominant phyla of the AOB community were Proteobacteria and a group of

Soil AOA and AOB Community Composition
The predominant phyla of the AOA community were Crenarchaeota and Thaumarchaeota, occupying more than 75% of the total sequences (Figure 5a). The genus Nitrososphaera (AOA) was only represented 0.31% and 0.61% in two growing seasons. (Figure 5c). The predominant phyla of the AOB community were Proteobacteria and a group of unidentified bacteria found in the environment, occupying more than 95% of the total sequences in spring, while Proteobacteria and, the group of unidentified bacteria found in the environment, and a group of unidentified bacteria occupied more than 90% in autumn (Figure 5b). The predominant genera of the AOB community were the group of bacteria from the norank in the environment, Nitrosospira, group of unidentified from Nitrosomonadales, and Nitrosovibrio, and the known genera accounted for 8.1% (Nitrosovibrio) and 19.8% (Nitrosospira) of the total sequences in spring. The predominant genera of the AOB community were Nitrosomonas, group of bacteria from the norank in the environment, Nitrosospira and, Nitrosovibrio, and the known genera accounted for 6.5% (Nitrosovibrio), 14.2% (Nitrosospira) and, 35.0% (Nitrosomonas) of the total sequences in autumn. Interestingly, Nitrosomonas was the rarest genus in all treatments in spring, while Nitrosomonas increased significantly in autumn. In contrast, the group of bacteria from the norank in the environment decreased significantly from spring to autumn. This indicated that seasonal variation had a significant effect on Nitrosomonas and group of bacteria from the norank in the environment (Figure 5d).

Relationships Between Soil AOA and AOB Communities and Soil Physical and Chemical Properties
The RDA analysis showed that the soil AOA community was significantly affected by moisture (p < 0.05) and NO − 3 -N (p < 0.05), and the soil AOB community was significantly affected by moisture (p < 0.01), pH (p < 0.01), and NO − 3 -N (p <0.01) in spring (Figure 6a,b). The soil AOA community was significantly affected by AK (p < 0.01), and the soil AOB community was significantly affected by soil AP (p <0.01) and AK (p < 0.01) in autumn (Figure 6c,d).

Co-Occurrence Network Analysis
We used the network graph to compare the co-occurrence relationship and community complexity of ammonia-oxidizing bacteria communities between spring and autumn ( Figure 7). The network topological properties of AOA and AOB were shown in Table 4. The results showed that the modularity index of each process was greater than 0.4, indicating that the co-occurrence network of each process was a typical modular structure. There were significant differences in the size of the AOA and AOB community networks in different seasons. The total number of edges, positive correlation edges, and graph density of the AOA network in spring were significantly higher than those in autumn ( Figure  7a,b), In comparison, the total number of edges, positive, and negative correlation edges and graph density of the AOB network in autumn were significantly higher than those in spring (Figure 7c,d). The analysis of the network graph between AOA and AOB showed that the total number of edges, positive correlation edges, and graph density of AOA and AOB in autumn were significantly higher than those in spring. The relationship between AOA and AOB was relatively simple to analyze in spring, but in autumn, the relationship between AOA and AOB was more complex, the AOB community played a leading role in the whole network (Figure 7e,f). These results indicated that seasonal variation regulated the co-occurrence relationship, correlation, and complexity of AOA and AOB, and changed their complex interactions.

Co-Occurrence Network Analysis
We used the network graph to compare the co-occurrence relationship and community complexity of ammonia-oxidizing bacteria communities between spring and autumn ( Figure 7). The network topological properties of AOA and AOB were shown in Table 4. The results showed that the modularity index of each process was greater than 0.4, indicating that the co-occurrence network of each process was a typical modular structure. There were significant differences in the size of the AOA and AOB community networks in different seasons. The total number of edges, positive correlation edges, and graph density of the AOA network in spring were significantly higher than those in autumn (Figure 7a,b), In comparison, the total number of edges, positive, and negative correlation edges and graph density of the AOB network in autumn were significantly higher than those in spring (Figure 7c,d). The analysis of the network graph between AOA and AOB showed that the total number of edges, positive correlation edges, and graph density of AOA and AOB in autumn were significantly higher than those in spring. The relationship between AOA and AOB was relatively simple to analyze in spring, but in autumn, the relationship between AOA and AOB was more complex, the AOB community played a leading role in the whole network (Figure 7e,f). These results indicated that seasonal variation regulated the co-occurrence relationship, correlation, and complexity of AOA and AOB, and changed their complex interactions. and autumn (f). A link represents a significant correlation (Spearman, |ρ| > 0.8, p < 0.05). The blue line represents positive correlation and the red line represents negative correlation. The node size is proportional to the number of connections. A and B after OTUs represent AOA and AOB, respectively. Only the groups with the top 50 average relative abundance of AOA and AOB were shown. Table 4. Topological properties of co-occurrence network of AOA and AOB in spring and autumn. Nodes  50  50  50  50  100  100  Edges  112  97  88  217  228  341  Positive correlation  edge numbers  106  89  81  180  197  283   Negative correlation edge  numbers  6  8  7  37  28

Discussion
Long-term continuous cropping results in excessive soil nutrient consumption and the deterioration of soil physical and chemical properties, while reasonable intercropping can effectively improve soil physical and chemical properties and promote nutrient recycling and transformation in agricultural production [36,37]. Previous studies have shown that intercropping system can affect soil physical and chemical properties [38,39]. Our results showed that compared with the cucumber monoculture, the rye-cucumber treatment significantly increased the soil NH + 4 -N in spring, which could be due to the increase of soil nitrogen content by pruning rye stubble [40]. The alfalfa-cucumber and chrysanthemum-cucumber treatment NO− 3 -N contents were significantly increased in autumn, Figure 7. Co-occurrence network of ammonia oxidizing bacteria communities in different seasons. The AOA community in spring (a) and autumn (b), AOB community in spring (c) and autumn (d), AOA and AOB communities in spring (e) and autumn (f). A link represents a significant correlation (Spearman, |ρ| > 0.8, p < 0.05). The blue line represents positive correlation and the red line represents negative correlation. The node size is proportional to the number of connections. A and B after OTUs represent AOA and AOB, respectively. Only the groups with the top 50 average relative abundance of AOA and AOB were shown. Table 4. Topological properties of co-occurrence network of AOA and AOB in spring and autumn.

Discussion
Long-term continuous cropping results in excessive soil nutrient consumption and the deterioration of soil physical and chemical properties, while reasonable intercropping can effectively improve soil physical and chemical properties and promote nutrient recycling and transformation in agricultural production [36,37]. Previous studies have shown that intercropping system can affect soil physical and chemical properties [38,39]. Our results showed that compared with the cucumber monoculture, the rye-cucumber treatment significantly increased the soil NH +4 -N in spring, which could be due to the increase of soil nitrogen content by pruning rye stubble [40]. The alfalfa-cucumber and chrysanthemumcucumber treatment NO − 3 -N contents were significantly increased in autumn, which could be due to the solid water in autumn that slowed moisture evaporation and reduced the soil NO − 3 -N leaching degree. On the other hand, soil nitrifying bacteria converted the unused NH +4 -N into NO − 3 -N and stored it in the soil under the appropriate soil temperature and pH rang, causing some intercropping systems to increase in NO − 3 -N contents in autumn [41]. In the intercropping systems, the AP content significantly decreased, which is consistent with previous research results [27], this may be due to the increase of plant diversity in intercropping, activate the difficult soluble phosphorus in soil, beneficial to crops on the absorption of phosphorus, thus resulting in a decrease of soil in AP [8], but the intercropping treatments did not significantly change the water, pH and AK (Table 1), this may be caused by the experiment of short-term [42,43].
The ammonia-oxidizing process is mainly conducted by ammonia-oxidizing microorganisms such as AOA and AOB, and PNR is one of the critical indicators of soil nitrification activity. We found that most intercropping treatments significantly reduced soil PNR in spring (Figure 2), which may be due to interspecific competition and synergy. Some studies have found that not all intercropping systems can improve soil ammonia-oxidizing microorganisms. Previous studies have found that compared with monoculture, intercropping of broad bean and maize significantly reduced soil AOA and AOB abundances and rhizosphere microbial abundances. The change of nitrification rate may be related to the competition among species, which may affect the nitrification rate of soil [10,28]. The PNR values of the mustard-cucumber, rape-cucumber, and chrysanthemum-cucumber intercropping treatments were significantly higher than the cucumber monoculture in autumn. This could be justified by the fact that most intercropping treatments have better water retention in autumn, and the soil moisture content of the cucumber monoculture was lower than the intercropping treatments. A previous study found that the soil moisture content was significantly positively correlated with the PNR, and a decreased moisture content can lead to a decrease of PNR [1].
Previous studies have shown that intercropping system can significantly change the abundances of ammonia-oxidizing microorganisms in soil [3,8,9]. Those findings are consistent with the results of this study. This may be related to the improvement of the soil environment by increasing plant diversity [1,2]. In addition, AOB were dominant players in the ammonia-oxidizing process until AOA were discovered and proven to be more abundant and diverse in population compared to AOB in some soils, and it was also found that AOA could play a leading role in nitrification [44,45]. According to the qPCR results, the AOA abundance was significantly higher than AOB in all treatments in the two growing seasons, which was consistent with most findings in agricultural soils [46]. In the comparison of the same treatment, the ratio of AOA to AOB in spring was 9.16 and 58.35, while in autumn, it was 9.63 and 31.50. With the change of seasons, the rate of AOA:AOB showed a downward trend, indicating that AOB rather than AOA may play a leading role in ammonia oxidation [47]. We also found that different growing seasons had a significant effect on AOA and AOB abundances. From spring to autumn, AOB abundance increased in all treatments ( Figure 3). The results of this study are consistent with other results, these studies indicate that seasonal climate change and increased plant diversity caused by intercropping may be the reasons for the change of soil AOA and AOB abundances [48][49][50], indicating that different intercropping patterns and seasonal variation could cause changes in AOA and AOB abundances. However, a study has shown that seasonal heterogeneity had a significant impact on the AOA abundance of the rhizosphere of Phragmites communis but had no significant effect on the AOB abundance [51], suggesting that nitrifying microbial functional genes have different responses to different ecosystems and land management practices.
A variety of environmental factors could influence aerobic ammonia-oxidizing bacteria through complex mechanisms of action, among which pH, moisture, NH +4 -N, and NO − 3 -N have significant effects on the abundances of AOA and AOB [7,15]. In this study, the AOA and AOB abundances were closely related to soil properties. For example, the diffusion capacity of O 2 decreased with an increase in the soil moisture, which significantly affected the microbial activities of AOA [52]. This may be the cause of the decrease in AOA abundance of in soil. Similarly, the factor analysis driving the variation in soil AOA and AOB abundances showed that the moisture and AK were negatively correlated to the AOA abundance, and the pH was negatively correlated to the AOB abundance in spring ( Table 2). The AP was negatively correlation to the AOA abundance in autumn. Studies have shown that the favorable pH for the growth and proliferation of AOB is 7.0-7.40 in soil, while nitrification is inhibited in soils at low pH [53,54]. In our study, the pH in spring ranged from 7.10-7.32, which was higher than that in autumn, ranging from 6.75-6.91. In addition, previous studies also showed that AOA generally dominate nitrification in acidic soils, while AOB dominate nitrification in alkaline soils [55,56], which was contrary to our results. Our results showed that the abundance of AOB decreased with the increase of pH, which may be because a suitable alkaline soil is conducive to the survival of AOB, and excessively high pH reduced AOB activity [57]. Previously, it was shown that only AOB was found to be involved in 13 CO 2 fixation and ammonia-oxidizing processes in the alkaline tidal soils of Fengqiu city in China, and AOB was the dominant microorganisms in the soil nitrification. [58]. Therefore, nitrification may have been dominated by AOB in our study since pH was the key environmental factor regulating the amount of AOB in spring. The reason why AOB cannot grow in a low pH environment may be the low affinity of AOB to substrates, their inability to adapt to an oligotrophic environment and their inability to metabolize in an environment with a pH lower than 6.50 [59]. However, recent studies have found that some strains and species of AOB have also been isolated from acidic agricultural soils with pH values between 5.0-7.50, and the Nitrosospira cluster 3a.2 of AOB isolated from acidic woodland soil could be self-oxidized with ammonia [60,61]. These results showed that the contribution of AOB to ammonia-oxidation in acidic soils could not be ignored and further research to understand its unknown function. AOB could survive well in grass soils with high-dose substrates, while AOA usually survive in soils with low substrate concentrations, proving that AOA are more suitable for growth under oligotrophic conditions and that high substrate concentrations could inhibit their abundance [62]. This could be why the AOA abundance was significantly negatively correlated with the soil AP and AK nutrients in our study.
The soil microbial diversity and ecosystem stability are closely related to nutrient transformation, and intercropping can affect various agro-ecosystem services by improving the soil microecology [36]. Our study showed that the results showed that the diversity index of AOA and AOB was increased by intercropping systems (Table 3), indicating that intercropping improved the diversity of microorganisms involved in the ammonia oxidation reaction in continuous-cropping soil, which may contribute to the enhancement of ammonia oxidation in agro-ecosystems. NMDS and Non-parametric multivariate statistical analysis showed that the three replicates of AOA and AOB were closely distributed, and there was no apparent separation between the treatments, and the community composition of different treatments was similar. However, the AOA and AOB community structures changed significantly from spring to autumn (Figure 4). Studies have shown that changes in seasons had significant effects on the community structures and activity of AOA and AOB in different soil types [10,63]. These results proved that AOA and AOB community structures were susceptible to seasonal variation, and the effect of seasonal variation was greater than intercropping.
The wheat-cucumber and alfalfa-cucumber treatments significantly increased the relative abundance of Nitrososphaera and Nitrosospira for AOA and AOB at the genus level in spring, respectively. From a seasonal perspective, the relative abundance of Nitrosospira and Nitrosovibrio decreased significantly in most treatments from spring to autumn. Interestingly, the most abundant Nitrosomonas in all systems increased significantly, but the group of bacteria from the norank in the environment decreased significantly at the genus level when the season changed from spring to autumn ( Figure 5). This phenomenon may have been caused by the large temperature difference between the two growing seasons in northeast China, which affected the availability of soil substrates and changed the activity of Nitrosomonas and other bacteria; therefore, the AOA and AOB community composition was changed in the two growing seasons [64]. Our results showed that the change in the AOA community was relatively stable compared with the cucumber monoculture, the composition of the AOB community was mainly affected by seasonal variation, and this also further proves that the ammonia-oxidation process of the soil might have been dominated by the AOB community.
Moisture, AP, AK, NO − 3 -N, and pH were related to alterations in the nitrifying bacterial community structure ( Figure 6). The AOA community was mainly affected by the moisture, NO − 3 -N, and AK. The AOB community was affected primarily by the moisture, AP, AK, NO − 3 -N, and pH, which also indicated that AOA were more tolerant to different environmental factors compared to AOB. A previous study has shown that moisture was the main factor affecting AOA and AOB community structures, which is consistent with the results of our study [65], however NO − 3 -N did not cause changes in community structure, which is different from our results [14]. The community structure of AOB was regulated by soil pH. Previous studies have shown that the soil pH is a crucial factor influencing changes in the AOB community structure [14,66]. Notably, the correlation between NH +4 -N and the AOB community was not significant, which is contrary to the results of our study, which indicated that the ammonia concentration had a significant effect on the AOB community [53]. Soil AOA and AOB communities were correlated with the soil nutrients, which may have been caused by the fact that AP and AK changed the proportions of the nutrient matrix in the soil, thus affecting the flora structure of ammonia-oxidizing bacteria, which further indicates that AOA and AOB had a strong dependence on the soil nutrients [67,68].
The complex interactions between species are one of the crucial indicators of community biodiversity [69]. The significant changes observed in the network of two seasons in this study showed that the intraspecific and interspecific interactions of AOA and AOB varied with seasons. The intraspecies interactions of AOA and AOB communities were stronger in spring than in autumn. The results of the interspecific co-occurrence network showed that the interaction (edges) in autumn was significantly stronger than that in spring, indicating that the interaction between AOA and AOB communities was stronger in autumn, and the whole network was dominated by the AOB community ( Figure 7, Table 4). In spring and autumn, the difference of soil ammonia-oxidizing microorganisms and their interaction intensity may be due to the lower soil temperature in spring and the lower activity of ammonia-oxidizing microorganisms in the soil. The relationship between the microbial communities was relatively simple, and the interaction between the microorganisms was less. However, the soil temperature in autumn was higher, and the activity of ammonia-oxidizing microorganisms was higher, which leads to more complex interaction between microbial communities. The results of this study on the interaction between communities caused by seasonal variation are similar to previous studies [70,71]. This result also further proved that with the seasonal variation, the ammonia oxidation in the soil might be dominated by the AOB community.

Conclusions
In summary, this study showed that intercropping significantly changed the soil NH +4 -N and NO − 3 -N contents, and most intercropping systems significantly increased soil PNR, AOA, and AOB abundances. The AOA abundance in all treatments was significantly higher than AOB, and the ratio between them decreased as the seasonal variation. The wheatcucumber and rape-cucumber treatments increased AOA community diversity. Seasonal variation was an important factor affecting the relative abundance of AOB. Seasonal variation had greater influence on the community structures of AOA and AOB than intercropping. Moisture, pH, NO − 3 -N, AP, and AK were essential factors to regulate the AOA and AOB communities. The regulatory effects of AP and AK on the AOA and AOB communities should not be ignored. The interaction between the AOA and AOB communities was simple in spring, but on the contrary in autumn, the AOB community might be the main driver of ammoxidation in autumn.