In Nitrate-Rich Soil, Fallopia x bohemica Modiﬁes Functioning of N Cycle Compared to Native Monocultures

: The e ﬀ ects of invasive species at the ecosystem level remain an important component required to assess their impacts. Here, we conducted an experimental study with labeled nitrogen in two types of soil (low and high nitrate conditions), investigating the e ﬀ ects of (1) the presence of Fallopia x bohemica on the traits of three native species ( Humulus lupulus , Sambucus ebulus , and Urtica dioica ) and (2) interspeciﬁc competition (monoculture of the invasive species, monoculture of the native species, and a mixture of invasive / native species) on nitriﬁcation, denitriﬁcation, and related microbial communities (i.e., functional gene abundances). We found that the species with the higher nitrate assimilation rate ( U. dioica ) was a ﬀ ected di ﬀ erently by the invasive species, with no e ﬀ ect or even an increase in aboveground biomass and number of leaves. F. x bohemica also decreased denitriﬁcation, but only in the soil with high nitrate concentrations. The impacts of the invasive species on nitriﬁcation and soil microorganisms depended on the native species and the soil type, suggesting that competition for nitrogen between plants and between plants and microorganisms is highly dependent on species traits and environmental conditions. This research highlights that studies looking at the impacts of invasive species on ecosystems should consider the plant–soil–microorganism complex as a whole. one-way on ﬀ plant on two Kruskal–Wallis tests. n 5 for and n 15 for The ﬀ as the percent in the relative this ﬀ


Introduction
Because an invasive species is a new element in a system, it has an impact (positive or negative) on this system simply by being present [1]. Beyond this observation, there is no consistency in research on the ecological impacts of invasive species at the ecosystem level [1], although several reviews and meta-analyses have tried to clarify patterns. Plant invasions can alter various components of Diversity 2020, 12 ecosystems, such as carbon and nitrogen (N) cycles [2][3][4]. The majority of studies suggest increased mineral and nutrient contents in soils [5]. A meta-analysis including studies on 130 species showed that invasive plants enhanced microbial activity and available N in soil (e.g., ammoniac or ammonium) and that the effects on nitrification (ammonium-NH 4 + oxidation to nitrate-NO 3 − ) were highly variable [6].
Higher rates of net N mineralization (conversion of organic N to inorganic N) and nitrification were reported by literature [7], and an inhibition of nitrification was also shown [8,9]. More recently, it has been shown that some of these invasive plants might play a role in denitrification (NO 3 − reduction into gaseous forms [10,11]). All of these studies suggest that invasive plants might have positive or negative feedbacks on soil microorganism functioning, facilitating their establishment in invaded ecosystems, and impacting plant-plant competition for N-substrates [12].
In 2014, Castro-Díez et al. [13] also found an increase in N pools and fluxes that were explained by climate and partly by functional differences between invasive and native species. Plant community composition is theoretically one driver affecting how invasive species will impact resident biota [14]. An invasive species that uses the same resources as the native species but uses these resources more efficiently will result in a scenario of competitive exclusion and is likely to directly impact the microbial community [1,15]. In contrast, an invasive plant species that uses different soil resources than those used by the native species can coexist with other plant species [15] but may have indirect impacts on the plant community by affecting soil properties and communities in a new way [1]. Finally, because invasive plant species can substantially alter belowground communities and nutrient cycling and availability and because the magnitude and direction of these impacts are difficult to predict due to the inherent complexity of ecosystems, there is a need for more studies integrating multiple trophic levels [14,16].
Japanese knotweeds s.l. are a tall herbaceous invasive species that form a polyploid hybrid complex including F. japonica, F. sachalinensis, and F. x bohemica [17]. The latter is the hybrid taxon, phenotypically and genetically more diverse than the two others, and it often shows higher invasiveness [18]. These species are listed as some of the "worst" invasive species based on their impacts. In Nentwig et al.'s list of the worst alien species in Europe [19], which was based on a generic impact scoring system (GISS), Fallopia japonica is ranked 37th out of 86 total species and 12th among plant species. Several studies (reviewed in [20]) have shown that Japanese knotweeds s.l. have major impacts on ecosystems, such as alterations of soil chemistry, leading to increases or decreases in the abundance of soil microorganisms (amoebas or bacteria, respectively), and a disturbance to nutrient cycling [10,11]. These effects have been proven to depend on soil characteristics [11,21]. However, to our knowledge, nothing is known about the influence of the presence of native species co-occurring with F. x bohemica, with potentially different nitrogen-use strategies, on soil N functioning (e.g., nitrification and denitrification). However, there is increasing evidence in the literature that plant nitrogen-use strategies influence the soil microbial community involved in the N cycle and their functioning [22][23][24]. In two contrasting soils with different nitrate levels, we aimed to identify (1) the impacts of F. x bohemica in mixture with native species on soil N functioning and associated microbial communities and (2) whether these impacts are driven by native/invasive plant interactions and by different N-use strategies for native species compared to Japanese knotweeds.
We hypothesized that F. x bohemica would (1) inhibit denitrification and (2) increase or not affect nitrification to maintain nitrate available in the soil for its growth, especially in low resource environments (following the niche construction hypothesis), and where the competitive effect of F. x bohemica on native plant species is the strongest.

Study Species, Soils, and Plant Regeneration
Three native species were chosen. Clonal (with a belowground perennial organ) and herbaceous species commonly occurring in the same habitats as F. x bohemica (F) were selected: Urtica dioica (U), Sambucus ebulus (S), and Humulus lupulus (H). Individuals from at least 6 populations of each species, separated by at least 3 km, were collected in the Auvergne-Rhône-Alpes region in France (Table S1) and were cultivated in a common garden for approximately 18 months. Two 2 mm-sieved standardized sandy loam soils were purchased at LUFA Speyer (Speyer, Germany) and chosen based on their physicochemical and textural similarities. The two soils have contrasting soil nitrate (NO 3 − ) concentrations containing 20 ± 1.6 µg NO 3 − ·g −1 dry soil and 9 ± 0.8 µg NO 3 − ·g -1 dry soil (hereafter named N + soil and N − soil, respectively). The soil ammonium concentrations were very low for the two soils (0.09 ± 0.02 and 0.04 ± 0.02 µg NH 4 + ·g −1 dry soil).
More than 100 pieces of rhizome or stolon per species and per soil were planted in 0.5 L pots (one individual per pot), which were filled with one of the two soils, and cultivated in controlled conditions. Because the species showed different regeneration times (U. dioica > S. ebulus > H. lupulus > F. x bohemica; measured before the experiment using the same plant material cultivated in common gardens, data not published), an interval of several days was allowed between the planting of the species based on their regeneration times. This interval prevented the priority effect from the faster regenerating species and enabled the study to focus on the competition of resource use.

Experimental Design
Ten weeks after the regeneration time, the height of the regenerated plants was measured. The individuals selected for this experiment were near the height mean of each species, regardless of the population, to reduce intraspecific size variability. For each soil, 1.5 L pots containing two individuals separated by 5 cm were constituted. Seven competition treatments were created: U-U, F-U, S-S, F-S, H-H, F-H, and F-F, with 8 pots for each combination type of soil × competition treatment (i.e., 8 repetitions, n = 112). Plants were grown for 6 weeks at 20 • C with 16 h light/8 h dark. Each pot was watered three times a week. After 6 weeks of plant growth, a 15 N labeled N solution was added to five replicates of each treatment to follow the fate of nitrate in soil and plants, while the other three replicates received unlabeled N solution to measure the natural isotopic composition in the soil and plants ( 15 N at% contr ).

N Addition and Labeling
To avoid any fertilizing effect in our experiment, the amount of added 15

Soil and Plant Harvest and Processing
The plant and soil samples were collected 9 hours after the N labeling treatment. To avoid any bias due to NO 3 − transformation caused by microbial activity in the soil between the first and last pot collections, all pots were kept at 4 • C before their harvest. The aboveground parts and the roots of each individual were oven-dried at 70 • C for 48 h and crushed with a ball mill mixer (Retch MM200 Eragny, France) for the 8 replicates. The 5 replicates with 15 N addition were used (1) for plants, to determine the plant traits and the 15 N concentrations in plants and (2) for rhizospheric soil (i.e., the soil under root influence), to analyze the microbial N activities and the microbial N abundances. The 3 replicates for the unlabeled pots were only used for the natural soil and plant 15 N values. For plants, individual aboveground and root parts were separately collected from each pot. The rhizospheric soil, considered as the remaining soil adhering to the roots after a manually shaking [25], was carefully sampled and was used for 15 N, soil (i.e., nitrate and ammonium concentrations) and microbial analysis (i.e., microbial activities and gene abundances). Nitrification and denitrification Diversity 2020, 12, 156 4 of 15 enzyme activity, as well as the nitrate and ammonium concentrations, were measured on fresh rhizospheric soil, and a subsample was immediately frozen at −20 • C for the DNA extractions and gene abundance analysis.

15 N Analysis
Two mg of plant material was weighed to the nearest 0.1 mg in tin capsules (i.e., the aboveground parts and the roots of each individual separately). The nitrogen concentration and isotope composition were measured using an elemental analyzer (Flash EA 1112, Thermo Scientific) coupled in a continuous flow with an isotope ratio mass spectrometer (Isoprime 100, Elementar, Lyon, France). Data were normalized using IAEA N-2, IAEA 305A and B international reference materials and expressed as 15 N atom% ( 15 N at% = 15 N / ( 14 N + 15 N) × 100) [26]. 15 N uptake (A 15 N in µg/g, hereafter called "nitrate assimilation") was calculated using:

Plant Traits
We considered some plant functional traits, defined as "the characteristics of an organism that are considered relevant to its response to the environment and/or its effects on ecosystem functioning" [27]. At the end of the experiment, the plant height of all individuals was measured, and the leaves were numbered. Before crushing for 15 N analysis, the dry aboveground and root parts of the plants were weighed to obtain aboveground biomass (AB) and root biomass (RB). In our study, the root system was preferred to the belowground system (i.e., the whole system with the rhizome, their buds, and the roots), because the root system is more relevant to assess the impacts on the soil nitrogen cycle and the associated microbial N communities [28].

Nitrification Enzyme Activity (NEA)
To measure nitrification enzyme activity (NEA), 3 g of equivalent dry soil was incubated with 6 mL of (NH 4 ) 2 SO 4 solution (50 µg N-NH 4 + ·g −1 dry soil), and distilled water was added to each sample to achieve 24 mL of total liquid volume in a flask [10]. The flasks were placed in a stirring oven at 28 • C and 140 rpm. One milliliter of soil slurry was sampled after 2, 4, 6, 8, and 10 h of incubation, filtered at 0.2 µm and stored in vials at −20 • C until the NO 2 − and NO 3 − concentrations were measured using an ion chromatograph (ICS 900, Thermo Scientific™ Dionex™, Salt Lake City, UT, USA). Nitrification was expressed in µg of N-NO 2 − + N-NO 3 − ·h −1 ·g −1 of dry soil. To measure denitrification enzyme activity (DEA), 10 g (equivalent dry weight) of fresh soil was placed in a flask hermetically sealed with a rubber stopper [29]. The atmosphere of the flasks was replaced by a 90% helium and 10% acetylene mixture to obtain anaerobic conditions and inhibit the nitrous oxide reductase that catalyzes the final step of denitrification, in which N 2 O is converted into N 2 . Distilled water (2 mL) containing KNO 3 (50 µg N-NO 3 − .g −1 dry soil), glucose (500 µg C-glucose.g −1 dry soil), and glutamic acid (500 µg C-glutamic acid g −1 dry soil) was added through the rubber stopper using a syringe to ensure non-limiting amounts of carbon and NO 3 − for denitrification activity. The flasks were incubated at 28 • C for 8 h. After 2 h, N 2 O accumulation in the flask was measured every hour using a gas chromatograph (µGC R3000, SRA Instrument, Marcy L'Etoile, France). Denitrification was expressed in µg of N-N 2 O g −1 dry soil h −1 .

Soil Nitrate and Ammonium Concentrations
Soil NO 3 − and NH 4 + concentrations were determined following extraction from soil subsamples with a solution of calcium chloride (CaCl 2 ) at 0.01 M, and 20 mL of calcium chloride was added to 5 g of equivalent dry mass soil. The soil was shaken at 140 rpm for 2 h at 10 • C. The suspensions were filtered (Acrodisc Supor ® PF with membrane prefilter 0.8/0.2 µm, VWR) and stored at −20 • C before the amounts of NO 3 − and NH 4 + were measured by ionic chromatography described in Section 2.7.1.

Microbial Abundance Measurements
Soil DNA was extracted from 0.5 g of frozen soil using the FastDNA R Spin Kit for Soil (MPbio, Irvine, CA, USA) following the manufacturer's instructions. DNA concentrations were determined with the Quant-iT™ PicoGreen dsDNA Assay Kit (Invitrogen, Carlsbad, CA, USA). The abundances of the four microbial groups (ammonium-oxidizing bacteria and archaea and nitrite reductase genes) were determined by quantitative PCR (qPCR), targeting group-specific genes encoding key enzymes of the denitrification and nitrification pathways.
Abundances of ammonium-oxidizing bacteria (AOB) and archaea (AOA) were quantified based on the gene copy numbers present per gram of dry soil of the ammonium monooxygenase gene amoA (for both bacteria and archaea). The amoA functional gene was amplified using the gene primers amoA_1F and amoA_2R for the AOB and CrenamoA616r and CrenamoA23f for the AOA, [30]. The abundances of denitrifying bacteria were measured based on the abundance of the nitrite reductase genes nirS and nirK. Fragments of the nirK and nirS genes encoding the copper and cd1 nitrite reductases, respectively, were amplified as described by [31]. Duplicate quantitative real-time PCR (qPCR) was carried out using SYBR green as the fluorescent dye on a Lightcycler 480 (Roche Diagnostics, Meylan, France). The efficiencies of all performed qPCRs were in the range of 90-99%. The R 2 values were always greater than 0.95. The possible inhibition of PCR was tested in advance, and appropriate dilutions were chosen (data not shown).

Statistical Analyses
For all analyses, we excluded the 3 replicates for the unlabeled pots and we used the 5 replicates for the 15 N labeled pots. We analyzed our data separately for the two types of soil to assess the impact of F. x bohemica on species and on the functioning depending on plant species. Differences in aboveground biomass, root biomass, nitrate assimilation and nitrate accumulation between species were analyzed with one-way ANOVAs on monocultures for both soil types. Effects of plant type (natives, invasive, and mixtures) on plant biomass, nitrification, and denitrification in the two types of soil were analyzed with Kruskal-Wallis tests. n = 5 samples for invasive and n = 15 samples for mixtures and native species. The effects of F. x bohemica on traits of H. lupulus, S. ebulus, or U. dioica are presented as the percent change in the mixture relative to the monoculture condition of this species (H-H, S-S, and U-U, respectively). Significant differences to 0 (i.e., 0%, no effect of F on native species in mixture) were analyzed by Student's t-tests and differences between native species were analyzed by one-way ANOVAs in each soil. The effects of FF-FH-HH, FF-FS-SS, and FF-FU-UU on microbial activities and functional gene abundance were also tested by one-way ANOVA. Significant differences between means were tested by Tukey's multiple comparison tests (p < 0.05). The normality of residuals and homogeneity of variance were tested by Shapiro and Levene's tests, respectively. When necessary, data for microbial activities were squared root-transformed. Statistical analyses were performed with Jump Pro ® software (v.13.2).
Other results beyond the scope of our article, i.e., a comparison of soil types on plant traits, or the effect of native species on F. x bohemica, are shown without comments in the supplementary materials (Table S2 and Figures S1 and S2).

Comparison between Species for Growth and N Management in Plants
In low nitrate conditions, the aboveground biomass of monocultures were not different (p = 0.255), while in high nitrate conditions, F obtained higher aboveground biomass than the three other species ( Figure 1A, p < 0.001). For roots, both in low (p = 0.009) and in high (p = 0.039) nitrate conditions, U showed the highest biomass and S showed the lowest ( Figure 1B). In the NO 3 − poor environment, the nitrate assimilation rate was not different between species (p = 0.163). In the NO 3 − rich environment, U showed the highest assimilation rate and S showed the lowest ( Figure 1C, p = 0.034). The accumulation of labeled nitrate was higher in U in NO 3 − poor soil (p = 0.008) and in F compared to S in NO 3 − rich soil ( Figure 1D).
Diversity 2020, 12, x FOR PEER REVIEW 6 of 16 In low nitrate conditions, the aboveground biomass of monocultures were not different (p = 0.255), while in high nitrate conditions, F obtained higher aboveground biomass than the three other species ( Figure 1A, p < 0.001). For roots, both in low (p = 0.009) and in high (p = 0.039) nitrate conditions, U showed the highest biomass and S showed the lowest ( Figure 1B). In the NO3 -poor environment, the nitrate assimilation rate was not different between species (p = 0.163). In the NO3rich environment, U showed the highest assimilation rate and S showed the lowest ( Figure 1C, p = 0.034). The accumulation of labeled nitrate was higher in U in NO3 -poor soil (p = 0.008) and in F compared to S in NO3 -rich soil ( Figure 1D).     biomass showed a large increase in mixture with F compared to S-FS and U-FU mixtures (p = 0.0095; Figure 2D). For morphological traits, no significant differences in the effect of F. x bohemica on height were found between species (Figure 2A), while for leaf number, in NO3 -rich soil, U was found to have more leaves in the mixture compared to its monoculture (p = 0.05; Figure 2B).

Impacts on Plant N Uptakes
No significant differences were found between the responses of each native species to their mixtures with F in terms of assimilation rates and nitrate accumulation ( Figure 3A,B). However, we observed a significant effect of the presence of F. x bohemica on species nitrate assimilation and accumulation compared to their monocultures (i.e., when the percentage is significantly different from 0). While nitrate assimilation of H. lupulus was not impacted compared to the other natives, its nitrate accumulation decreased in low nitrate conditions (p = 0.0413, −50%), following its decrease in biomass (Figure 2). In low nitrate conditions, U and S were found to have a significant decrease for assimilation rates (p = 0.019, −26% on average and p < 0.0001, −44%, respectively for U and S). For NO3 -rich soil, only significant effects were found in nitrate accumulation for S (p = 0.0006, −88%). While H and S were negatively affected compared to their monocultures, U showed higher aboveground biomass in both NO 3 − poor and NO 3 − rich soils compared to other native species (p = 0.018 and p = 0.05 for low and high nitrate conditions compared to H and S, respectively; Figure 2C). For root biomass, significant changes between native species were found only for NO 3 − rich soil. H root biomass showed a large increase in mixture with F compared to S-FS and U-FU mixtures (p = 0.0095; Figure 2D). For morphological traits, no significant differences in the effect of F. x bohemica on height were found between species (Figure 2A), while for leaf number, in NO 3 − rich soil, U was found to have more leaves in the mixture compared to its monoculture (p = 0.05; Figure 2B).

Impacts on Plant N Uptakes
No significant differences were found between the responses of each native species to their mixtures with F in terms of assimilation rates and nitrate accumulation ( Figure 3A,B). However, we observed a significant effect of the presence of F. x bohemica on species nitrate assimilation and accumulation compared to their monocultures (i.e., when the percentage is significantly different from 0). While nitrate assimilation of H. lupulus was not impacted compared to the other natives, its nitrate accumulation decreased in low nitrate conditions (p = 0.0413, −50%), following its decrease in biomass ( Figure 2). In low nitrate conditions, U and S were found to have a significant decrease for assimilation rates (p = 0.019, −26% on average and p < 0.0001, −44%, respectively for U and S). For NO 3 − rich soil, only significant effects were found in nitrate accumulation for S (p = 0.0006, −88%).

Impacts on Plant N Uptakes
No significant differences were found between the responses of each native species to their mixtures with F in terms of assimilation rates and nitrate accumulation ( Figure 3A,B). However, we observed a significant effect of the presence of F. x bohemica on species nitrate assimilation and accumulation compared to their monocultures (i.e., when the percentage is significantly different from 0). While nitrate assimilation of H. lupulus was not impacted compared to the other natives, its nitrate accumulation decreased in low nitrate conditions (p = 0.0413, −50%), following its decrease in biomass ( Figure 2). In low nitrate conditions, U and S were found to have a significant decrease for assimilation rates (p = 0.019, −26% on average and p < 0.0001, −44%, respectively for U and S). For NO3 -rich soil, only significant effects were found in nitrate accumulation for S (p = 0.0006, −88%).

Effects on Plant Biomass
With low NO 3 − concentrations, no change in root biomass was found between invasive, native, and mixed cultures (Table 1a). With high nitrate conditions, significant differences were found (1) between natives and mixtures for root biomass with highest root biomass for invasive and mixtures, and (2) in aboveground parts (p < 0.0001, Table 1b) with the lowest aboveground biomass for native species.

Effects on Nitrification and Denitrification
Plant type had significant effects on NEA, although these effects varied with soil type ( Figure 4A). While NEA increased in the mixture compared to that in the native species in soil with low NO 3 − Diversity 2020, 12, 156 9 of 15 concentrations (X 2 = 6.25; p = 0.04), NEA decreased in the mixture and native species groups compared to that in the invasive species in soil with high NO 3 − concentrations (X 2 = 8.12; p = 0.02). We found significant negative effects on the DEA of invasive species monocultures or mixtures compared to the native species monocultures ( Figure 4B) but only in soil with high NO 3 − concentrations (X 2 = 16.9; p = 0.0002).

Effects on Nitrification and Denitrification
Plant type had significant effects on NEA, although these effects varied with soil type ( Figure  4A). While NEA increased in the mixture compared to that in the native species in soil with low NO3concentrations (Χ² = 6.25; p = 0.04), NEA decreased in the mixture and native species groups compared to that in the invasive species in soil with high NO3 -concentrations (Χ² = 8.12; p = 0.02). We found significant negative effects on the DEA of invasive species monocultures or mixtures compared to the native species monocultures ( Figure 4B) but only in soil with high NO3 -concentrations (Χ² = 16.9; p = 0.0002).

Effects of Competition Treatments on Soil Functioning (Nitrification and Denitrification) in Terms of Microbial Activities and Gene Abundances
Only a few effects on soil functioning were found in soils with low NO 3 − concentrations. No change in nitrification or denitrification was found with the comparison of F and H activities ( Figure 5A,D) and gene abundances (Table 2), and the same results were found for the comparison between F and U ( Figure 5C,F and Table 2). However, significant effects were found on NEA in the mixture of F and S (p = 0.025; Figure 5B), AOA (p = 0.023, Table 2), nirK (p = 0.033; Table 2) and nirS (p = 0.004; Table 2). Soil with F-F and S-S had a significantly lower nitrification rate than that of the soil planted with F-S (p = 0.036 and p = 0.051, respectively, Figure 5B). The same pattern was found for AOA (p = 0.029 and p = 0.056 for comparisons of F-S vs. F-F and F-S vs. S-S, respectively) and for nirS (p = 0.005 and p = 0.013 for F-S vs. F-F and F-S vs. S-S, respectively). For nirK, gene abundance was significantly higher under F-S than under S-S (p = 0.03). Stronger effects on soil functioning were found in soils with high nitrate concentrations. NEA was reduced when the invasive species competed with two native species ( Figure 5A, F-F and F-H comparison, p = 0.0018 and Figure 5C, F-F and F-U comparison, p = 0.005), but this pattern was not observed for S ( Figure 5B). The results on the nitrifying gene abundances showed the same pattern as NEA only for U, where AOA and AOB were more abundant in F-F than in F-U (p = 0.011 and p = 0.01, respectively). In this case, NEA and the associated functional genes were significantly correlated (linear regression between NEA and AOA, R 2 = 0.38, p = 0.014 and linear regression between NEA and AOB, R 2 = 0.34, p = 0.023). The AOB was found to be significantly lower in the S-S treatment than in the F-F treatment ( Table 2, p = 0.013). For DEA, only S-S and U-U showed a tendency to increase compared to F-F ( Figure 5E,F; p = 0.078 and p = 0.065, respectively). For H, DEA was higher in H-H than in F-H ( Figure 5D, p = 0.008). No significant effect was found on denitrifying genes in the comparison between invasive and native species. 0.01, respectively). In this case, NEA and the associated functional genes were significantly correlated (linear regression between NEA and AOA, R² = 0.38, p = 0.014 and linear regression between NEA and AOB, R² = 0.34, p =0.023). The AOB was found to be significantly lower in the S-S treatment than in the F-F treatment ( Table 2, p = 0.013). For DEA, only S-S and U-U showed a tendency to increase compared to F-F ( Figure 5E and 5F; p = 0.078 and p = 0.065, respectively). For H, DEA was higher in H-H than in F-H ( Figure 5D, p = 0.008). No significant effect was found on denitrifying genes in the comparison between invasive and native species.   Table 2. Effects of competition treatments on nitrifying (copy numbers of ammonia-oxidizing bacteria (AOB) and archaea (AOA)) and denitrifying (copy numbers of Cu nitrite reductors nirK and nirS) gene abundances for two soils (a: with low nitrate concentrations and b: high nitrate concentrations). Mean and standard errors are shown; n = 5. The results of one-way ANOVA are presented with p-values when p < 0.05 (ns: non-significant).

Impact of Knotweed on Competing Species: From Little Negative Effect to Facilitation
The comparison of monocultures of the four species showed contrasting strategies between the species, especially between F. x bohemica and U. dioica. More than the other species, F. x bohemica was able to produce high aboveground biomass in a nitrate-rich environment, and thus, although its assimilation rate was not significantly different than that of the other species, it was able to accumulate the highest amount of labeled nitrate in NO 3 − rich soil. These results in natural soils are consistent with recent findings [32], who showed that Japanese knotweeds did not have a better nitrogen uptake efficiency than that of native species but had a better nitrogen use efficiency (i.e., the biomass produced due to the uptake of nitrogen). On the other hand, U. dioica showed the highest root biomass, especially in low nitrate soil, and the highest nitrate assimilation rate. This resulted in a better nitrate accumulation in poor environments. From these findings, we could expect that F. x bohemica affected the other species by better competing for resources, especially in NO 3 − rich soil; in contrast, U. dioica was more efficient in poor soil. The results indeed showed that the impact of the presence of F. x bohemica on the native species performance in the mixture depended on the species. S. ebulus was either negatively affected or not affected. In U. dioica, we observed more aboveground biomass and a greater number of leaves in mixtures with F. x bohemica than that in monocultures. H. lupulus was globally negatively affected, except for an increase in root biomass only in high nitrate conditions. To our knowledge, this study is the first observation of facilitation induced by F. x bohemica in vascular species. However, these results should be considered with caution, given the number of published studies (experiments and in situ) that have found only negative effects of Fallopia spp. on the biomass of competing species (see [20]). We did not find a satisfactory explanation of the pattern on root biomass in H. lupulus. However, the literature already noted that U. dioica (1) could outcompete F. japonica (even if the relative yield of U. dioica was decreased by competition with F. japonica) [33], (2) is a strong competitor [34] and (3) is one of the few species able to coexist in Fallopia spp. stands [35,36]. Moreover, U. dioica was shown to be able to store nitrogen in its belowground organs (i.e., rhizomes and roots) as nitrate and/or amino acids (i.e., asparagin and arginin) forms [37]. U. dioica also exhibited a high nitrate reductase activity, as well as Fallopia japonica [38]. All these results may indicate that F. x bohemica and U. dioica have competing nitrogen acquisition strategies, but that their niches are sufficiently different to allow coexistence, depending on environmental conditions. This hypothesis would require more investigations to be confirmed.

Impacts of knotweed on Soil N-Functions and Microbial Communities: The Surprising Case of NO 3 − Poor Soil
Unexpectedly, only a few significant impacts of F. x bohemica in mixtures with native species were found on nitrifying and denitrifying communities and soil functioning under low nitrate conditions. Generally, soils with low mineral-N content are expected to be strongly influenced by plant-microorganism competition. Cantarel et al. [28] showed that an increase in mineral-N soil content (by fertilization) can partially release the effects of plants on microbial N activities and supported the hypothesis that plants impact microbial N communities via competition for N. Moreover, F. x bohemica is known to be (1) a very good competitive species [18,39] and (2) tolerant to various environmental conditions, from low to high nutrient soils [40]. We supposed that this tolerance to low nitrate conditions came from its ability to drive the nitrogen cycle resulting in more available nitrates and, more particularly, its ability to inhibit denitrification (BDI [11]). Moreover, in its native area, Fallopia spp. has been shown to influence the soil N cycle in increasing the concentration of soil N, during the local growth of the patch, accumulating in organic and mineral forms [41,42]. The conservation of nitrogen appears to be an adaptation to its native ecosystem. All of these results suggest that the most significant results on the N cycle would occur in low nitrate conditions. However, contrary to our expectations, only one significant difference occurred in low nitrate conditions: nitrification and copy number of ammonia-oxidizing archaea (AOA) increased in the F-S mixture compared to the S. ebulus and F. x bohemica monocultures. AOA is predominant in our soils, as often found in the rhizosphere [43], because of the higher AOA affinity for oxygen or carbon dioxide [44]. We also found that denitrifying gene abundances (i.e., nirK and nirS) were more abundant in the F-S mixture than in the F and S monocultures, without a change in denitrification. This result could reflect differences in microbial nitrite reductase community efficiency or phenotypic changes with faster or higher microbial turnover for denitrifying communities but lower efficiency for their nitrate/nitrite metabolism [45]. Overall, these results showed that competition between plants and nitrifying communities seemed to be released when S. ebulus was present in the mixture, suggesting an increase in nitrate produced by nitrification. However, denitrification was not affected, suggesting that some of the produced nitrates were used by plants. S. ebulus was strongly impacted by the presence of F. x bohemica in the mixture, with a large decrease in nitrate absorption and accumulation (more than 50%) in favor of Japanese knotweed assimilation and accumulation ( Figure S2). One of the key results in NO 3 − poor soil was the finding of no effect of F. x bohemica and/or mixtures on denitrification. Japanese knotweed is known to have the capacity to inhibit soil denitrification [10,11] by the production of plant metabolites (i.e., procyanidins; [46,47]). However, under low nitrate conditions, we did not find any change in denitrification rates between F. x bohemica and native species. Procyanidins are known to be complex polymeric polyphenols [48] that are potentially costly for F. x bohemica in limited nitrate conditions and are thus produced less in stressful conditions or are dependent on the developmental stage of the plants. Another explanation could be that produced procyanidins were ineffective in our experimental conditions because, for instance, they were sorbed into the soil or degraded.

Impacts of Knotweed on Soil N-Functions and Microbial Communities: NO 3 − Rich Soil, Ideal Playground for Knotweed?
In NO 3 − rich soil, nitrification was generally found to be reduced under mixtures and monocultures of natives, except for S. ebulus. This decrease was followed by a decrease in the nitrifying bacterial abundance of H. lupulus and in the nitrifying bacterial and archaeal abundances of U. dioica. All these modifications of the nitrification process and nitrifier communities could be seen as a mechanism implied in the niche construction of the native species following the arrival of an invasive species, which induced changes in the niche construction of native plant species and of the nitrifying community. In our experiment, denitrification was always reduced in mixture conditions compared to native monoculture conditions, suggesting that F. x bohemica can set up the biological inhibition of denitrification (BDI), regardless of the native plant associated with F. x bohemica. BDI is known to increase soil NO 3 − availability in soil [46,49]. Here, we found that only the growth and N accumulation of knotweed were impacted by the BDI in mixtures, with a higher biomass of F. x bohemica than that of native species. These results highlighted that the effect of procyanidins was beneficial for plant procyanidin producers and that F. x bohemica took advantage of soil nitrate not used by microorganisms following the BDI. The decrease in DEA was never coupled to the abundances of denitrifiers (nirK and nirS). These results suggest that the denitrification modifications induced by Japanese knotweed are weakly mediated by changes in the abundance of the denitrifying community, as already found by [10]. We provide evidence that changes in microbial activity can occur independently of changes in associated functional gene abundances, as found for nitrifying and denitrifying microbial composition [50]. Because of the relative stability of the microbial N community structure, investigations in RNA-based analyses would be required for future studies in invaded/non-invaded regions, as microbial functioning indicators of plant-microbe competition for N-substrates.

Conclusions
Our results highlight that in NO 3 − rich soil conditions, the invasive species F. x bohemica showed higher aboveground biomass and N accumulation than the three native species that were studied. In the presence of F. x bohemica (in its monoculture or in mixtures), denitrification was reduced compared to the soils with native monocultures, suggesting that F. x bohemica captures a part of the soil nitrate for its own growth to the detriment of the denitrifying microbial community and native plant species growth. In contrast, very few impacts on soil N functioning were found in NO 3 − poor soil, suggesting that its control on the activity and abundance of the microorganisms involved in denitrification may be less pronounced than previously expected or hidden by other processes.
Supplementary Materials: The following are available online at http://www.mdpi.com/1424-2818/12/4/156/s1, Table S1: Localization of plant populations, Table S2: Results on the effect of soil type on plant traits, Figure S1: Effects of native species on traits of F. x bohemica expressed as a percent change in the mixture relative to the invasive monoculture, Figure S2: Effects of mixtures on N uptake of F. x bohemica expressed as a percent change in the mixture relative to the invasive species in monocultures. Funding: This research was funded by FEDER-Plan Loire Grandeur Nature, Région Auvergne-Rhône-Alpes, Département de la Loire, Agence de l'eau Loire-Bretagne, and the FR BioEnvis -Université Lyon 1.