Skip to main content

Potential distribution and ecological impacts of Acmella radicans (Jacquin) R.K. Jansen (a new Yunnan invasive species record) in China

Abstract

Background Acmella radicans

(Jacquin) R.K. Jansen is a new invasive species record for Yunnan Province, China. Native to Central America, it has also been recently recorded invading other parts of Asia. To prevent this weed from becoming a serious issue, an assessment of its ecological impacts and potential distribution is needed. We predicted the potential distribution of A. radicans in China using the MaxEnt model and its ecological impacts on local plant communities and soil nutrients were explored.

Results

Simulated training using model parameters produced an area under curve value of 0.974, providing a high degree of confidence in model predictions. Environmental variables with the greatest predictive power were precipitation of wettest month, isothermality, topsoil TEB (total exchangeable bases), and precipitation seasonality, with a cumulative contribution of more than 72.70% and a cumulative permutation importance of more than 69.20%. The predicted potential suitable area of A. radicans in China is concentrated in the southern region. Projected areas of A. radicans ranked as high and moderately suitable comprised 5425 and 26,338 km2, accounting for 0.06 and 0.27% of the Chinese mainland area, respectively. Over the 5 years of monitoring, the population density of A. radicans increased while at the same time the population density and importance values of most other plant species declined markedly. Community species richness, diversity, and evenness values significantly declined. Soil organic matter, total N, total P, available N, and available P concentrations decreased significantly with increasing plant cover of A. radicans, whereas pH, total K and available K increased.

Conclusion

Our study was the first to show that A. radicans is predicted to expand its range in China and may profoundly affect plant communities, species diversity, and the soil environment. Early warning and monitoring of A. radicans must be pursued with greater vigilance in southern China to prevent its further spread.

Peer Review reports

Background

Invasive alien species are defined as non-native species which have become naturalized through establishing self-sustaining populations in semi-natural or natural ecosystems [1]. Invasive species generally are species that may spread and colonize relatively large geographic areas [2]. The spread of invasive alien plant species has resulted in biodiversity loss, environmental problems, and great economic harm [3, 4]. In a review of biotic and economic impacts, Rai et al. [5] summarized the overall impacts as altering the composition but also structure and function, leading to ecological, economic, and social impacts. Invasive plants can reduce species richness and abundance of native species, alter soil chemistry and soil microbial communities, and disrupt community structure and ecosystem processes [6,7,8]. The effects of invasive plants on the biota are often apparent via comparisons with native plant species, over which they frequently have competitive advantages and fewer natural enemies present in the recipient environment [9, 10]. Rapid growth rates by comparison to native species frequently lead to the invasive plants dominating a given invaded environment at the expense of native communities [11]. The invasion process often includes four sequential stages: transport, introduction, establishment, and spread [12]. Early detection of invasive species new to an area, including assessing their potential impact, is the most cost-effective management strategy, proactively managing a species that has not yet spread into large areas of suitable habitat.

China, as the world’s third largest country, includes five climate zones: old-temperate, temperate, warm-temperate, subtropical, and tropical, which results in a rich biodiversity [13]. Almost any potential invasive alien species originating in one of these diverse climate zones throughout the world may find a suitable niche in China. Thus, with its varied topography, climatic conditions, and ecosystems, China is highly vulnerable to invasive alien species. Invasion by non-native species has accelerated via recent rapid increases in economic development, international trade, global tourism, and climate change in China [14, 15]. Many of the problematic invaders have been members of the Asteraceae family, such as Mikania micrantha, Ageratina adenophora, and Eupatorium odoratum [16].

Another more recent invasive plant invader in the Asteraceae family is Acmella radicans (Jacquin) R.K. Jansen. The plant is native to Central America but has recently been discovered in parts of Asia, including several new occurrence records in India [17,18,19]. In China, A. radicans was first found to be naturalized in Anhui Province in 2014 [20]. In 2017, A. radicans was first discovered as a potential issue in Yunnan Province during a survey of invasive alien plants. It was recorded and collected from Banhong Township, Cangyuan County. This plant was identified as A. radicans after a thorough survey of literature and expert opinion and the specimen was deposited in the Agricultural Environment and Resource Research Institute, Yunnan Academy of Agricultural Sciences, Kunming, China. Further field investigations and monitoring throughout the province showed A. radicans was widely distributed in Baoshan City and Lincang City as a serious invasive species, primarily invading farmland, tea gardens, orchard land, roadsides, and ditches. Even so, A. radicans is relatively early in its spread in the area, i.e., early in the fourth invasion stage, and thus is a good target for an “early detection, rapid response” approach. To carry out this approach, we need to understand its potential impacts and potential for further spread as examined in the current study. To date, literature predicting its ability to spread, and its potential ecological impacts is scarce, with the literature in Asia mostly limited to records of its occurrence without reference to impacts [17,18,19].

Species Distribution Models (SDMs) are valuable tools for studying the potential suitable distribution of invasive species [21] such as A. radicans. Commonly used SDMs include Maximum Entropy Model (MaxEnt) [22], Genetic Algorithm for Rule-set Prediction Model (GARP) [23], Biological Climatic Model (BIOCLIM) [24], and Match Climate and Compare Location Model (CLIMEX) [25]. Among these models, MaxEnt is the most widely used SDM at present. MaxEnt determines predicted species distributions through calculation of parameters with maximum entropy as a product of interactions between species and their environment [22]. MaxEnt and spatial analysis techniques have been used to predict the potential distribution of invasive plants such as Xanthium italicum [26], Alternanthera philoxeroides [27], Erigeron canadensis [28], Lantana camara [29], and many others. Many studies have confirmed the importance of MaxEnt model optimization, constrained by species distribution pattern, model conditions, validity of model assumptions and applicability of optimization measures, but model optimization may not necessarily improve the prediction ability [30,31,32]. Although MaxEnt model performs well in discrimination, due to many variables affecting species distribution, statistical models based on available distribution data often cannot reflect the complexity of the ecological demands of a given species [21, 33, 34]. As in any model system then, MaxEnt models provide a reasonable approximation of predicted distribution but are limited in the absence of detailed ecological information. Since SDMs focus on species distribution patterns, the relationship between predicted results and species ecological characteristics, such as population density, number of breeding populations, reproductive success, and other factors warrant further attention [34, 35]. Thus it is useful in studies such as ours to evaluate the model results in terms of the biology and ecology of the study species.

The objectives of the present study were to explore the potential distribution of A. radicans using a species distribution model (MaxEnt), examine ecological impacts of this new invasive plant species in China, and provide useful information to help minimize and mitigate invasion by A. radicans.

Results

Model optimization and evaluation of prediction precision

Based on 396 A. radicans distribution records and 19 environmental predictor variables, the potential distribution in China of A. radicans was predicted using MaxEnt. With MaxEnt default parameter settings, RM = 1, FC = LQPTH, and delta AICc = 30.74. When RM = 1, FC = LPTH, delta AICc = 0, the model is optimal and omission rate at 5% is lower than the model under the default parameters (Table 1), which is 13.04% lower than the default value. Therefore, we set RM = 1 and FC = LPTH as the final parameters of the model. The AUC value of the simulated training under this parameter was 0.974 (Fig. 1), indicating a high level of confidence in model predictions.

Table 1 Evaluation results of MaxEnt model under different parameter settings
Fig. 1
figure 1

Receiver operating characteristic (ROC) test of the MaxEnt model (AUC is the area under curve)

Evaluation of environmental variables

The Jackknife test was utilized to assess the degree of influence of environmental variables on the predicted outcome and thus determine how different variables contributed to the predicted potential distribution of A. radicans. The results of the Jackknife test (Fig. 2) showed that the top 8 environmental variables ranked from highest to lowest were isothermality (Bio3), precipitation of wettest month (Bio13), topsoil silt fraction (T-SILT), temperature annual range (Bio7), topsoil TEB (total exchangeable bases) (T-TEB), annual mean temperature (Bio1), topsoil1 CEC (T-CEC-CLAY), and topsoil organic carbon (T-OC) (i.e., ranked by regularization training gain when using only individual variables). Contribution rate and permutation importance value are also important indicators to assess the degree of influence of each environmental variable (Table 2). Precipitation of wettest month (Bio13), isothermality (Bio3), topsoil TEB (T-TEB), precipitation seasonality (Bio15), topsoil CEC (T-CEC-CLAY), precipitation of warmest quarter (Bio18), topsoil silt fraction (T-SILT), and annual mean temperature (Bio1), constituted the eight environmental variables contributing most to the potential suitability distribution of A. radicans, as evaluated via contribution rate and permutation importance values. The cumulative contribution rate reached 89.60% and the cumulative permutation importance value reached 82.10%. Combining the Jackknife test results, contribution rate, and permutation importance value analysis, we concluded that precipitation of wettest month (Bio13), isothermality (Bio3), topsoil TEB (T-TEB), and precipitation seasonality (Bio15) were the dominant variables affecting the potential suitability distribution of A. radicans (Table 2).

Table 2 Contribution rate and permutation importance value of each environmental variable (%)
Fig. 2
figure 2

Jackknife test of the environmental variables

Potential distribution of Acmella radicans

The potential suitable geographic area for A. radicans in China is mainly concentrated in the southern region (Fig. 3). The area characterized as slightly suitable covers 282,332 km², accounting for 2.94% of Chinese mainland. The slightly suitable area is chiefly distributed in eastern Sichuan, western Chongqing, southern Guizhou, northern Guangxi, all of Fujian, and a small amount in Hunan, Jiangxi, Zhejiang, and Guangdong.

Fig. 3
figure 3

Potential suitable habitats of Acmella radicans predicted by the MaxEnt model. a Global potential geographical distribution (excluding Antarctica). b Potential geographical distribution in China

The moderately suitable area of A. radicans in China covers 26,338 km², comprising 0.27% of the Chinese mainland. It is principally distributed in the western, southwestern, and south-central areas of Yunnan. Other potential distribution locales include southeastern Tibet, a small area encompassing southern Sichuan and northern Yunnan, the junction of southwestern Guizhou and northern Guangxi, south-central Hunan, central Hainan, and the eastern coastal region of Zhejiang.

The area deemed highly suitable for A. radicans in China covers 5425 km², accounting for 0.06% of the Chinese mainland. These highly suitable areas for A. radicans were concentrated in Yunnan and scattered in the west, southwest and south-central areas of Yunnan.

Plant species and densities

The total area infested by A. radicans was over 5560 ha in Yunnan, primarily distributed in farmland, orchard land, roadsides, and wasteland in Baoshan City, and Lincang City. A total of 17 plant species belonging to 17 genera and 7 families were recorded within the study quadrats (Table 3). All 17 species were herbaceous and included 14 annual plants and 3 perennial plants, accounting for 82% and 18% of all species, respectively. Within communities of A. radicans, population densities of A. radicans were 19.81, 23.83, 40.41, 44.13, and 60.66 corresponding to study years 1, 2, 3, 4, and 5, respectively. Population density of most plant species declined markedly as the duration of the presence of A. radicans increased, however population density of A. radicans itself, increased significantly with increasing A. radicans growing period within five study years. Eight plant species A. radicans, Ageratum conyzoides, Bidens pilosa, Borreria latifolia, Chloris virgata, Cynodon dactylon, Digitaria sanguinalis, and Setaria plicata occurred over all five years and were dominant in A. radicans communities, having higher population densities than other species.

Table 3 Plant species and densities under different surveyed years from 2018 to 2022 (individual/m2) (mean ± SD)

Effects of Acmella radicans on plant importance values

There was variation in the response of individual species to increased duration of time over which A. radicans was present in the system. Importance values of A. conyzoides, B. pilosa, B. latifolia, C. virgata, Commelina benghalensis, C. dactylon, D. sanguinalis, Echinochloa hispidula, Eleusine indica, and S. plicata were higher within A. radicans communities, and declined markedly as the duration of the A. radicans infestation increased. Importance values of A. radicans, increased significantly the longer A. radicans was present. No trends were apparent for other species because of a high degree of variability across the A. radicans categories (Table 4).

Table 4 Plant importance values under different surveyed years (mean ± SD) (%)

Effects of Acmella radicans on species diversity

The results demonstrated that the first and second years when A. radicans was present exhibited the highest species richness (17.00 and 16.75), Simpson index (0.864 and 0.854), Shannon-Wiener index (2.196 and 2.205), and Pielou index (0.773 and 0.775), followed by the subsequent three years of the A. radicans infestation (Table 5). Overall, species richness, diversity, and evenness parameters very rarely exhibited significant differences within the first two years. Over time, species richness, diversity, and evenness values significantly declined indicating that invasion of A. radicans resulted in decreased species richness, diversity indices, and evenness index of local plant species communities.

Table 5 Diversity indices of plant communities for different surveyed years (mean ± SD)

Effects of Acmella radicans on soil characteristics

Among different A. radicans cover classes, soil characteristics varied significantly (Table 6). As the % cover of A. radicans increased, the pH, total K, and available K of the soil significantly increased; however, the organic matter, total N, total P, available N, and available P content of the soil declined significantly. It is evident that A. radicans may change the soil environment through absorbing more soil nutrients than other plants, which may facilitate its growth and invasion.

Table 6 Soil properties (i.e., pH, organic matter, total N, total P, total K, available N, available P, and available K) of different Acmella radicans % cover classes (mean ± SD)

Discussion

Acmella radicans is known to be a serious invasive plant in Bangladesh, Cuba, Curaçao, India, Tanzania, and Thailand [17,18,19, 36], but is rarely reported in China. We were the first to record the occurrence of A. radicans in Yunnan Province, China. Our field survey found that the total area infested by A. radicans was over 5560 ha in Yunnan, primarily distributed in farmland, orchard land, roadsides, and wasteland in Baoshan City and Lincang City. Our MaxEnt model revealed that the potentially suitable habitat in China for A. radicans was concentrated in southern China. The extent of areas rated as high and moderate suitability habitat areas for A. radicans in China were 5425 and 26,338 km2, accounting for 0.06 and 0.27% of the Chinese mainland area, respectively. The survival and expansion of invasive plant species are often favored by high rainfall and temperature levels in invaded regions [6, 8]. The study area is characterized by a subtropical climate with heavy rainfall and high temperature, which is suitable for plant growth of A. radicans. Moreover, the range of habitat suitability would likely increase under climate change.

The MaxEnt model approach has been extensively applied in species distribution simulation studies [37, 38] due to its advantages of being less sensitive to covariates in environmental variables [39], remaining stable with small sample sizes [40], and being able to fit complex variable relationships [41]. However, since the complex functional relationships in the MaxEnt model tend to lead to overfitting, the MaxEnt model is usually optimized by adjusting the parameters [42, 43]. Two important parameters that affect the MaxEnt model are the RM and FC. Currently, the ENMeval package [44], the SDMtune package [45], the Kuenm package [46], Wallace, Dismo and other open source software packages are effective in ascertaining the best parameter settings for the MaxEnt model. The Kuenm package allows MaxEnt to be used to establish a detailed and repeatable niche model, providing detailed model selection and calibration unavailable in other data package options [47]. In this study, we used the Kuenm package to find the optimal parameter settings of the MaxEnt model based on the mission rate at 5% and delta AICc. The fact the AUC value of the optimal model of 0.974 exceeded 0.9 indicated that the prediction performance of the MaxEnt model was improved by finding the optimal parameter settings, and that the possibility of over fitting the model was weakened.

It is important to evaluate potential spread of invasive plants based on local landscape factors [48]. Based on our observations in its present range in China, habitats A. radicans tends to invade include farmland, orchard land, roadsides, and wasteland. Throughout both moderately and highly suitable regions in its potential range in China, the landscape includes extensive areas of farmland producing similar crops to those found in its current range. Thus, the land area deemed moderately and highly suitable by our modelling is also vulnerable to the same kinds of impacts on plant community diversity and soil quality. Mountainous areas and other areas with extensive forest in the potential invaded range are likely the most serious barriers to spread; however, given the ability of A. radicans to grow on roadsides its ability is to invade these areas is highly probable as well. If A. radicans is likely to spread further in Yunnan and other parts of southern China as the SDM modelling predicts, it is also important to evaluate the seriousness of the impact on invaded plant communities. Invasive alien plants can reduce species richness, species diversity indices, evenness index, and native community structure and function [9, 49, 50]. The tendency of invasive plants to rapidly gain a competitive advantage over neighboring plants is attributed to their rapid growth, reproductive propagation ability, adaptation to a broad range of habitats, and strong allelopathic profiles [51,52,53,54]. By the second year of our five-year field survey, A. radicans had become the dominant and most densely populated species within the plant communities we monitored. Population density and importance values of many plant species, i.e., (A) conyzoides, (B) pilosa, B. latifolia, (C) virgata, C. benghalensis, C. dactylon, (D) sanguinalis, (E) hispidula, E. indica, and S. plicata declined significantly as A. radicans increased over time. Our field study found that the plant biomass, root number and depth, stem length, leaf area, and seed production number of A. radicans were markedly higher than those of neighboring plants, showing A. radicans achieved its competitive advantage via its extensive roots and larger aboveground size. Moreover, A. radicans may also have potential allelopathic compounds that could result in the inhibition of seed germination and seedling growth in neighboring plants [55]. Clearly, A. radicans is a strong competitor in this environment.

In addition to outcompeting neighboring plants, invasive alien plants can alter plant community structure and function in infested areas [9, 49]. Species richness and diversity of local plant communities was reduced 41% and 16% by L. camara invasion, respectively [56]. Ageratina adenophora caused a 68% reduction in species richness of understory vegetation in a Pinus yunnanensis forest, displacing countless native species in the process [57]. Similarly, the presence of M. micrantha was shown to decrease plant community parameters such as species richness and diversity [58]. Our results showed that the longer A. radicans was present in a plant community, the more species richness, diversity, and evenness values declined, likely threatening the stability of the system [59]. Thus, management of A. radicans is required in infested areas before it is too late.

Increasing numbers of studies have suggested that some invasive alien plants can effectively alter soil conditions, modifying soil chemistry, nutrient content, and nutrient availability, modifications which may actually lead to further invasion [60,61,62]. Soil nutrient concentrations have been observed to either increase or decrease in invaded habitats as compared with uninvaded habitats, depending on various factors [63,64,65]. For example, organic matter and concentrations of nitrogen, phosphorus, and potassium were lower in Parthenium hysterophorus-invaded soil than in non-invaded sites [64]. Mikania micrantha can likewise deplete soil nutrients, when growing either in monoculture or mixed culture with other plant species [50]. Our current study showed that the pH, total K, and available K in the soil increased significantly with increasing percent cover of A. radicans, but the organic matter, total N, total P, available N, and available P content of the soil declined significantly. Clearly A. radicans has the potential to change the soil environment through absorbing more soil nutrients, which may facilitate its growth and invasion.

Conclusions

In conclusion, our results indicated that the MaxEnt model fit (AUC = 0.974) was robust and the most important environmental variables affecting the potential suitability of a habitat included precipitation of the wettest month, isothermality, topsoil TEB, and precipitation seasonality. Within China, potential suitable areas for A. radicans are mainly located in southern regions. Most suitable habitats for A. radicans are within Yunnan, whereas the moderately suitable habitats are distributed in parts of Yunnan as well as in some other provinces in southern China. In addition to evaluating the ability of A. radicans to expand its range, we looked at its potential impacts on invaded plant communities by evaluating species diversity within the current invaded range. Species richness, species diversity, and evenness of local communities were also reduced with increasing duration of A. radicans invasion. The importance values and population density of most plant species as well as species richness, diversity, and evenness values of local communities decreased with increasing duration of A. radicans infestation. Soil nutrients declined markedly with increasing A. radicans cover percentage. Hence, our study is a wakeup call for managers to be more aware of this species and develop increased monitoring and mitigation efforts before A. radicans becomes more widely distributed. Furthermore, given the relatively few studies on the biology and impacts of this species, we recommend pursuing more detailed research on the competitive mechanisms that result in the reduced plant community diversity observed in southern China. Because there is evidence that A. radicans produces phytochemicals harmful to other plants (Shen et al., unpublished), it would be particularly useful to research possible allelopathic mechanisms.

Materials and methods

Study site

Baoshan City (24°08′ N-25°51′ N, 98°05′ E-100°02′ E) and Lincang City (23°05′ N-25°03′ N, 98°40’-100°32’E) are neighboring and located in the west boundary area of Yunnan Province, Southwest China (Fig. 4). Baoshan shares a boundary of 167.78 km with Burma in the south and northwest, and Lincang shares a boundary of 290.79 km with Burma in the southwest [66, 67]. The study area is part of the “Gold Cross” for biological diversity of the earth [68, 69]. Baoshan and Lincang belong to the “Bridgehead-Golden Port” of Chinese southern Indian Ocean strategy and important “windows and gateways” of performance of reform and open policy along the frontier [70]. Two cities have a same subtropical low-latitude mountain plateau monsoon climate, characterized by warm winter, cool summer, abundant rainfall, and wet and dry seasons. The annual average temperature of Baoshan is 14–17 °C and the annual rainfall is 746.6–2095.2 mm. The annual average temperature and annual rainfall of Lincang are 16.5–19.6 °C and 1485.7 mm, respectively [66, 67]. Due to particular geographical and ecological conditions, Baoshan and Lincang are suitable for the growth of many tropical and subtropical invasive alien plants, including some new invasive species entering in Yunnan through long border with Burma.

Fig. 4
figure 4

Distribution sampling site of Acmella radicans in Yunnan

Study species

Acmella radicans, is an annual erect herb in the Asteraceae family, native to Central America [17]. Some plant characteristic of A. radicans such as plant habit, soil surface seed bank, inflorescence, stem, and root system are shown in Fig. 5. Flowering and fruiting of A. radicans occurs between November and March and each plant can produce up to 14,300 seeds. This plant prefers moist habitats including riparian areas, ditches, and cultivated fields. Native to Central America introduced populations of A. radicans have been found in Bangladesh, Cuba, Curaçao, India, Tanzania, and Thailand [17,18,19, 36]. This plant is frequently used as a vegetable and in traditional medicine, particularly in its native range in meso-America [17, 55, 71]. Populations of A. radicans typically grow in riparian zones along streams producing many seeds, fostering continued range expansion through seed dispersal by water [20].

Fig. 5
figure 5

The plant habit (a, b), soil surface seed bank (c), inflorescence (d), stem (e), and root system (f) of Acmella radicans

Occurrence data ofAcmella radicans.

In our field survey, 12 new distribution records were obtained in Yunnan Province, Southwest China (Fig. 4). Their specific distribution locations are listed in Table 7.

Table 7 Distribution locations of Acmella radicans in Yunnan

Data for 759 A. radicans distribution records were collected from the Global Biodiversity Information Facility (http://www.gbif.org/), the Chinese Virtual Herbarium (http://www.cvh.ac.cn/), and China National Knowledge Infrastructure (https://www.cnki.net/), including 5 distribution records from eastern China. The aforementioned data was combined with the 12 distribution records obtained from our field survey, and data considered invalid were removed. Additionally, to reduce model overfitting and improve model prediction precision, redundant data were removed using the function of removing duplicate occurrences (RDOc) in ENMtools [72] to ensure that only a single distribution point existed in the same grid, which resulted in 396 distribution records available for MaxEnt modeling.

Bioclimatic and soil variables

The initial pool of environmental variables included 19 bioclimatic variables, 1 topographic variable (elevation), and 14 soil variables. The 14 soil variables comprised topsoil fractions for gravel, sand, silt, and clay as well as topsoil USDA texture classification, topsoil reference bulk density, topsoil organic carbon, topsoil pH (H2O), topsoi1 CEC (clay), topsoil CEC (soil), topsoil base saturation, topsoil TEB, topsoil sodicity (ESP), and topsoil salinity (Elco). Global bioclimatic variables and elevation data were from WorldClim (https://worldclim.org/) using Historical climate data (1970–2000), and global soil data from HWSD (https://www.fao.org/soils-portal/soil-survey/soil-maps-and-databases/harmonized-world-soil-database-v12/en/). All environmental variables were at a spatial resolution of 2.5 min. We utilized the world administrative map from the National Earth System Science Data Center, National Science and Technology Infrastructure of China.

To avoid overfitting the model predictions due to covariance among environmental variables, correlation analysis of environmental variables was required prior to using them to construct our ecological niche model [73]. Therefore, we did Pearson correlation tests in SPSS 27 for the bioclimatic variables and soil variables (Fig. 6). The environmental variables with higher contribution rates were retained in combination with the Jackknife test results in the MaxEnt model, and 19 environmental variables from the initial pool were selected for MaxEnt modeling (Table 8).

Table 8 Variables used in the model prediction
Fig. 6
figure 6

Pearson correlation coefficient. a Correlation coefficient matrix of 19 climate variables. b Correlation coefficient matrix of 14 soil variables

MaxEnt model calibration, settings, and evaluation

The prediction performance of MaxEnt is influenced by two parameters: the regularization multiplier (RM) and feature combination (FC). We used the kuenm package in R (version 3.6.3) to adjust the parameters to achieve the detailed calibration, parameter selection, candidate model evaluation, and final model establishment of MaxEnt [46]. The RM value was set from 0.5 to 4.0, with each increment of 0.5, for a total of 8 selections. Additionally, 29 FCs of 5 features, including linear (L), quadratic (Q), product (P), threshold (T), and hinge (H), were selected for testing, including L, Q, P, T, H, LQ, LP, LT, LH, QP, QT, QH, PT, PH, TH, LQP, LQT, LQH, LPT, LPH, QPT, QPH, QTH, PTH, LQPT, LQPH, LQPH, LQTH, LPTH, LPTH, and LQPTH. Kuenm selected significant models from 232 candidate models with an omission rate of ≤ 5%, and then selected the model with delta AICc < 2 as the optimal model according to akaike information criterion (AICc) [46].

After parameterization according to the optimal model evaluated by kuenm, the species distribution point data and 19 environmental variables were imported by MaxEnt (version 3.4.4). A random sample of 75% of the distribution was selected as the training data set for modeling, and the remaining 25% of the distribution was used as the test data set to validate the model. The maximum number of background points was set to 10,000 and the maximum number of iterations to 500. The calculation results were then repeated 10 times using the Bootstrap method for the final average and we chose the Logistic output type. In the environment parameter setting, we evaluated the weight of each environmental variable using the Jackknife method and determined the dominant environment variable by combining the contribution rate and replacement important value of each environment variable. We used the area under curve (AUC) value of the receiver operating characteristic curve (ROC) to reflect the combined model sensitivity level and specificity [74]. An AUC value exceeding 0.8 indicates that the model predicts well, and an AUC value exceeding 0.9 indicates that the model predicts very well [75].

Finally, based on the global suitability threshold of A. radicans predicted by the MaxEnt model, and Janks natural breakpoint method was used to divide the suitability index of A. radicans into four classes: unsuitable (0 < P < 0. 05), slightly suitable area (0.05 < P < 0. 18), moderately suitable area (0.18 < P < 0. 41), and highly suitable area (P > 0.41), and obtained the potential distribution of A. radicans within China.

Ecological impact assessment of Acmella radicans

To explore the ecological impacts of A. radicans on local plant communities and soil nutrient characteristics, thirty 1 m×1 m quadrats were randomly selected at each kind of habitats with four replicates in the same village in Changning County every year for five years since 2018. Climate and altitude did not differ among quadrats. Thus, a total of six hundred quadrats for five different habitats were surveyed each year for five years (2018 - year one up to 2022 – year five). Information collected from each quadrat consisted of plant species, plant density plant cover, frequency, and plant height for all species. Meanwhile, fifty soil samples were taken randomly from different A. radicans cover percentages (0%, 1–25%, 26–50%, 51–75%, and 76–100%) at the last field surveyed (2022) and then combined and treated as a single composite sample. The soil samples were then ground up, sifted through a 2 mm sieve, and air-dried at room temperature in preparation to analyze soil chemical characteristics. The pH, soil organic matter, total N, total P, total K, available N, available P, available K of soil nutrients were examined at the Soil Analysis and Detection Center of the Agricultural Environment and Resource Research Institute, Yunnan Academy of Agricultural Sciences, China.

Lastly, plant density, importance values and plant diversity parameters were calculated. The importance value (IV) for a given species was calculated as IV = (relative density + relative cover + relative frequency)/3 (Relative values were obtained via dividing species specific values by the sums of the densities, cover proportions and frequencies of all species in a plot, respectively) [52]. Simpson diversity index (D) [76] was calculated as D = 1-∑[Ni(Ni − 1)/N(N − 1)], with Ni being the total number of individuals from i species, N the total number of individuals for all species in a plot. D values range from 0 to 1, from lowest to highest diversity. Shannon-Wiener diversity index (H) [77] was calculated as H = -∑PilnPi, where Pi is the proportion of species i relative to the total number of species per plot. Pielou evenness index (J) [77] was calculated as J = H/lnS, where S is the species richness of each plot.

An analysis of variance (one-way ANOVA) was used to compare treatments. Duncan’s multiple range tests were used to compare treatments at a 5% significance level.

Data availability

The raw data supporting the conclusions of this article will be made available by the corresponding author Fudou Zhang, without reservation.

References

  1. Axmacher JC, Sang W. Plant invasions in China – challenges and chances. PLoS ONE. 2013;8:e64173.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  2. Richardson DM, Pyšek P, Rejmánek M, Barbour MG, Panetta FD, West CJ. Naturalization and invasion of alien plants: concepts and definitions. Divers. Distrib. 2000;6:93–107.

    Google Scholar 

  3. Ricciardi A. Assessing species invasions as a cause of extinction. Trends Ecol Evol. 2004;19:619.

    Article  Google Scholar 

  4. Porté AJ, Lamarque LJ, Lortie CJ, Michalet R, Delzon S. Invasive Acer negundo outperforms native species in non-limiting resource environments due to its higher phenotypic plasticity. BMC Ecol. 2011;11:28.

    Article  PubMed  PubMed Central  Google Scholar 

  5. Rai RK, Shrestha L, Joshi S, Clements DR. Biotic and economic impacts of plant invasions. In: Clements DR, editor. Global plant invasions. Springer Nature Switzerland AG; 2022. pp. 301–15.

  6. Hejda M, Pyšek P, Pergl J, Sádlo J, Chytrý M, Jarošík V. Invasion success of alien plants: do habitats affinities in the native distribution range matter? Global Ecol. Biogeogr. 2009;18:372–82.

    Article  Google Scholar 

  7. Hulme PE. Trade, transport and trouble: managing invasive species pathways in an era of globalization. J Appl Ecol. 2009;46:10–8.

    Article  Google Scholar 

  8. Powell KI, Chase JM, Knight TM. A synthesis of plant invasion effects on biodiversity across spatial scales. Am J Bot. 2011;98:539–48.

    Article  PubMed  Google Scholar 

  9. Vilà M, Espinar JL, Hejda M, Hulme PE, Jarošík V, Maron JL, et al. Ecological impacts of invasive alien plants: a meta-analysis of their effects on species, communities and ecosystems. Ecol Lett. 2011;14:702–8.

    Article  PubMed  Google Scholar 

  10. Keane RMR, Crawley MJM. Exotic plant invasions and the enemy release hypothesis. Trends Ecol Evol. 2002;17:164–70.

    Article  Google Scholar 

  11. Ortega YK, Pearson DE. Weak vs. strong invaders of natural plant communities: assessing invasibility and impact. Ecol Appl. 2005;15:651–61.

    Article  Google Scholar 

  12. Theoharides KA, Dukes JS. Plant invasion across space and time: factors affecting nonindigenous species success during four stages of invasion. New Phytol. 2007;176:256–73.

    Article  PubMed  Google Scholar 

  13. Shrestha BB, Witt ABR, Shen S, Khuroo AA, Shrestha UB, Naqinezhad A. Plant invasions in Asia. In: Clements DR, editor. Global plant invasions. Springer Nature Switzerland AG; 2022. pp. 89–127.

  14. Xie Y, Li Z, William PG, Li D. Invasive species in China – an overview. Biodivers Conserv. 2001;10:1317–41.

    Article  Google Scholar 

  15. Jiang H, Fan Q, Li J, Shi S, Li S, Liao W, et al. Naturalization of alien plants in China. Biodivers Conserv. 2011;20:1545–56.

    Article  Google Scholar 

  16. Fang Y, Zhang X, Wei H, Wang D, Chen R, Wang L, et al. Predicting the invasive trend of exotic plants in China based on the ensemble model under climate change: a case for three invasive plants of Asteraceae. Sci Total Environ. 2021;756:143841.

    Article  CAS  PubMed  Google Scholar 

  17. Bagga J, Deshmukh UB. Acmella radicans (Jacquin) R.K. Jansen (Asteraceae) – A new distributional plant record for Jharkhand State (India). J New Biol Rep. 2018;7:24–27.

  18. Rahman MM, Khan SA, Hossain GM, Jakaria M, Rahim MA. Acmella radicans (Jacq.) R.K. Jansen (Asteraceae) – a new angiosperm record. J Biol Sci. 2016;5:87–93.

  19. Maity D, Sardar A, Dash SS. Acmella radicans (Asteraceae), and American weed new to Eastern India. Nelumbo. 2017;59:54–7.

    Article  Google Scholar 

  20. Wang Z, Yan X, Li H, Ma J. Acmella radicans var. Debilis (Kunth) R.K. Jansen (Asteraceae), a newly naturalized plant in China. J Trop Subtrop Bot. 2015;23:643–6.

    Google Scholar 

  21. Qin Z, Zhang J, DiTommaso A, Diez JM, Zhao Y, Wang F. Predicting the potential distribution of three allergenic invasive Ambrosia (ragweed) species in Asia. J Environ Inf. 2022;391:49–66.

    Google Scholar 

  22. Phillips SJ, Anderson RP, Schapire RE. Maximum entropy modeling of species geographic distribution. Ecol Model. 2006;190:231–59.

    Article  Google Scholar 

  23. Cheng L, Gao J. Prediction of spatial distribution of Eupatorium adenophorum Sprengel based on GARP model: a case study in Longitudinal Range-Gorge Region of Yunnan Province. Chin J Ecol. 2008;27:1531–6.

    Google Scholar 

  24. Botkin DB, Saxe H, Araujo MB, Betts R, Bradshaw RHW, Cedhagen T, et al. Forecasting the effects of global warming on biodiversity. Bioscience. 2007;57:227–36.

    Article  Google Scholar 

  25. Poutsma J, Loomans AJ, Aukema B, Heijerman T. Predicting the potential geographical distribution of the harlequin ladybird, Harmonia axyridis, using the CLIMEX mode. Biocontrol. 2008;53:103–25.

    Article  Google Scholar 

  26. Zhang Y, Tang J, Ren G, Zhao K, Wang X. Global potential distribution prediction of Xanthium italicum based on Maxent model. Sci Rep. 2021;11:16545.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  27. Yan H, Feng L, Zhao Y, Feng L, Wu D, Zhu C. Prediction of the spatial distribution of Alternanthera philoxeroides in China based on ArcGIS and MaxEnt. Glob Ecol Conserv. 2020;21:e00856.

    Google Scholar 

  28. Yan H, Feng L, Zhao Y, Feng L, Zhu C, Qu Y, et al. Predicting the potential distribution of an invasive species, Erigeron canadensis L. in China with a maximum entropy mode. Glob Ecol Conserv. 2020;21:e00822.

    Google Scholar 

  29. Bushi D, Mahato R, Nimasow OD, Nimasow G. MaxEnt-based prediction of the potential invasion of Lantana camara L. under climate change scenarios in Arunachal Pradesh, India. Acta Ecol Sin. 2022;34:674–83.

    Google Scholar 

  30. Phillips SJ, Dudík M, Elith J, Graham CH, Lehmann A, Leathwick J, et al. Sample selection bias and presence-only distribution models: implications for background and pseudo-absence data. Ecol Appl. 2009;19:181–97.

    Article  PubMed  Google Scholar 

  31. Warton DI, Renner IW, Daniel R. Model-based control of observer bias for the analysis of presence-only data in ecology. PLoS ONE. 2013;8:e79168.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  32. Fourcade Y, Engler JO, Rödder D, Secondi J. Mapping species distributions with MAXENT using a geographically biased sample of presence data: a performance assessment of methods for correcting sampling bias. PLoS ONE. 2014;9:e97122.

    Article  PubMed  PubMed Central  Google Scholar 

  33. Jiménez-Valverde A, Lobo JM, Hortal J. Not as good as they seem: the importance of concepts in species distribution modeling. Divers Distrib. 2008;14:885–90.

    Article  Google Scholar 

  34. Westbrook AS, Nikkel E, Clements DR, DiTommaso A. Modeling and managing invasive weeds in a changing climate. In: Ziska LH, editor. Invasive species and global climate change. GB: CABI; 2023. pp. 282–306.

    Google Scholar 

  35. Lobo JM, Jiménez-Valverde A, Real R. AUC: a misleading measure of the performance of predictive distribution models. Global Ecol Biogeogr. 2008;17:145–51.

    Article  Google Scholar 

  36. Panyadee P, Inta A. Taxonomy and ethnobotany of Acmella (Asteraceae) in Thailand. Biodiversitas. 2022;23:2177–86.

    Article  Google Scholar 

  37. Holder AM, Markarian A, Doyle JM, Olson JR. Predicting geographic distributions of fishes in remote stream networks using maximum entropy modeling and landscape characterizations. Ecol Model. 2020;433:109231.

    Article  Google Scholar 

  38. Alatawi AS, Gilbert F, Reader T. Modelling terrestrial reptile species richness, distributions and habitat suitability in Saudi Arabia. J Arid Environ. 2020;178:104153.

    Article  Google Scholar 

  39. Elith J, Phillips SJ, Hastie T, Dudík M, Chee YE, Yates CJ. A statistical explanation of Maxent for ecologists. Divers Distrib. 2011;17:43–57.

    Article  Google Scholar 

  40. Guisan A, Zimmermann NE, Elith J, Graham CH, Phillips S, Peterson AT. What matters for predicting the occurrences of trees: techniques, data, or species. Characteristics? Ecol Monogr. 2007;77:615–30.

    Article  Google Scholar 

  41. Elith J, Graham CH, Anderson RP, Dudik M, Ferrier S, Guisan A, et al. Novel methods improve prediction of species’ distributions from occurrence data. Ecography. 2014;29:129–51.

    Article  Google Scholar 

  42. Radosavljevic A, Anderson RP. Making better Maxent models of species distributions: complexity, overfitting and evaluation. J Biogeogr. 2014;41:629–43.

    Article  Google Scholar 

  43. Morales NS, Fernández IC, Baca-González V. MaxEnt’s parameter configuration and small samples: are we paying attention to recommendations? A systematic review. PeerJ. 2017;5:e3093.

    Article  PubMed  PubMed Central  Google Scholar 

  44. Muscarella R, Galante PJ, Soley-Guardia M, Boria RA, Kass JM, Uriarte M, et al. ENMeval: an R package for conducting spatially independent evaluations and estimating optimal model complexity for Maxent ecological niche models. Methods Ecol Evol. 2014;5:1198–205.

    Article  Google Scholar 

  45. Vignali S, Barras AG, Arlettaz R, Braunisch V. SDMtune: an package to tune and evaluate species distribution models. Ecol Evol. 2020;10:11488–506.

    Article  PubMed  PubMed Central  Google Scholar 

  46. Cobos ME, Peterson AT, Barve N, Osorio-Olvera L. Kuenm: an R package for detailed development of ecological niche models using Maxent. PeerJ. 2019;7:e6281.

    Article  PubMed  PubMed Central  Google Scholar 

  47. Escobar LE. Ecological niche modeling: an introduction for veterinarians and epidemiologists. Front Vet Sci. 2020;7:519059.

    Article  PubMed  PubMed Central  Google Scholar 

  48. Nikkel E, Clements DR, Anderson D, Williams JL. Regional habitat suitability for aquatic and terrestrial invasive plant species may expand or contract with climate change. Biol Invasions. 2023;25:3805–3822.

    Article  PubMed  PubMed Central  Google Scholar 

  49. Ehrenfeld JG. Ecosystem consequences of biological invasions. Annu Rev Ecol Evol Syst. 2010;41:59–80.

    Article  Google Scholar 

  50. Shen S, Xu G, Clements DR, Jin G, Chen A, Zhang F, et al. Suppression of the invasive plant mile-a-minute (Mikania micrantha) by local crop sweet potato (Ipomoea batatas) by means of higher growth rate and competition for soil nutrient. BMC Ecol. 2015;15:1.

    Article  PubMed  PubMed Central  Google Scholar 

  51. De Cauwer B, Devos R, Claerhout S, Bulcke R, Reheul D. Seed dormancy, germination, emergence and seed longevity in Galinsoga parviflora and G. Quadriradiata. Weed Res. 2014;54:38–47.

    Article  Google Scholar 

  52. Shen S, Xu G, Clements DR, Jin G, Liu S, Zhang F, et al. Effects of invasive plant Mikania micrantha on plant community and diversity in farming systems. Asian J Plant Sci. 2015;14:27–33.

    Article  Google Scholar 

  53. Damalas CA. Distribution, biology, and agricultural importance of Galinsoga parviflora (Asteraceae). Weed Biol Manag. 2008;8:147–53.

    Article  Google Scholar 

  54. Shen S, Xu G, Li D, Yang S, Jin M, Liu S, et al. Adventitious roots support population expansion of the invasive plant Mikania micrantha Kunth. Physiol Plant. 2021;173:911–9.

    Article  CAS  PubMed  Google Scholar 

  55. Jirovetz L, Buchbauer G, Abraham GT, Shafi MP. Chemical composition and olfactoric characterization of Acmella radicans (Jacq.) R.K. Jansen var. Radicans from southern India. Flavour Frag J. 2006;21:88–91.

    Article  CAS  Google Scholar 

  56. Singh HP, Batish DR, Dogra KS, Kaur S, Kohli RK, Negi A. Negative effect of litter of invasive weed Lantana camara on structure and composition of vegetation in the lower Siwalik Hills, northern India. Environ Monit Assess. 2014;186:3379–89.

    Article  PubMed  Google Scholar 

  57. Fu D, Wu X, Huang N, Duan C. Effects of the invasive herb Ageratina adenophora on understory plant communities and tree seedling growth in Pinus yunnanensis forests in Yunnan, China. J Res. 2018;23:112–9.

    Article  CAS  Google Scholar 

  58. Kaur R, Malhotra S, Inderjit. Effects of invasion of Mikania micrantha on germination of rice seedlings, plant richness, chemical properties and respiration of soil. Biol Fert Soils. 2012;48:481–8.

    Article  CAS  Google Scholar 

  59. Tilman D, Reich PB, Knops JM. Biodiversity and ecosystem stability in a decade-long grassland experiment. Nature. 2006;441:629–32.

    Article  CAS  PubMed  Google Scholar 

  60. Callaway RM, Newingham B, Zabinski CA, Mahall BE. Compensatory growth and competitive ability of an invasive weed are enhanced by soil fungi and native neighbours. Ecol Lett. 2001;4:429–33.

    Article  Google Scholar 

  61. Huang Q, Li X, Huang F, Wang R, Lu B, Shen Y, et al. Nutrient addition increases the capacity for division of labor and the benefits of clonal integration in an invasive plant. Sci Total Environ. 2018;643:1232–8.

    Article  CAS  PubMed  Google Scholar 

  62. Shen S, Xu G, Yang Y, Yu X, Li D, Yang S, et al. Increased suppressive effect of Ipomoea batatas (sweet potato) on Mikania micrantha (mile-a-minute) under high fertilization levels. Manag Biol Invasion. 2020;11:560–75.

    Article  Google Scholar 

  63. Kaur A, Kaur S, Singh HP, Batish DR, Kohli RK. Phenotypic variations alter the ecological impact of invasive alien species: lessons from Parthenium hysterophorus. J Environ Manag. 2019;241:187–97.

    Article  CAS  Google Scholar 

  64. Timsina B, Shrestha BB, Rokaya MB, Münzbergová Z. Impact of Parthenium hysterophorus L. invasion on plant species composition and soil properties of grassland communities in Nepal. Flora. 2011;206:233–40.

    Article  Google Scholar 

  65. Zhao M, Lu X, Zhao H, Yang Y, Hale L, Gao Q, et al. Ageratina adenophora invasions are associated with microbially mediated differences in biogeochemical cycles. Sci Total Environ. 2019;677:47–56.

    Article  CAS  PubMed  Google Scholar 

  66. Long Z, Gao C. On spatial distribution characteristics and influencing factors of the hot spring tourist destinations in Baoshan City. J Southwest China Normal Univ (Nat Sci Ed). 2021;46:140–9.

    Google Scholar 

  67. Liu F, Zhang X, Ren J, Liu Y, Wang W. Construction of ecological security pattern based on ecosystem service evaluation and minimum resistance model: a case study of Lincang City. Environ Sci Surv. 2021;40:78–83.

    Google Scholar 

  68. Hou M, Zhu W, Li R, Qing H, Wei Y. Presenting and analysing of the Gold Cross in biological diversity of the earth. J Kunming Univ Sci Technol (Sci Technol). 2006;31:1–7.

    Google Scholar 

  69. Dong W, Guo X. The present status and conservation of species diversity in Yunnan Province. China Sci Tech Inf Process. 2008;14:19–23.

    Google Scholar 

  70. Zhang A. Reflect on illegal residence in Yunnan and Burma border area. J Yunnan Police Officer Acad. 2013;100:38–42.

    Google Scholar 

  71. Rios-Chavez P, Ramirez-Chavez E, Armenta-Salinas C, Molina-Torres J. Acmella radicans var. Radicans: in vitro culture establishment and alkamide content. Vitro Cell Dev Biol -Plant. 2003;39:37–41.

    Article  CAS  Google Scholar 

  72. Warren DL, Glor RE, Turelli M. ENMTools: a toolbox for comparative studies of environmental niche models. Ecography. 2010;33:607–11.

    Article  Google Scholar 

  73. Stuhldreher G, Fartmann T. Threatened grassland butterflies as indicators of microclimatic niches along an elevational gradient-implications for conservation in times of climate change. Ecol Indic. 2018;94:83–98.

    Article  Google Scholar 

  74. Gebrewahid Y, Abrehe S, Meresa E, Eyasu G, Abay K, Gebreab G, et al. Current and future predicting potential areas of Oxytenanthera abyssinica (A. Richard) using MaxEnt model under climate change in Northern Ethiopia. Ecol Process. 2020;9:6.

    Article  Google Scholar 

  75. Sun S, Zhang Y, Huang D, Wang H, Cao Q, Fan P, et al. The effect of climate change on the richness distribution pattern of oaks (Quercus L.) in China. Sci Total Environ. 2020;744:140786.

    Article  CAS  PubMed  Google Scholar 

  76. Simpson EH. Measurement of diversity. Nature. 1949;163:688.

    Article  Google Scholar 

  77. Ma K, Liu Y. Measurement of biotic community diversity. I. α diversity (part 2). Chin Biodivers. 1994;2:231–9.

    Google Scholar 

Download references

Acknowledgements

We wish to thank Pengfei Liu from the Agricultural Environmental Protection and Monitoring Station of Changning County for his tremendous field support.

Funding

This research was supported by the Key Research and Development Program of Yunnan Province (202103AF140007 and 202203AE140008), National Key R&D Program of China (2021YFC2600400), Ten Thousand Talent Program (Young Top-Notch Talent) of Yunnan Province (YNWR-QNBJ-2018-201), and Special Funds of Major Science and Technology Project in Yunnan Province (202102AE090003).

Author information

Authors and Affiliations

Authors

Contributions

SS conceived and designed the experiments, performed the experiments, analyzed the data, and wrote the draft of the manuscript. FZ (Fengping Zheng), WZ, GX, DL, SY, GJ, AC, YC, ZF, and LY performed the experiments. DC and EN analyzed the data and wrote the draft of the manuscript. FZ (Fudou Zhang) conceived and designed the experiments and performed the experiments. All authors reviewed the manuscript.

Corresponding authors

Correspondence to Lun Yin or Fudou Zhang.

Ethics declarations

Ethics approval and consent to participate

All the steps/ methods were in compliance with relevant institutional, national, and international guidelines and legislation.

Consent for publication

Not applicable.

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Rights and permissions

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://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Shen, S., Zheng, F., Zhang, W. et al. Potential distribution and ecological impacts of Acmella radicans (Jacquin) R.K. Jansen (a new Yunnan invasive species record) in China. BMC Plant Biol 24, 494 (2024). https://doi.org/10.1186/s12870-024-05191-5

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s12870-024-05191-5

Keywords