Fallopia japonica and Impatiens glandulifera are colonized by species-poor root-associated fungal communities but have minor impacts on soil properties in riparian habitats

Fallopia japonica and Impatiens glandulifera are major plant invaders on a global scale that often become dominant in riparian areas. However, little is known about how these species affect interactions in soil–plant systems. The aim of this study was to investigate the impact of both species on abiotic and biotic soil properties, with a special focus on fungi. We investigated eight sites along small streams invaded by F. japonica and I. glandulifera, respectively, and compared each with nearby sites dominated by the native species Urtica dioica. Three different types of samples were collected: bulk soil, rhizosphere soil and roots from invasive and native stands at each site. Bulk soil samples were analysed for soil physicochemical, microbial properties (soil microbial respiration and ergosterol) and soil arthropod abundance (Acari and Collembola). Soil respiration was also evaluated in rhizosphere samples. The fungal community composition of both bulk soil and roots were analysed using a metabarcoding approach. Soil physicochemical properties as well as soil microbial activity, fungal biomass and soil fungal operational unit taxonomic unit (OTU) richness did not differ between invaded and native riparian habitats, indicating only minor belowground impacts of the two invasive plant species. Soil microbial activity, fungal biomass and soil fungal OTU richness were rather related to the soil physicochemical properties. In contrast, Acari abundance decreased by 68% in the presence of F. japonica, while Collembola abundance increased by 11% in I. glandulifera sites. Moreover, root-associated fungal communities differed between the invasive and native plants. In F. japonica roots, fungal OTU richness of all investigated ecological groups (mycorrhiza, endophytes, parasites, saprobes) were lower compared to U. dioica. However, in I. glandulifera roots only the OTU richness of mycorrhiza and saprobic fungi was lower. Overall, our findings show that F. japonica and I. glandulifera can influence the abundance of soil arthropods and are characterized by lower OTU richness of root-associated fungi.

microbial respiration and ergosterol) and soil arthropod abundance (Acari and Collembola). Soil respiration was also evaluated in rhizosphere samples. The fungal community composition of both bulk soil and roots were analysed using a metabarcoding approach. Soil physicochemical properties as well as soil microbial activity, fungal biomass and soil fungal operational unit taxonomic unit (OTU) richness did not differ between invaded and native riparian habitats, indicating only minor belowground impacts of the two invasive plant species. Soil microbial activity, fungal biomass and soil fungal OTU richness were rather related to the soil physicochemical properties. In contrast, Acari abundance decreased by 68% in the presence of F. japonica, while Collembola abundance increased by 11% in I. glandulifera sites. Moreover, root-associated fungal communities differed between the invasive and native plants. In F. japonica roots, fungal OTU richness of all investigated ecological groups (mycorrhiza, endophytes, parasites, saprobes) were lower compared to U. dioica. However, in I. glandulifera roots only the OTU richness of mycorrhiza and saprobic fungi was lower. Overall, our findings show that F. japonica and I. glandulifera can influence the abundance of soil arthropods and are characterized by lower OTU richness of root-associated fungi.

Introduction
Invasive alien plants are widely accepted as a major cause of biodiversity loss in native communities, which, in turn, can negatively affect the entire ecosystem (Hooper et al. 2012). Besides effects on aboveground compartments, invasive plants also impact belowground compartments Abgrall et al. 2019) altering resources and ecosystem functions and processes (Crooks 2002;Ehrenfeld 2003;Waller et al. 2020).
Invasive plants can impact soil physicochemical properties, e.g. by increasing pH and water content (Leucanthemum vulgare, Ahmad et al. 2020) or lowering pH and increasing water content (Carpobrotus edulis, Novoa et al. 2014). Invasive plants also have various effects on the soil nutrient status. For example, invasion by Spartina alterniflora can increase nitrogen and organic carbon levels  or increase phosphorus storage while having no effect on C and N (Wang et al. 2019). Soil functions can be affected in terms of decreased litter decomposition and increased microbial activity Vilà et al. 2011) or reduced soil biological activity (Pehle and Schirmel 2015). Soil microbial communities, including fungi, are influenced by plant invasion with respect to biomass (e.g. increasing biomass with woody plant invasion; Liao and Boutton 2008) and shifts in the community composition and diversity (Ravit et al. 2003;Wolfe and Klironomos 2005;Si et al. 2013;Sapsford et al. 2022). Plant invasions can further influence the soil fauna (Belnap et al. 2005;Biederman and Boutton 2009;Abgrall et al. 2019), which support important soil functions such as decomposition (de Deyn et al. 2003;Tresch et al. 2019). For example, Collembola and Acari, often the dominant arthropod groups in soil (Petersen and Luxton 1982;Deharveng 1996;Behan-Pelletier 2003) respond to changes in the vegetation, including plant invasion (Salamon et al. 2004;Wissuwa et al. 2012). Altering the vegetation composition by replacing native plant species, could thus negativly affect soil arthropod abundance and diversity (Salamon et al. 2004;Wissuwa et al. 2012;Rusterholz et al. 2014). However, impacts of invasive plants on soil fauna are highly variable and depend on the habitat type (Abgrall et al. 2019), since many soil arthropod species are also dependent on abiotic soil properties (Cassagne et al. 2003;Bedano 2004).
Fallopia japonica (Houtt.) Ronse Decr., Polygonaceae and Impatiens glandulifera Royle, Balsaminaceae are highly abundant invasive plant species in many areas worldwide (Beerling and Perrins 1993;Beerling et al. 1994). Although their effects on soil physicochemical and biological properties are well studied, their impacts on fungal communities under field conditions are not well understood, especially due to the lack of studies using modern molecular analyses (metabarcoding). Soil fungi have important roles in ecosystem functioning (Öpik et al. 2006). Symbionts such as mycorrhizal species interact with plants through their roots by providing water and nutrients, and are essential for survival and health of most plant species (Smith and Read 2008;Kivlin et al. 2011;Tedersoo et al. 2014). In contrast to mycorrhizal fungi, root endophytes only interact with plant roots during a part of their life cycle. Mycorrhizal and endophytic fungi are especially beneficial under stressful environmental conditions and can enhance the performance, stress tolerance and disease resistance of plants (Smith and Read 2008;Rodriguez et al. 2009;Kivlin et al. 2011;Redman et al. 2011;Riess et al. 2014;Tedersoo et al. 2014). The community of mycorrhizal, endophytic as well as parasitic fungi, which also occur on living tissue, depends on the plant species and may thus be modified by invasion (Knogge 1996;Christian et al. 2016). In contrast to mycorrhizal and endophytic fungal groups, saprobic fungi are involved in the decomposition of particulate organic matter (Duguay and Klironomos 2000) and, hence, may be less affected by plant invasion.
The Asian knotweed, F. japonica, is a perennial, rhizome-forming species and one of the most invasive plants in Europe due to its ability to quickly produce high biomass (Beerling et al. 1994). It can affect soil characteristics (e.g., increasing mineral modifying soil nutrient stocks; Dassonville et al. 2007;Stefanowicz et al. 2017), C and N dynamics in soil (Koutika et al. 2007) or increase soil pH compared to adjacent non-invaded sites (Čerevková et al. 2019). Fallopia japonica shows allelopathic effects (Murrell et al. 2011;Parepa and Bossdorf 2016) with significant impacts on the soil food web structure (Abgrall et al. 2018), decreasing abundance and diversity of soil fungi, such as mycorrhizal fungi . Another highly invasive plant in Europe is the Himalayan balsam, I. glandulifera (Beerling and Perrins 1993). This annual herb has strong dispersal capabilities with masses of floatable seeds. Impatiens glandulifera is known for its impact on native communities with positive effects on soil fauna abundance (Tanner et al. 2013;Rusterholz et al. 2014), modifying soil fungal community composition by increasing fungal diversity and lowering arbuscular mycorrhizal fungal root colonization (Pattison et al. 2016;Gaggini et al. , 2019. In this study, we investigated riparian sites invaded by both species (F. japonica or I. glandulifera) and corresponding control sites dominated by the native common nettle, Urtica dioica L., Urticaceae. We aimed to investigate the chemistry of soil and effects on soil bacterial activity, soil arthropod abundance as well as soil fungal biomass, diversity and composition in both soil and root samples. For the fungi, we used a metabarcoding approach where sequences were assigned to operational taxonomic units (OTU). We assessed the following research questions: (1) Are soil physicochemical parameters (organic carbon content, pH value, water content) different between riparian sites invaded by F. japonica or I. glandulifera compared to native sites with U. dioica? (2) How does plant invasion by both species and soil physicochemical parameters affect soil microbial activity and soil arthropod abundance (Acari and Collembola)?
(3) How is the biomass, OTU richness, and composition of fungi in the bulk soil affected by F. japonica and I. glandulifera invasion and related to soil physicochemical parameters? (4) How does the richness and composition of root-associated fungal OTUs differ in the presence of invasive and native plant species?

Study area and target plant species
Study sites were located along small streams in forests around Landau (study area ~ 1000 km 2 ), Rhineland-Palatinate, Germany (Online Resource 1). Sixteen riparian sites containing either sites dominated by Fallopia japonica and Urtica dioica (N = 8 site pairs) or Impatiens glandulifera and U. dioica (N = 8 site pairs) were selected (Online Resource 2). The perennial herb Urtica dioica was selected as a reference species because it is a highly abundant, naturally occurring species at riparian floodplains in Central Europe (Edwards et al. 1998). Stands between the invasive and native plant species of each site pair were in direct vicinity (maximum distance up to 10 m), to minimize differences in soil conditions. All sites were at least 10 m 2 and characterized by homogenous vegetation and a cover of > 80% by their respective target species. Other common plants within the sites were Glechoma hederacea and Rubus subg. Rubus. In the sites dominated by F. japonica, an average of 3.4 plant species were identified, while in corresponding U. dioica sites, the average was 3.9. In sites dominated by I. glandulifera, the average number of species was 4.8 and in corresponding U. dioica sites, the average was 3.3. The forests around the target sites can be assigned to the phytocoenology classes Alnetea glutinosae and Quercetea robori-petraea according to Schubert. (2010), with Alnus glutinosa as the dominant tree species. Details of the vegetation at each site can be found in Online Resource 1.

Sampling
In each site, plant-specific samples were collected, namely bulk soil, rhizosphere soil and plant roots. Sampling was performed during September to October 2019. Bulk soil samples were randomly collected at the topsoil (0 to − 10 cm) using a bulb planter from three plant-specific stands and subsequently pooled and homogenized to one composite sample. Roots were taken from four randomly chosen individuals of each species using a shovel and homogenized to a composite sample by species and site. The soil attached to the roots after shaking the collected plant was separated from the plants by hand and collected as rhizosphere soil, in total 4 replicates, which were subsequently homogenized to a composite sample per species and site. Roots were further rinsed first with tap water, and finally with autoclaved water. A part of each soil type was frozen at − 20 °C and stored for analysis of soil microbial activity. The remaining bulk soil samples, rhizosphere soil samples and individual root samples were subsequently dried at 50 °C for 48 h.

Soil physicochemical properties
The bulk soil samples were used for analyses of soil physicochemical properties. Soil pH was measured using a pH meter (pH 3110 SET 2, Xylem Analytics, Weilheim, Germany) in a 0.01 molar CaCl 2 solution (DIN EN 15933:2012-11) with field moist soil. Water content in bulk and rhizosphere soil was determined gravimetrically according to ISO 11465:1993. Samples were further analysed for organic carbon (DIN ISO 10694:1996-08) by dry burning and total nitrogen (VDLUFA I,A2.2.5:2012) via thermal conductometry. These analyses were performed by Speyer Agricultural Research Institute (Speyer, Germany).

Soil microbial activity
Soil microbial activity was investigated separately for the bulk soil and the rhizosphere soil. Basal and glucose-induced soil microbial activity were assessed by the MicroResp™ method (Campbell et al. 2003). The evolved CO 2 after incubation is indicative for soil respiration, which is an expression of soil microbial activity (Creamer et al. 2016a, b). The frozen samples were thawed at 6 °C and sieved to 2 mm. Then, they were adjusted to 45% of water holding capacity, following the recommendations of the producer. The substrates (water or glucose) were added to the soil samples followed by an incubation of 6 h at 21 °C. CO 2 was measured by a shift of the cresol red colour in the deep well plate system at 572 nm (Nanophotometer Infinite® M200, Tecan, Switzerland). Quantitative analysis was done using known concentrations of sodium bicarbonate and an excess of hydrochloric acid (Campbell et al. 2003) by interpolation of data in the calibration curve.

Soil arthropod abundance
We focussed on the highly abundant and functionally important soil arthropod taxa Acari and Collembola. We combined two methods, pitfall traps and Berlese funnels, for sampling. Four pitfall traps (plastic cups with an opening diameter of 6.5 cm) were placed one meter apart in each sampling site. Traps were set between 24th August and 15th September 2019 and were emptied and refilled every ten days. A 1:4 propylene glycol:water solution with some drops of detergent was used as a trapping liquid. Captured individuals were stored in 70% ethanol. Soil organisms were additionally extracted from the soil bulk samples using Berlese funnels. 200 ml of fresh soil per site from the collected bulk soil was filled into a Berlese funnel for one week (4-11th September). A constant light source was attached to the top of the funnel to induce downward movement for the soil organisms. Organisms were collected at the bottom of the funnel in vials containing 20 ml of 70% ethanol and the collected invertebrates were identified to class level. The abundance of Acari and Collembola individuals from pitfall traps and Berlese funnels were combined for the analyses.

Soil fungal biomass
Ergosterol was used as a proxy of soil fungal biomass in bulk soil samples, based on the method of Gong et al. (2001). 4 g of air-dried and milled soil (Planetary Micro Mill PULVERISETTE 7 premium line, Fritsch, Idar-Oberstein, Germany) was extracted with 12 ml of methanol for 60 min on a horizontal shaker (Kreisschüttler 3015, GFL, Burgwedel, Germany). The mixture was sonicated for 10 min (DT 514H, Bandelin electronics, Berlin, Germany) and centrifuged for 10 min at 2000 g (Universal 320, Hettich Lab Technology, Tuttlingen, Germany). The supernatant was ultra-centrifuged for 3 min at 7270 g (Micro centaur, MSE, London, UK) prior to highperformance liquid chromatography (HPLC) analysis. An aliquot of 20 μL of the supernatant was analysed via HPLC with UV detection at 282 nm (HPLC 1200 series, Agilent technologies, Santa Clara, CA, USA), equipped with a C18 LiChrospher® column (LiChrospher RP-18e, 5 μm, 100 Å, 250 × 4.6 mm, Merck, Darmstadt, Germany). Methanol (flow rate 1.7 ml min −1 ) was used as a mobile phase under isochratic condition. The temperature of the column was set at 38 °C. The limit of detection (LOD) of the method was 0.06 mg kg −1 (Meyer et al. 2021).

Soil and root fungal OTU richness and ecological groups
Fungal operational taxonomic unit (OTU) richness was analysed separately for bulk soil and root samples of the investigated plant species. Next-generation sequencing was carried out by the company Advanced Identification Methods (AIM) (Munich, Germany). For detailed processing of sequences, see Online Resource 3. In short, the internal transcribed spacer 2 region of the nuc-rDNA (ITS2) was amplified using the primers ITS3tagmix/ITS4ngs (Tedersoo et al. 2014(Tedersoo et al. , 2015. Quality filtering reduced the initial 3,287,492 full sequences to a final dataset with 241,258 unique, non-singleton sequences. Clustering with a threshold of 2% and de novo chimera detection resulted in 6,468 OTUs (Online Resource 4). OTUs were matched against NCBI GenBank database (ncbi. nlm.nih.gov; October 2019; see Online Resource 5). Only 5,012 OTUs matched with sequences assigned to Fungi with more than 95% identity and were used for further analysis. In general, there is a large discrepancy between the described and the estimated number of species (Hawksworth and Rossman 1997). Hence, some OTUs cannot be assigned to any taxon or only at a higher taxonomic level. To increase the quality of database labelling, ambiguously attributed OTUs (e.g., "uncultured eukaryote" or "fungal endophyte") were checked manually and assigned to phyla using MycoBank (mycobank.org; Robert et al. 2013) as a reference database.
Prior to the statistical analysis, four main fungal ecological groups were identified: (1) Mycorrhiza is a widespread well-known symbiotic formation supporting both fungal and plant partner (Smith and Read 2008). In our study area, two main mycorrhizal groups are common: the arbuscular mycorrhizal fungi (AMF) comprising only Glomeromycota (Parniske 2008) and the ectomycorrhizal fungi (ECM) including mainly Basidiomycota, e.g. the large genus Cortinarius (Garnica et al. 2016), and Ascomycota, e.g. the species-rich order Helotiales (Tedersoo et al. 2009). In Europe, AMF typically occurs in herbaceous plants whereas ECM is usually found in woody species.
(2) Root endophytes include phylogenetically diverse fungi and occurs ubiquitously in terrestrial plants across a variety of ecosystems. At some stage of their complex life cycle, they inhabit plant tissues without causing any obvious disease symptoms (Rodriguez et al. 2009). Species-rich endophytic fungal taxa are for example Trichoderma (Harman et al. 2004) and many members of the Sebacinales (Oberwinkler et al. 2013). (3) Saprobic fungi break down energy-rich organic substances from tissue of dead organisms or soil, where they play a key role beside bacteria (Lebreton et al. 2021 Rodrigues and Fonseca 2003) and Ustilago (Bauer et al. 1997). Annotation of fungal OTUs to the ecological groups was realized in a two-step process: First, a literature search for taxa with well-known ecology was carried out. This is reasonable on genus or higher taxonomical range (Põlme et al. 2020). In a subsequent second step, a screening for OTU ecology was carried out using information deposited in metadata in the NCBI database (Online Resource 6). This approach allowed us to assess the fungal OTU richness of soil and root-associated fungi as well as the OTU richness in the respective ecological group (Zanne et al. 2020).

Data analysis
The data was analysed using R version 4.0.5 (R Core Team 2021). The package glmmTMB version 1.0.2.1 was used to create generalized linear mixed models (Brooks et al. 2017). The models were run separately for the site pairs of F. japonica / U. dioica and I. glandulifera / U. dioica, respectively. For the comparison of soil properties (organic carbon, total nitrogen, pH value, water content) between invaded and native sites, only plant species (factor with the two levels "invasive" [either F. japonica or I. glandulifera] and "native" [U. dioica]) was included as an explanatory variable and a Gaussian distribution was used. The study site was included as a random effect to account for our paired study design. The models for questions 2-4 incorporated plant species (factor with two levels), pH value (continuous), organic carbon content (continuous), and water content (continuous) as explanatory variables, and study site as a random effect. We decided on these variables by investigating the intercorrelation of the different explanatory variables (Online Resource 7). For example, organic carbon and total nitrogen content (r = 0.89) were highly correlated. This allowed us to simplify the model by excluding nitrogen content. The response variables used were the soil microbial activity in bulk and rhizosphere soil, soil organism abundance, fungal biomass, total fungal OTU richness, and the OTU richness of the different fungal ecological groups (analysed separately for soil fungi and plant-associated fungi). Depending on the analysed variables, different family functions for the models were used. For the fungal OTU richness and the abundance of soil organisms, a negative binomial model for overdispersed count data was applied (Lindén and Mäntyniemi 2011). For the fungal biomass and both soil microbial activities, a gaussian model was used. The model results are presented by showing the estimate (Est.), z-value and P-value.
To reveal differences in the fungal OTU composition between invaded (either F. japonica or I. glandulifera) and native (U. dioica) sites and between fungi associated with invasive and native plants, permutational multivariate analyses of variances (PER-MANOVA) were performed (command adonis2 in the package 'vegan', Oksanen et al. 2019). The parameters organic carbon content, pH value and water content were used as further environmental variables to correspond to the previous models. Because of the paired study design, the study sites were used as strata in the formula. The PERMANOVA were based on presence-absence data for fungal OTUs occurring in more than one site or plant and the Jaccard index was used as a similarity measure. PER-MANOVA were performed for soil fungal OTU composition (question 3) and root associated fungal OTU composition (question 4). Variation of fungal OTU compositions were visualized using nonmetric multidimensional scaling (NMDS) with the command 'metaMDS' in R package 'vegan'.

Soil physicochemical properties between invaded and native riparian vegetation
Soil properties (organic carbon content, nitrogen content, pH value and water content) did not significantly differ between sites invaded by Fallopia japonica or Impatiens glandulifera and the paired sites dominated by native Urtica dioica (Table S1).

Relations of soil microbial activity to riparian plant invasion and soil physicochemical properties
Water content and soil organic carbon were significantly related to bulk soil microbial activity in F. japonica site pairs. The basal soil respiration decreased with increasing water content (Est. = − 0.04, z = − 2.91, P = 0.004) and increased with increasing soil organic carbon (Est. 0.59, z = 2.86, P = 0.004). However, the glucose-induced respiration was not affected by soil physicochemical parameters (Table S2). We did not find significant effects in basal and glucose-induced soil respiration for bulk soils from the I. glandulifera site pairs. Moreover, no significant effects were found for rhizosphere soils from both F. japonica or I. glandulifera site pairs (Table S2).

Soil fungi in relation to plant invasion and soil physicochemical properties
Analysis of bulk soil samples from the site pairs of F. japonica / U. dioica revealed no differences in total fungal biomass measured as ergosterol. Fungal biomass was positively correlated with organic carbon content (Est. = 0.88, z = 2.56, P = 0.010) while pH and water content had no influence (Table S3).
Regarding the site pairs of I. glandulifera and U. dioica, neither fungal biomass nor fungal OTU richness were different between the invaded and uninvaded sites (Table S3). However, fungal biomass increased with increasing organic carbon content (Est. = 1.59, z = 6.15, P < 0.001) but decreased with higher soil water content (Est. = − 0.11, z = − 2.52, P = 0.012). In terms of the fungal OTU richness, only soil water content was found to have effects (Table S3). Specifically, the OTU richness of endophytes (Est. = − 0.02, z = − 2.06, P = 0.039), saprobes (Est. = − 0.02, z = − 2.66, P = 0.008) and parasites (Est. = − 0.02, z = − 2.17, P = 0.030) decreased with higher water content. The soil fungal community did not differ between I. glandulifera and U. dioica sites and none of the soil parameters had a significant effect (Table 1, Fig. 2b). For details of the most common fungi taxa found in soil samples of both F. japonica and I. glandulifera site pairs, see Online Resource 8.

Discussion
In this study we analysed the impact of riparian plant invasions on soil, arthropods, and fungi by comparing sites invaded by Fallopia japonica or Impatiens glandulifera to corresponding sites inhabited by the native Urtica dioica. We found that both invasive species had only minor effects on soil physicochemical properties, soil microbial activity, soil arthropod abundance and soil fungal biomass, OTU richness, and composition. We further investigated root-associated fungal communities of the three plant species. As a key finding, we observed strong reductions in

Riparian plant invasion has low impact on soil physicochemical properties and soil microbial activity
None of the investigated soil physicochemical parameters (organic carbon content, pH, water content) nor soil microbial activity were significantly affected by the presence of an invasive plant invader investigated here. While soil organic carbon content seems to be less influenced (Aguilera et al. 2010), previous studies have found that micronutrient levels can increase in presence of F. japonica, especially under poor nutrient conditions ). We studied riparian and relatively nutrient-rich habitats along streams, which could be one reason why we did not detect an impact of F. japonica on soil physicochemical parameters. Moreover, impacts of plant invasions often change over time (Strayer et al. 2006;Dostál et al. 2013), but information about the age of our studied plant stands was not available. For I. glandulifera, our results conform with findings of Čuda et al. (2017) who also found only minor impact of this species on soil physiochemical properties. Since root traits such as root exudation can influence soil functions (Bardgett 2018), the lack of an effect on soil microbial activity (measured by soil respiration) was unexpected. Both F. japonica and I. glandulifera have also been shown to negatively impact soil microbial activity ). However, Fallopia spp. may mainly affect anaerobic respiration through root exudates (for example procyanidins), while aerobic respiration is less impacted (Bardon et al. 2016(Bardon et al. , 2017. This observation might be directly linked to the lack of effects in the fungal community and biomass among sites (see below), as fungi are considered a main driver of microbial activity in soils (Žifčáková et al. 2016;Xu et al. 2019).
Fallopia japonica and I. glandulifera affect soil arthropod abundance The analysed soil arthropods Acari and Collembola were strongly affected by the presence of the invasive plants, though in contrasting ways. The abundance of Collembola significantly increased in I. glandulifera stands but was not affected by F. japonica. While the impact of I. glandulifera was previously linked to its annual life cycle, which leads to a higher availability of food for collembolan herbivores (Tanner et al. 2013), the lack of an effect by F. japonica was unexpected due to its known allelopathic impacts (Kato-Noguchi 2021) on soil microfauna (Abgrall et al. 2018). However, detritivores are often less affected by invasive plants than other trophic groups, which might be explained by often increased plant Relations between soil fungal OTU composition and environmental variables were tested with permutational multivariate analyses of variances, see Table 1 productivity in invaded habitats Ehrenfeld 2010;Schirmel et al. 2016). In our study, invasion by F. japonica negatively impacted Acari abundance while I. glandulifera had no significant effect. For I. glandulifera, this is in accordance with earlier reports, where its presence had differential effects on different groups of Acari, which lead to no effect in total abundances . The lower Acari abundances in F. japonica sites may be explained by its phenolic compounds in the leaf litter, which have negative impacts on saprobic species (Skubala and Mierny 2009;Abgrall et al. 2018). Moreover, reduced Acari abundances under F. japonica might be a result of changes in the soil food web (Abgrall et al. 2018), although this was not demonstrated in our study.
In the case of F. japonica, we uncovered a positive and negative relationship between Collembola abundances and organic carbon contents and pH, respectively. For I. glandulifera sites, the opposite patterns for both variables were observed. This supports previous studies that have documented the importance of soil organic matter for shaping soil fauna communities (Rendoš et al. 2016). The influence of pH might also be indirect, by altering the soil microbial community (Fierer et al. 2009). Acari abundances were also affected by soil pH. The abundances increased with pH value in F. japonica site pairs but decreased in I. glandulifera site pairs. The effect of pH value on Acari abundance may be explained by the fact that pH differences are associated with different microhabitats, which influence Acari abundance and communities (Nielsen et al. 2010;Wissuwa et al. 2013).
Soil physicochemical properties are more influential for soil fungi than the presence of invasive plants The presence of invasive plant species had little effect on the soil fungal biomass, OTU richness and composition in the context of this study. We only found an increase in saprobic fungi richness in soils from F. japonica sites compared to those from U. dioica. Saprobic fungi may be promoted by increased amounts of litter introduced by F. japonica which is degraded by saprobes (Lebreton et al. 2021). This finding is unexpected and contrasts with earlier studies, where plant invasion induced changes in fungal diversity and function (Wolfe and Klironomos 2005;Si et al. 2013). One explanation might be that the plant invasions did not significantly alter physicochemical properties compared to the uninvaded sites, which are known to strongly influence soil fungi abundances and diversity (for example Lauber et al. 2008;Wakelin et al. 2016;Yang et al. 2017). This indicates that plant invasion may only affect the soil fungi community if it alters the physicochemical properties of the habitat. Consequently, the soil physicochemical properties were much more important in explaining soil fungi biomass and OTU richness. It is well known that water availability, organic carbon and pH have major effects on the fungal species composition (Li et al. 2015;Xiao et al. 2018). Here, the fungal biomass was affected by the soil organic carbon content providing more resources for saprobic fungal species supporting their development (Bossio and Scow 1998;Yao and Wilson 2000;Drenovsky et al. 2004). However, increased organic carbon also leads to a reduction in overall as well as mycorrhizal OTU numbers in the F. japonica site pairs in our study. Contrary to this observation, Zhu et al. (2020) showed that mycorrhizal species richness increases with organic carbon content. Other studies observed a decrease in fungal diversity along a gradient of the same variable (Liu et al. 2015). This might be because the interaction of mycorrhizal taxa and organic carbon in the soil is dependent on the local conditions (Hanson et al. 2008;Wei et al. 2019).
At F. japonica site pairs, pH value was associated with an increase of fungal OTU richness in all ecological groups except mycorrhiza (endophytes, parasites, and saprobes). In contrast, at I. glandulifera site pairs, pH had no significant influence on fungal OTU richness. Moreover, we did not observe an effect of pH on fungal community composition. Our findings are thus surprising, because pH value is known to affect especially mycorrhizal species richness and composition in various habitats (Porter et al. 1987;An et al. 2008Davison et al. 2015;Řezáčová et al. 2019). One possible explanation for this contradictory result might be that the differences in the pH values across our study sites were too small (on average 5.7-6.2; Table S1) to impact the fungal OTU richness and composition. In other studies, these pH ranges were much wider, for example between 4.8 and 7.5 (Porter et al. 1987) or 4.0 and 7.6 (Řezáčová et al. 2019). Increasing soil water content negatively affected mycorrhizal fungal OTU richness in the F. japonica site pairs and fungal biomass as well as fungal OTU richness of endophytes, parasites, and saprobes in the I. glandulifera site pairs. Fungi have a lower moisture optimum than other soil microbes and have a lower competitive ability under high soil moisture conditions (Kouyeas 1964). This has also been shown by Cavagnaro (2016), where higher soil moisture negatively impacted mycorrhizal growth and colonization of plants, and Deepika and Kothamasi (2015) who found lower levels of arbuscular mycorrhizal fungi (AMF) diversity in flooded soils.
Invasive plant species are characterised by low OTU richness of root fungal communities The root-associated fungal communities were composed of low phylogenetic OTU richness in the invasive species, and were remarkably distinct to the native counterpart U. dioica. Irrespective of the fungal ecological group, significantly fewer root-associated OTUs were found for F. japonica compared to U. dioica. The second species investigated, I. glandulifera, was also characterized by reduced rootassociated fungal OTUs, however significant reductions compared to U. dioica were observed only for mycorrhiza and saprobic taxa. The lower endophyte and mycorrhizal OTU richness in F. japonica and I. glandulifera roots were likely caused by changes in host-symbiont interactions. Symbionts are often not host specific but functional variability between different plant-mycorrhiza combinations have been observed (Klironomos 2000). Therefore, an invasive plant species may not be able to benefit from local soil fungi which might reduce mycorrhizal and root endophytic diversity (Čuda et al. 2020). This is in line with Callaway et al. (2011) who reported that AMF from a native range stimulated stronger positive feedbacks in invasive Robinia pseudoacacia than AMF from an invasive range. On the other hand, root endophytes of Centaurea stoebe were reported to directly increase the competitive effects of the introduced plant species (Aschehoug et al. 2012). Root endophytes of two Centaurea species (one invasive and one native) were also found to be different, indicating that endophytes may not be able to form associations with a variety of plants (Geisen et al. 2017). Garnica et al. (2022) reported a strong colonization rate for the root endophytic species Serendipita herbamans in F. japonica in a lab study, especially in low-nutrient conditions. Secondary metabolites inhibit the growth of especially mycorrhizal fungi (Pinzone et al. 2018), which might also explain the lower number of rootassociated fungi found in F. japonica compared to I. glandulifera. Moreover, I. glandulifera is able to associate to some degree with AMF (Tanner et al. 2014) and according to Pattison et al. (2016), I. glandulifera can affect the AMF and fungal endophyte diversity.
The lower number of mycorrhizal OTUs in roots of both F. japonica and I. glandulifera compared to native U. dioica indicates reduced symbiotic interactions. The root associated fungi in the invasive species include no AMF, which indicates that these species may not be able to associate well with symbiotic fungi in the invasive range. Some Glomeromycota were however detected in U. dioica roots (Rhizophagus spp., see Online Resource 8). Root endophytic OTUs were lower in F. japonica, which might be caused by the absence of the most common and diverse genus Trichoderma (Harman et al. 2004). One reason might be due to the secondary metabolites utilized by F. japonica. We also found reduced numbers of parasitic fungal OTUs in root samples of F. japonica compared to U. dioica. The most common genus in U. dioica, Fusarium, was absent in F. japonica. The absence of directly associated major fungal parasites inhibiting its growth in the invasive range might contribute to the invasion success of F. japonica. This is in accordance with the enemy release hypothesis, which states that invasive species have increased success because of lower numbers of natural enemies (Mitchell and Power 2003). Root associated OTUs classified as saprobic were lower in both invasive species compared to the native species. Saprobic fungi may enter tissue to break down dead matter (Lebreton et al. 2021). The reduced number of OTUs in the invasive species might indicate that a smaller subset of species is able to metabolize the tissue of the two invasive plant species. Mincheva et al. (2014) showed that F. japonica is a low-quality food source and had a different fungal composition compared to litter from native sites.
Like the soil fungal OTU abundances, root associated fungi were affected by some soil parameters as well. The most important factor was soil water content, having a predominantly negative effect on the overall fungal OTU richness as well as the ecological groups mycorrhiza and parasites in the I. glandulifera site pairs. In F. japonica site pairs, only the number of mycorrhizal OTUs was negatively associated with increasing water content, while higher water content was positively correlated with OTU richness in endophytes and parasites. The negative effect of higher moisture on mycorrhiza has previously been shown by Cavagnaro (2016) and Deepika and Kothamasi (2015). Tree root associated saprotrophic and parasitic community structure has also previously been shown to correlate with soil moisture, where higher soil moisture benefited only some species (Burke et al. 2009).

Conclusion
We observed only minor effects of the investigated plant invasive species Fallopia japonica and Impatiens glandulifera on soil physicochemical properties, soil microbial activity, soil fungal biomass and OTU richness. In contrast, the abundances of Acari and Collembola were impacted, indicating that both species can have direct impacts on soil biota. However, these impacts were not uniform between F. japonica and I. glandulifera, which implies that these effects are species-specific and may depend on the local soil biota composition. Looking at tissue-entering fungi, F. japonica clearly had lower fungi OTU richness while I. glandulifera was characterized only by fewer mycorrhizal and saprobic OTUs. This was further related to soil characteristics, especially water content and pH decreased fungal OTU richness. Overall, our findings indicate that F. japonica and I. glandulifera invasions may have various and partly contrasting impacts to previous studies on belowground ecosystem compartments such as soil invertebrates and fungi. Therefore, effects of these plant invaders seem to be highly context-dependent and site-specific. Furthermore, our results are in line with the enemy release hypothesis, indicating that these invasive species, especially F. japonica, may be less associated with fungal parasites in the invasive range, allowing them to perform better than native species.

Conflicts of interest
The authors declare that they have no conflict of interest.

Data availability
The authors declare that the data supporting the findings of this study are available within the article and its supplementary information files.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.