Early Bolting, Yield, and Quality of Angelica sinensis (Oliv.) Diels Responses to Intercropping Patterns

Intercropping is a sustainable method for cultivating medicinal herbs since it requires lower dependence on chemical fertilizers than a sole cropping system. In this study, we compared the effects of sole cropping and intercropping on early bolting, yield, and the chemical composition of Angelica sinensis (Oliv.) Diels. Field experiments were conducted, in 2018 and in 2019, using different cropping systems including sole cropping of A. sinensis (AS), sole cropping of Vicia faba (VF), and intercropping (without fertilization) at three ratios: one row of A. sinensis + three rows of V. faba, AS/VF (1:3), two rows of A. sinensis + two rows V. faba, AS/VF (2:2), three rows of A. sinensis + one row V. faba, AS/VF (3:1). The effect of each cropping system was evaluated by measuring the dry biomass of V. faba and the dry biomass, ferulic acid content, and essential oil content and composition of A. sinensis. The early bolting rate of A. sinensis was significantly lower in the intercropping system as compared with that in a sole cropping system. The AS/VF (3:1) intercropping pattern resulted in an optimal yield and the highest ferulic acid content of A. sinensis, highest dry biomass of V. faba, and highest land equivalent ratio (LER). Additionally, the A. sinensis was more aggressive (the aggressivity value of A. sinensis was positive, and its competitive ratio was >1) under AS/VF (3:1) intercropping pattern, and it dominated over V. faba (which had negative aggressivity values and a competitive ratio of <1) under AS/VF (3:1) intercropping pattern. Ligustilide was the most dominant component of the essential oil of A. sinensis, regardless of the cropping system; however, the chemical component of essential oil was not influenced by intercropping patterns. Overall, the AS/VF (3:1) intercropping pattern without fertilization was the most productive, with the highest LER and ferulic acid content. These data indicate that intercropping can serve as an alternative for reducing the use of chemical fertilizers and intercropping also decreases the early bolting rate of A. sinensis, thus, enabling its sustainable production.


Introduction
Angelica sinensis (Oliv.) Diels (Apiaceae) is a biennial medicinal herb that grows in shrub meadows, forest margins, hillsides, and river valleys in China [1]. As a famous traditional Chinese medicine, the root of A. sinensis has been used for thousands of years. In the past, A. sinensis has been traditionally prescribed for treating gynecological diseases because of its ability to "replenish and invigorate the blood" [1,2]. In recent years, indepth research on the chemical constituents and pharmacological action of A. sinensis has prompted its use for managing cancer [3], cardiovascular disease [4], and Alzheimer's disease [5]. Meanwhile, it has been used as an ingredient in health care products and cosmetics [6]. Therefore, the increase in application range has led to an increased demand for A. sinensis. In past years, exploitation and habitat destruction have seriously threatened the resources of all natural production areas of A. sinensis. Investigations of wild A. sinensis In this study, we evaluated the effects of intercropping A. sinensis with V. faba on the yield, early bolting rate and essential oil content, and the yield and composition of A. sinensis, and assessed the benefit of intercropping and sole cropping by comparing the LER, RCC, A, and CR values in these cropping systems. We hypothesized that intercropping A. sinensis with V. faba would: (1) provide partial shade to A. sinensis plants, and thus, decreasing the early bolting rate; (2) decrease the yield of A. sinensis because of no fertilizer application, (3) result in higher quality of A. sinensis than sole cropping. There was no significant effect of cropping system on organic matter and K contents both in 2018 and in 2019 ( Table 1). The available N and available P were influenced by the cropping system in 2018. The available N content in AS/VF (2:2) was significantly higher than sole cropping of A. sinensis, while the available P content in AS/VF (3:1) was significantly higher than sole cropping of A. sinensis (Table 1). However, the cropping system had no significant effect on available N and available P in 2019. Meanwhile, there was no significant difference in soil chemical properties between 2018 and 2019 (Table 1). Data are means of three triplicates ± SD. Means in a row with different superscripts differ significantly (p < 0.05, two-way ANOVA, Tukey-HSD). AS, sole cropping of A. sinensis; AS/VF (3:1), three rows of A. sinensis + one row V. faba; AS/VF (2:2), two rows of A. sinensis + two rows V. faba; AS/VF (1:3), one row of A. sinensis + three rows of V. faba.

Correlations among Soil Nutrients and Plant Growth and Essential Oil Production
The results of Spearman's correlation coefficient values among soil organic matter, and available N, P, and K and biomass yield, ferulic acid content, essential oil content, and essential oil yield are presented in Figure 1. On the one hand, the results showed that the content of available K was negatively correlated with biomass yield, essential oil content, essential oil yield, and chloroplast content, both in 2018 and in 2019. On the other hand, in 2018, there was a significant negative correlation between biomass yield, essential oil yield, chloroplast content, and available N, while there was a significant positive correlation between available P and available N ( Figure 1). These results for 2019 were the same as for 2018, except for the significant positive correlation between soil organic carbon and available N. In addition, biomass yield, essential oil content, essential oil yield, and ferulic acid content of A. sinensis increased with an increase in chlorophyll content, both in 2018 and in 2019. Correlations among soil organic matter, available N, P, and K, and biomass yield, ferulic acid content, essential oil content, essential oil yield: (A) in 2018; (B) in 2019 (Note: The criterion for statistical significance was set at p < 0.05; spaces indicate no correlation or insignificant correlation; AK represents available K; AN represents available K; AP represents available K; SOC represents soil organic carbon; BY represents biomass yield of A. sinensis; EOC represents essential oil content of A. sinensis; EOY represents essential oil yield of A. sinensis; FC represents ferulic acid content of A. sinensis; CT represents chloroplast content of A. sinensis.).

Early Bolting Rate of A. sinensis
Intercropping significantly decreased the early bolting rate of A. sinensis as compared with the sole cropping system in 2018 and in 2019 (p < 0.05) ( Table 2). In 2018, the early bolting rate of A. sinensis in the AS/VF (3:1), AS/VF (2:2), and AS/VF (1:3) intercropping patterns were 17.37%, 22.58%, and 26.30%, respectively, which were lower than that of the sole cropping system (36.78%) by 52.77%, 38.60%, and 29.22%, respectively. In 2019, the early bolting followed a similar pattern as in 2018. However, the early bolting rate of A. sinensis in 2018 was slightly higher than that in 2019.

Early Bolting Rate of A. sinensis
Intercropping significantly decreased the early bolting rate of A. sinensis as compared with the sole cropping system in 2018 and in 2019 (p < 0.05) ( Table 2). In 2018, the early bolting rate of A. sinensis in the AS/VF (3:1), AS/VF (2:2), and AS/VF (1:3) intercropping patterns were 17.37%, 22.58%, and 26.30%, respectively, which were lower than that of the sole cropping system (36.78%) by 52.77%, 38.60%, and 29.22%, respectively. In 2019, the early bolting followed a similar pattern as in 2018. However, the early bolting rate of A. sinensis in 2018 was slightly higher than that in 2019.

Chlorophyll Content
The chlorophyll content of V. faba was significantly affected by the intercropping pattern in both years (Table 3). As compared with V. faba monoculture, the AS/VF (3:1) and AS/VF (2:2) intercropping patterns significantly increased the chlorophyll content in the leaves of V. faba in both years. In A. sinensis, although the chlorophyll content was decreased by intercropping, the magnitude of reduction was not significant in both years (p > 0.05) ( Table 3). In addition, there was no significant difference of the chlorophyll contents in leaves of V. faba and A. sinensis between 2018 and 2019.

A. sinensis Root Biomass and V. faba Biomass Yield
The results demonstrated that there was a significant difference in A. sinensis root biomass and V. faba biomass yield under different cropping systems (Table 3). In 2018 and in 2019, the highest biomass yields (5042.8 kg/hm 2 and 5116.5 kg/hm 2 ) of A. sinensis were obtained under sole cropping system. As the number of planting rows of A. sinensis decreased, the biomass yield of A. sinensis also decreased. The lowest biomass yields (525.00 kg/hm 2 and 638.5 kg/hm 2 ) were acquired in the AS/VF (1:3) intercropping pattern in 2018 and in 2019, respectively (Table 3). Intercropping significantly enhanced the V. faba biomass yield. In 2018 and in 2019, the highest (10763.5 kg/hm 2 and 11110.8 kg/hm 2 ) and lowest (4319.6 kg/hm 2 and 4433.3 kg/hm 2 ) V. faba biomass yield were obtained under the AS/VF (3:1) intercropping system and the sole cropping system, respectively ( Table 3). The results also indicated that there was no significant difference in the biomass yield of A. sinensis and V. faba biomass yield between 2018 and 2019.

Ferulic Acid Content
The ferulic acid content of A. sinensis roots was analyzed using HPLC ( Figure S1). The results in relation to the ferulic acid content demonstrated significant differences among cropping systems in both years. On the one hand, the highest ferulic acid contents (0.16% and 0.14%) were recorded under the AS/VF (3:1) intercropping system which demonstrated about 13.22% and 8.07% increase over the sole cropping of A. sinensis in 2018 and in 2019, respectively (Table 4). On the other hand, the lowest ferulic acid contents (0.10% and 0.096%) were found in samples under the AS/VF (1:3) intercropping system in 2018 and in 2019, respectively. Overall, the ferulic acid content had no significant difference in both years.

Essential Oil Content, Yield, and Composition
As compared with the sole cropping system, the essential oil content of A. sinensis was higher in the AS/VF (2:2) intercropping pattern and lower in the AS/VF (3:1) and AS/VF (1:3) patterns, however, the difference in the essential oil content between sole cropping and intercropping systems was not significant (p > 0.05) ( Table 4). Similar results were obtained in both years. These data suggest that intercropping has no significant impact on the essential oil content of A. sinensis. The essential oil yield of A. sinensis was significantly higher in the sole cropping system than in the intercropping system both in 2018 and in 2019 (Table 4), that is, intercropping reduced the essential oil yield of A. sinensis. Meanwhile, there was no significant difference in essential oil and essential oil yield of A. sinensis between 2018 and 2019.
Analysis of the composition of A. sinensis essential oil revealed that ligustilide was the most abundant ingredient, and the content of ligustilide was not affected by intercropping (p > 0.05) in 2018 and 2019 ( Table 5). The amount of senkyunolide A and butylphthalide was lower in the intercropping systems than in the sole cropping system, but the difference was not significant (Table 5). Similar results were obtained in both 2018 and 2019. Meanwhile, as for as senkyunolide H, senkyunolide I, levistolide A and butylidenephthalide, cropping system has no significant influence on them both in 2018 and in 2019 (Table 5). In addition, there was no significant difference in the chemical composition of essential oil between 2018 and 2019.

Competition Indices under Different Intercropping Patterns
The value of LER was greater than one in the AS/VF (3:1) and AS/VF (2:2) intercropping patterns (

Competition Indices under Different Intercropping Patterns
The value of LER was greater than one in the AS/VF (3:1) and AS/VF (2:2)

Early Bolting Rate of A. sinensis
In this study, intercropping significantly decreased the early bolting rate of A. sinensis, as expected. Previous studies have shown that early bolting in A. sinensis was affected by many factors, such as genetics, ecological conditions (light, climate, soil, and water),

Early Bolting Rate of A. sinensis
In this study, intercropping significantly decreased the early bolting rate of A. sinensis, as expected. Previous studies have shown that early bolting in A. sinensis was affected by many factors, such as genetics, ecological conditions (light, climate, soil, and water),

Early Bolting Rate of A. sinensis
In this study, intercropping significantly decreased the early bolting rate of A. sinensis, as expected. Previous studies have shown that early bolting in A. sinensis was affected by many factors, such as genetics, ecological conditions (light, climate, soil, and water), nutrients (organic acid, free amino acids, and soluble sugars), and cultivation practices (sowing date, harvest date, type and method of fertilization, and plant density). Lin (2010) [39] found that early bolting rate decreased by 5% and 5.5% when covered with 75% and 50% sun shade net, respectively. Han (2018) [9] showed that the ratio of organic fertilizer and chemical fertilizer applied to A. sinensis had a significant influence on early bolting; increasing the amount of organic fertilizer and decreasing the amount of chemical fertilizer significantly reduced early bolting. Liu (2003) [24] reported that photoperiod was potentially one of the main reasons affecting early bolting in A. sinensis. Intercropping of V. faba with A. sinensis reduced the natural light intensity through shading (sole crop of A. sinensis, 59.55 klx light intensity, AS/VF (3:1), 42.12 klx; AS/VF (2:2), 44.35 klx; AS/VF (1:3), 47.53 klx), which affected the photoperiod and consequently the early bolting rate of A. sinensis. Therefore, the low early bolting rate of A. sinensis in intercropping patterns observed in this study could be attributed to no fertilization and low light intensity due to shading by V. faba plants.

Dry Root Biomass of A. sinensis
In 2018 and in 2019, the highest dry root yields of A. sinensis were observed in the sole cropping system (5042.8 kg/hm 2 and 5116.5 kg/hm 2 , respectively), followed by the AS/VF (3:1) intercropping pattern (4196.5 kg/hm 2 and 4284.00 kg/hm 2 , respectively). Han (2018) [9] showed that the dry root yield of A. sinensis grown in the Gansu Province ranged from 2608.30 kg/hm 2 (no fertilization) to 6250.00 kg/hm 2 (A2B1C0: organic fertilizer, 36 kg; diammonium phosphate: 1.8 kg; and microbial fertilizer: 0 kg) under different fertilization treatments. Another study showed that the highest dry root yield of A. sinensis (3200 kg/hm 2 ) was obtained in the pure organic fertilizer treatment in Gansu Province [11]. An experiment carried out at the same experimental facility in Qinghai Province showed that the dry root yield of A. sinensis ranged from 5749 to 7287.5 kg/hm 2 under different fertilization treatments (Supplementary Figure S2). In this study, the dry root yield of A. sinensis in the sole cropping system differed from that in the AS/VF (1:3) and AS/VF (2:2) intercropping patterns, and the dry root yield was not affected in the AS/VF (3:1) pattern, as the value was within the range reported previously. Additionally, previous studies have shown that an increase in the leaf chlorophyll content led to an increase in the photosynthetic activity and biomass of medicinal plants [40][41][42]. In this study, although the chlorophyll content of A. sinensis was not affected by the intercropping pattern, the increased trend of biomass was consistent with the increased trend of chlorophyll content, indicating that chlorophyll content was one of the factors affecting the biomass yield of A. sinensis. Therefore, higher A. sinensis biomass yield in the AS/VF (3:1) intercropping pattern may be due to the following reasons: (1) higher chlorophyll content of A. sinensis, leading to higher photosynthetic activity; (2) low early bolting rate; (3) efficient use of resources [43], due to the effect of V. faba on N availability for A. sinensis, or facilitative interaction between V. faba and A. sinensis plants in the intercropping system [44,45]. Moreover, water, radiation, nutrients, and other resources were probably more efficiently used because of differences in the morphological features and growth environments of intercropped plants [46,47]. Furthermore, intra-and interspecific competition was less in the AS/VF (3:1) intercropping pattern than in the sole cropping systems of A. sinensis and V. faba.
As compared with A. sinensis monoculture, the yield of A. sinensis per unit area in intercropping system decreased with a reduction in the proportion of A. sinensis plants.
The results of this study are consistent with those of previous studies on the intercropping of dill/common bean [26], Moldavian/faba bean [27], black cumin/faba bean [28], and faba bean/peppermint [30].

Ferulic Acid Content
According to previous reports, the content of ferulic acid in A. sinensis ranges from 0.0211 to 0.143%, depending on the method of extraction [48][49][50][51][52][53][54]. In this study, the ferulic acid content of A. sinensis varied from 0.095 to 0.160%. In 2018 and in 2019, the highest ferulic acid contents (0.160% and 0.144%) were observed in the AS/VF (3:1) intercropping pattern, which were significantly higher than that obtained in the other two intercropping patterns. Additionally, the ferulic acid content obtained in the AS/VF (3:1) intercropping pattern was greater than the maximum values obtained in previous studies. Thus, the quality of A. sinensis was affected by the intercropping pattern, as the amount of ferulic acid in the AS/VF (3:1) intercropping pattern was higher than that reported previously. Several factors affect the biosynthesis of plant secondary metabolites, such as light, drought, mineral nutrients, and high salinity [55]. In previous studies, the contents of organic carbon (C), total N, total K, available K, calcium (Ca), copper (Cu), zinc (Zn), nickel (Ni), and manganese (Mn) in the soil were positively correlated with the ferulic acid content of A. sinensis [2]. The results in this study were inconsistent with those of previous studies, which showed that the ferulic acid content of A. sinensis was negatively correlated with the available K in the soil (Figure 1). Therefore, the reasons for the high ferulic acid in the AS/VF (3:1) intercropping pattern need further study.

Essential Oil Content, Yield, and Composition
The essential oil contents of A. sinensis roots ranged from 0.45 to 0.96% in this study, which was higher than 0.4% stipulated in the Chinese Pharmacopoeia Commission (2020). Rong et al. (2011) [56] showed that the essential oil content varied from 0.2152 to 1.6231% in different areas. Therefore, the essential oil contents of A. sinensis in different cropping systems in this study were within the range reported previously. The essential oil content of A. sinensis was not influenced by the intercropping pattern, which was supported by the results of Vafadar-Yengeje et al. (2019) [57], who reported no significant difference in the essential oil content of Moldavian balm between its monoculture and Moldavian balm/V. faba intercropping. It must be noted that N plays a vital role in determining the photosynthetic rate and the development and differentiation of essential oil-related cells [58,59]. Therefore, the reason why no difference was observed in the essential oil content of A. sinensis between intercropping and sole cropping systems was probably the absence of exogenous N supply [41].
Our finding that intercropping decreases the essential oil yield of A. sinensis was inconsistent with the results of Machiani et al. (2018) [30], Maffei and Mucciarelli (2003) [60], and Weisani et al. (2015) [26]; all three studies investigated legume/medicinal plant intercropping systems, which enhanced the essential oil yield of the medicinal plant as compared with its sole cropping system. Previous studies have shown that essential oil yield was influenced by crop biomass, biosynthesis rate, oil storage glands, and glandular trichome number and size [40]. In this study, the essential oil content of A. sinensis was not affected by intercropping, suggesting that the essential oil yield depends on the biomass yield of A. sinensis. This result was supported by a correlation analysis; the Spearman's correlation coefficients of essential oil yield and biomass yield were significantly high (Figure 1), which demonstrated that any factor that enhanced the biomass yield could increase the essential oil yield [40]. Higher biomass yield is associated with higher levels of photosynthesis. Therefore, higher chlorophyll content of A. sinensis in the sole cropping system increased its level of photosynthesis, leading to higher essential oil yield [61].
In this study, ligustilide was identified as the principal component of the essential oil of A. sinensis, which was in agreement with previous reports [62][63][64][65]. Additionally, the content of ligustilide was not affected by intercropping (p > 0.05). However, the content of ligustilide in different cropping systems was higher than that in samples reported from previous studies [66,67], but comparable with a recent study using the same method as this study, reported by Gui and Zheng (2018) [65]. Although the contents of senkyunolide H and butylidenephthalide decreased in the intercropping system, their levels were within the range reported previously [65]. Hence, the chemical composition of A. sinensis essential oil was not influenced by the intercropping pattern. This further implied that the essential oil composition of A. sinensis obtained in the intercropping system without fertilization was comparable to that obtained in different fertilization treatments [55,63,65]. Some studies have shown that intercropping legumes with medicinal plants decreased the content of the main essential oil components and increased that of minor essential oil components [26,27], while other studies have indicated that intercropped plants contained a significantly higher content of most of the major essential oil components [28]. In addition, Maffei and Muccia-relli (2003) [60] reported that intercropping systems did not have a significant influence on the chemical profile of the essential oil, consistent with the current study.

Competition Indices
Values of the LER in the AS/VF (3:1) intercropping pattern in 2018 and in 2019 were 1.24 and 1.25, which implied that 24% and 25% more land area would be needed in the sole cropping system to reach the same yield as that obtained in the intercropping system, respectively [36]. These results suggest that intercropping offers a yield advantage over monoculture because of the efficient use of land and other environmental resources (solar radiation, nutrients, and water) for plant growth [15]. Our results were in agreement with the findings of Rostaei et al. (2018) [28], Fallah et al. (2018) [27], and Amani et al. (2018a) [29], who studied dill-faba bean, dragonhead/faba bean, and peppermint/faba bean intercropping, respectively.
Aggressivity was used to express how much higher the increase in A. sinensis production was as compared with V. faba in the intercropping system [33]. Considering all planting patterns, A. sinensis was the dominant species in the AS/VF (3:1) intercropping pattern. CR gives a better evaluation of the competitive ability among intercropping components and is more advantageous with respect to the aggressivity and RCC indices [33]. Generally, the CR index indicates the ratio of partial LER of two plants and takes into consideration the proportion of the plants with which they are initially sown [37]. The CR values of A. sinensis were less than one and higher than that of V. faba in the AS/VF (3:1) intercropping pattern, indicating an absolute yield advantage of A. sinensis over V. faba in intercropping patterns. This was due to the higher aggressivity of A. sinensis and its capability in optimizing the use of available resources as compared with V. faba, leading to a dominant position. Similar results were reported by Yilmaz et al. (2014) [36] in barley and vetch intercropping and by Lithourgidis et al. (2011) [33] in pea/cereal intercropping.

Plant Material
The A. sinensis seeds were obtained from the seed base of the Group for Medicinal Plant Resources and Vegetation Restoration of Tibetan Plateau of Northwest Institute of Plateau Biology, Chinese Academy of Sciences. A specimen with the voucher number Zhou2017115, identified by Professor Guoying Zhou, was collected from the seed base at the flowering stage. Seeds were sown in 2017 and in 2018, according to the standardized technical specifications for A. sinensis. Seedlings were extracted from the soil at 16 weeks after sowing, and then preserved.

Field Trial
The life cycle of A. sinensis is generally completed within 3 years (first year, seedling stage; second year, pharmaceutical period; and third year, bolting period). This study mainly focused on the pharmaceutical period.  Figure 6. The weather data presented in Figure 6 revealed some heterogeneity in the average monthly precipitation between the two growing seasons. In 2018, precipitation was lower in April, May, and June than that in 2019, while precipitation was higher in July, August, September, and October than that in 2019. Soil characteristics, measured prior to the experiment using normal agrochemical procedures, were as follows: ripe soil type, yellow loam; pH = 7.5; organic matter = 12.73 g/kg; total nitrogen (N) = 1.49 g/kg; total phosphorus (P) = 0.83 g/kg; total potassium (K) = 19.27 g/kg; available N = 96.78 mg/kg; available P = 25.59 mg/kg; available K = 160 mg/kg. Prior to the experiment, potato (Solanum tuberosum L.) was grown on the test field. characteristics, measured prior to the experiment using normal agrochemical procedures, were as follows: ripe soil type, yellow loam; pH = 7.5; organic matter = 12.73 g/kg; total nitrogen (N) = 1.49 g/kg; total phosphorus (P) = 0.83 g/kg; total potassium (K) = 19.27 g/kg; available N = 96.78 mg/kg; available P = 25.59 mg/kg; available K = 160 mg/kg. Prior to the experiment, potato (Solanum tuberosum L.) was grown on the test field. The field experiment was carried out in a randomized block design, with one factor (cropping system) and three replications. Five treatments were performed in this experiment: sole cropping of V. faba (VF), sole cropping of A. sinensis (AS), and three intercropping ratios of A. sinensis/V. faba, i.e., one row of A. sinensis + three rows of V. faba, AS/VF (1:3); two rows of A. sinensis + two rows of V. faba, AS/VF (2:2); and three rows of A. sinensis + one row of V. faba, AS/VF (3:1). The soil was plowed to a depth of 30 cm, and the field was divided into blocks and plots (2.6 m × 5 m). The soil was covered with a black plastic sheet (1.2 m) before transplanting A. sinensis and V. faba seedlings. Each plot contained two columns and eight rows per column. A. sinensis seedlings of uniform size were transplanted into the experimental plots on 8 April 2018 and 12 April 2019. In both the sole cropping and intercropping plots, seedlings of A. sinensis were transplanted at 25 cm rowto-row spacing and 20 cm plant-to-plant spacing, whereas V. faba seeds were planted at 25 cm row-to-row spacing and 25 cm plant-to-plant spacing. No chemical fertilizers or pesticides were used during the experiment, and weeds were controlled manually. The experiment was carried out under natural conditions without additional irrigation.

Determination of V. faba Biomass Yield and A. sinensis Biomass Yield
The A. sinensis and V. faba plants were harvested from a 1 m 2 area in each plot (with 100% planting density) at maturity (185 and 154 days after transplanting, respectively) in 2018 and in 2019. The root samples were meticulously washed with tap water, and then dried in the shade at room temperature for a month. The dried root samples of A. sinensis were weighed, and biomass yield (kg/hm 2 ) was calculated. V. faba seeds were threshed and dried at 45 ℃ for 48 h, and V. faba dry biomass yield was calculated per unit area. The field experiment was carried out in a randomized block design, with one factor (cropping system) and three replications. Five treatments were performed in this experiment: sole cropping of V. faba (VF), sole cropping of A. sinensis (AS), and three intercropping ratios of A. sinensis/V. faba, i.e., one row of A. sinensis + three rows of V. faba, AS/VF (1:3); two rows of A. sinensis + two rows of V. faba, AS/VF (2:2); and three rows of A. sinensis + one row of V. faba, AS/VF (3:1). The soil was plowed to a depth of 30 cm, and the field was divided into blocks and plots (2.6 m × 5 m). The soil was covered with a black plastic sheet (1.2 m) before transplanting A. sinensis and V. faba seedlings. Each plot contained two columns and eight rows per column. A. sinensis seedlings of uniform size were transplanted into the experimental plots on 8 April 2018 and 12 April 2019. In both the sole cropping and intercropping plots, seedlings of A. sinensis were transplanted at 25 cm row-to-row spacing and 20 cm plant-to-plant spacing, whereas V. faba seeds were planted at 25 cm row-to-row spacing and 25 cm plant-to-plant spacing. No chemical fertilizers or pesticides were used during the experiment, and weeds were controlled manually. The experiment was carried out under natural conditions without additional irrigation.

Determination of V. faba Biomass Yield and A. sinensis Biomass Yield
The A. sinensis and V. faba plants were harvested from a 1 m 2 area in each plot (with 100% planting density) at maturity (185 and 154 days after transplanting, respectively) in 2018 and in 2019. The root samples were meticulously washed with tap water, and then dried in the shade at room temperature for a month. The dried root samples of A. sinensis were weighed, and biomass yield (kg/hm 2 ) was calculated. V. faba seeds were threshed and dried at 45°C for 48 h, and V. faba dry biomass yield was calculated per unit area.

Chlorophyll Content Analysis
The chlorophyll contents of the A. sinensis and V. faba plants were evaluated using the SPAD meter (SPAD 502, Minolta Ltd., Osaka, Japan). Ten leaves were randomly selected per plot, and the middle section of the leaf blade was used for chlorophyll content measurements; the measurements were carried out at the same time (rapid growing stage) in 2018 and in 2019.

Soil Sampling and Analysis
After harvesting the A. sinensis plants, soil samples (0-30 cm depth) were collected and homogenized from each plot in 2018 and in 2019. Soil samples were sieved through a 2 mm mesh to remove any stones and roots, and then dried at room temperature. The air-dried soil samples were used for analysis to evaluate the chemical properties. The available N, available P, and available K contents were measured using the Kjeldahl method [68], Mo-VF colorimetric method [69], and ammonium acetate flame photometer method [69], respectively. The potassium dichromate oxidation spectrophotometric method was used to measure soil organic carbon [68].

Essential Oil Extraction
The dried roots of the A. sinensis plants were crushed and passed through a 60 mesh sieve. The ground root tissue (50 g) was added to a 1 L round-bottom flask containing 500 mL of distilled water. Then, the flask was connected to a Clevenger-type apparatus, as described by the Chinese Pharmacopoeia Commission (2020). Hydrodistillation was carried out for 5 h in nine replicates (n = 9), and the essential oil content (% (v/w)) and yield (gm −2 ) were calculated using the following equations: Essential oil(%) = Amount of essential oil extracted Amount of ground dried root × 100% Essential oil yield = A. sinensis biomass yield × Percent essential oil content

Sample Preparation
Fresh root samples of A. sinensis (1 g) were ground in liquid nitrogen. Then, for each sample, 1 mL of 70% methanol was added to the ground sample and vortexed for 3 min. The sample was shaken every six times at 30 min intervals and stored at 4 • C overnight. After overnight incubation, the sample was centrifuged at 12,000 rpm for 15 min, and 300 µL of the supernatant was filtered through a 0. The supernatant was separated using an analytical column (Waters ACQUITY UPLC HSS T3, USA), with a particle size of 1.8 µm and column size of 3.0 mm × 100 mm. The column temperate was maintained at 40 • C during separation. Solvent A (0.1% formic acid in water) and solvent B (0.1% formic acid in acetonitrile) were used as the mobile phases. A 3 µL injection volume was used, and the flow rate was maintained at 0.5 mL/min. Sample separation was carried out using the following gradient conditions: 0.5 min, 10% B; 5 min, 90% B; 5.1 min, 10% B; 7 min stop.
The mass spectra were acquired using the QTRAP4500 system with a Duo Spray source in positive and negative ESI modes using the following parameters: ion spray voltage, −4500 V (negative ion mode) and 5500 V (positive ion mode); collision gas, medium; Turbo V spray temperature, 600 • C; heater gas (Gas 2), 60 psi; nebulizer gas (Gas 1), 50 psi; curtain gas, 20 psi. Data was acquired using the multiple reaction monitoring (MRM) mode, and the retention time corresponding to each MRM transition was automatically set using the MassLynx software. The Peak View Software 2.2 (SCIEX, Foster City, CA, USA) was used for data analysis.
4.9. Ferulic Acid Analysis 4.9.1. Sample Preparation Ground A. sinensis roots (200 mg) were dissolved in 20 mL of 70% methanol, and the samples were weighed. The conical flasks were heated to reflux at 80 • C for 3 min. Then, the samples were cooled and weighed again, and 70% methanol was added to the flasks to compensate for the lost weight. After shaking the samples and allowing them to stand still, the supernatant was collected and filtered. The filtrate (10 mL) was passed through a 0.45 µm organic membrane. Ferulic acid in each sample was analyzed by high performance liquid chromatography (HPLC).

HPLC Analysis
To conduct the HPLC analysis, standard ferulic acid (Batch number 17092501) was purchased from Chengdu Pufeide Biotechnology Co. Ltd., which was diluted to a concentration of 0.012 mg/mL using methanol. Quantitative and qualitative analyses of ferulic acid were performed using the Agilent 1260 (Agilent USA) Infinity II Quaternary Pump (G7111A). Data acquisition was performed using the Agilent HPLC software. Ferulic acid was separated on a unitary C18 chromatographic column (4.6 mm × 250 mm, 5 µm). Octadecyl silane bonded silica gel was used as a filler. The HPLC analysis was conducted using acetonitrile: 0.085%/phosphoric acid (20:80) solution as the mobile phase, with the following parameters: detection wavelength = 316 nm, column temperature = 35 • C, flow rate = 1 mL/min, injection volume = 10 µL. The calibration curve was constructed using serial dilutions of ferulic acid standard in methanol, with concentration ranging from 0.00324 to 0.054 mg/mL. The LER index was adopted to compare the A. sinensis/V. faba intercropping system with the sole cropping system of each crop. The value of LER was calculated according to the following equation: where Yvm and Yam represent the yield of V. faba and A. sinensis, respectively, in the sole cropping system, and Yvi and Yai represent the yield of V. faba and A. sinensis, respectively, in the intercropping system. LER > 1.0 indicates a yield advantage of intercropping over monoculture, whereas LER < 1.0 indicates no yield advantage of intercropping over monoculture [70,71].

Aggressivity (A)
The competitive relationship between two crops was evaluated by the aggressivity (A) as proposed by Willey (1979) [72] using the following formula: where Zai is the sown proportion of A. sinensis in intercropping with V. faba; and Zvi is the sown proportion of V. faba in intercropping with A. sinensis.

Crowding Coefficient (RCC)
In the intercropping system, the relative dominance of one crop over the other was calculated with the crowding coefficient based on the following equation: When the values of the two coefficients (RCCv and RCCa) are >1, there is a yield advantage in intercropping as compared with monoculture; when they are equal to 1, there is no yield advantage, and when they are <1, there is a disadvantage in yield.

Competitive Ratio (CR)
Another index, the competitive ratio (CR), was used to evaluate competitive ability among intercropping components. This index was more advantageous than other indices due to the fact that the CR gives stronger competitive ability to the species. The CR implies the ratio of the individual LER of the intercrop component in which they were initially sown proportionally. We used the following equation suggested by Dhima (2007) [37] to calculate the CR index: where CRv is the crowding ratio of V. faba and CRa is the crowding ratio of A. sinensis.

Statistical Analysis
All statistical analyses were performed using SPSS 21.0 (SPSS inc., Chicago, IL, USA). Prior to all statistical analyses, we tested the heterogeneity of variances using Levene's test, and the original data were normalized using log transformation or standardization prior to statistical analysis when necessary. A two-AVOVA with post hoc tests was used to determine the effect of treatment, year, and treatment × year, on all variables, including soil properties (organic matter, available nitrogen, available phosphorus, and available potassium), early bolting rate, chlorophyll content, biomass, ferulic acid content, essential oil content, essential oil yield, and essential oil composition. Then, the R software (version 4.1.0) was used for the data analysis and production of figures. We used the packages "corrplot" and "psych" in the R software to calculate correlation coefficient and p-value of significance test.

Conclusions
Based on the outcomes of the field cultivation of A. sinensis using an intercropping system without fertilization, we reached the following conclusions as described below: (1) Intercropping significantly reduces the early bolting rate of A. sinensis.
(2) The AS/VF (3:1) intercropping pattern results in an optimal yield and the highest ferulic acid content of A. sinensis, highest biomass yield of V. faba, and highest land equivalent ratio (LER), thus, the AS/VF (3:1) intercropping pattern without fertilization is the most productive with high quality. These data indicate that intercropping can serve as an alternative for reducing the use of chemical fertilizers and can also decrease the early bolting rate of A. sinensis, thus, enabling its sustainable production.
Author Contributions: L.Y., data curation and writing-original draft preparation; J.L., investigation and methodology; Y.X., conceptualization, methodology, and software; G.Z., writing-reviewing and editing. All authors have read and agreed to the published version of the manuscript. Data Availability Statement: Not applicable.