Next Article in Journal
Establishment of Axenic Cultures for Two Marine Microalgae and the Effects of Antibiotics on Their Growth
Previous Article in Journal
Impact of Extreme Weather Events on the Surface Energy Balance of the Low-Elevation Svalbard Glacier Aldegondabreen
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Response Patterns of Non-Native Rotifer Species Kellicottia bostoniensis (Rousselet, 1908) to Environmental Factors and Its Relationship with Potential Competitors and Predators

1
Department of Environmental Science and Engineering, Kyung Hee University, Yongin 17104, Republic of Korea
2
EcoResearch Inc., Gongju 32588, Republic of Korea
*
Authors to whom correspondence should be addressed.
Water 2025, 17(2), 273; https://doi.org/10.3390/w17020273
Submission received: 2 December 2024 / Revised: 9 January 2025 / Accepted: 17 January 2025 / Published: 19 January 2025

Abstract

:
Kellicottia bostoniensis is a rotifer species originating from North America, with numerous reports documenting its spread to Europe and South America. Meanwhile, in Asia, the occurrence of this non-native species has been scarcely documented, but its presence was recently reported in Korea in 2020. In Korea, the close geographic proximity of artificial reservoirs within a small area is expected to make the region vulnerable to the spread of K. bostoniensis. Therefore, it is deemed necessary to study the spatio-temporal distribution of K. bostoniensis, which appears to have extensively proliferated across the region, and its relationship with environmental factors. We analyzed the response of K. bostoniensis occurrence to environmental factors and aimed to compare its response patterns with those of potential competitor and predator groups, which were classified based on a functional group approach. This study, conducted on 12 small lakes located under similar climatic conditions, revealed that the response of K. bostoniensis to variables representing the trophic state of the lakes differed from those of the potential competitor and predator groups. The study sites, encompassing a wide range of trophic states from oligotrophic to hypereutrophic, across various lakes and survey periods, were well-suited for examining these relationships. Furthermore, when TSI was calculated to compare the abundance of K. bostoniensis and potential competitors across trophic conditions, the abundance of potential competitors was proportional to TSI. In contrast, K. bostoniensis exhibited higher abundances in mesotrophic and lower eutrophic environments, which distinguished it from its competitors. Our results indicate that the impact of introduced K. bostoniensis on native communities is minimal, while the resistance of native communities through the biological filter based on biological interactions is ineffective against this non-native rotifer species. This provides key evidence for explaining the widespread distribution of K. bostoniensis across continents.

1. Introduction

Non-native species are important at the community level because they directly or in-directly influence native species through biological interactions [1]. In invasive species studies, it is essential to observe their spreading across the new regions and to investigate the biotic and abiotic interactions that occur in these new environments. It is most of all necessary to monitor the abundance and distribution of such species and its dynamics, and to compare it with native species to see if the potentially invasive species causes extinction of local groups of organisms, especially taxonomically close ones [2,3,4]. Furthermore, information on the interactions between non-native species and the biotic and abiotic factors of their habitats provides valuable insights for predicting the likelihood of further introductions into new habitats and assessing potential changes to existing ecosystems [5].
Kellicottia bostoniensis (Rousselet, 1908) belongs to a freshwater zooplankton species (Rotifera, Monogononta, Ploima, Brachionidae), originating from North America. Since 1943, its presence has been reported in many European countries [6,7], and subsequently, its dispersion was also noted in South America during the 1990s [8]. In Asia, the occurrence of the genus Kellicottia was historically reported to be restricted to K. longispina until relatively recently [9]. In East Asia, although environmental conditions appear suitable for the establishment of K. bostoniensis, it is presumed that its population has not yet stably established itself [10]. In the East Asia in Palearctic realm, the occurrence of K. bostoniensis was first reported as a rare species in Japan in 1999 [11]. However, no further occurrences of K. bostoniensis have been recorded since then, and research on K. bostoniensis as a non-native species in Asia remains limited compared to studies conducted in North America, Europe, and South America.
The second recorded occurrence of K. bostoniensis in Asia was recently confirmed in Korea [12]. In Korea, artificial reservoirs constructed for agricultural and drinking water supply purposes are distributed nationwide, with many located in geographical proximity. This geographic arrangement suggests that human activities may continuously facilitate the spread of non-native species. Furthermore, many reservoirs are eutrophicated [13], and the lack of significant variation in geographical and climatic factors may result in low spatiotemporal heterogeneity of environmental factors within and between habitats. This can facilitate the spread of introduced non-native species [14]. Considering the temporal gaps between studies on the occurrence of K. bostoniensis in Asia and the limited information which increases the potential for misidentification with the morphologically similar K. longispina, there appears to be a significant lack of data regarding the distribution of K. bostoniensis and its interactions with other species within zooplankton communities in Asia, including Korea. Therefore, research on the spatio-temporal distribution of this species at the community level and its response patterns to environmental factors is considered necessary.
From the perspective of biological interactions, K. bostoniensis has a relatively long body length, including its spines (250–450 μm). It possesses long anterior and posterior spines, which lengthen in response to predation pressure, serving as a defensive mechanism against predators [15,16,17]. Thus, it is known to effectively respond to predation pressure from macro-invertebrate predators (e.g., Chaoborus larvae) and predatory rotifer, Asplanchna spp. [15,18]. Additionally, bacteria, nanoflagellates, and nano- and picophytoplankton have been studied as its primary food sources [19,20,21]. From the perspective of interactions with environmental factors, variables such as water temperature, dissolved oxygen concentration, and the trophic state index (TSI) of lakes have been identified as related to K. bostoniensis [14,15,17,22]. However, K. bostoniensis appears to respond differently to environmental factors depending on the study. For instance, Branco [15] described a strong negative correlation with water temperature, whereas Shurganova [23] reported a positive correlation with the same factor. Similarly, conflicting findings exist regarding the relationship between the trophic state of lakes and this species [14,17], suggesting that its responses vary depending on the complex environments and the range of environmental factors in the studied habitats. Based on current research findings, it is evident that further studies are needed to validate the relationships of K. bostoniensis with various biotic and abiotic factors across diverse environments. Moreover, studies interpreting population dynamics as the result of the combined effects of these factors remain insufficient.
In this study, a functional group-based approach was used to understand the population relationships between K. bostoniensis and other zooplankton species with which it interacts within the food web. Species with similar feeding characteristics within the functional group were identified as potential competitors, while species known as general predators of rotifers, including K. bostoniensis, were identified as potential predators. The response patterns of these groups to environmental factors were compared with those of K. bostoniensis. Additionally, the spatial distribution of K. bostoniensis within the Geumgang River basin in Korea, part of the East Asian region where regional studies on K. bostoniensis are limited, was investigated. This study aimed to elucidate the relationships between K. bostoniensis and environmental factors in a new environment, examine its fluctuations with potential competitor and predator groups, and explore the potential impact of this non-native rotifer species on native zooplankton communities.

2. Materials and Methods

2.1. Study Site and Data Collection

To investigate the spatial distribution of Kellicottia bostoniensis and its relationships with environmental factors and zooplankton community composition, a study was conducted on 12 small lakes (maximum surface area < 3 km2) located in close geographic proximity. These lakes, within a 50 km radius, belong to the same watershed (Geumgang River) and share similar physical characteristics and climatic conditions (Figure 1, Table A1). Surveys of zooplankton and measurements of environmental parameters were carried out in May, June, September, and November of 2023.

2.2. Data Collection

Zooplankton samples were collected using a plankton net (30 cm diameter, 60 μm mesh) at the center of each lake. Considering the vertical distribution of the zooplankton community, sampling was performed as follows: in lakes with depths more than 5 m, vertical towing covering the entire water column; and in lakes with depths less than 5 m, oblique towing was conducted from 5 m (towing speed: 1 m/s). The collected zooplankton samples were preserved with formalin (final concentration 5%), and a 1 to 5 mL subsample was extracted and examined in an S-R chamber under a microscope (Olympus BX51, Tokyo, Japan for identification to genus and species levels. Zooplankton abundance (density) was calculated by determining the volume of water filtered based on the net area and towing depth and was expressed as individuals per liter (ind./L).
Environmental factors that could influence zooplankton community were selected: water temperature (WT, °C), dissolved oxygen (DO, mg/L), suspended solids (SS, mg/L), total nitrogen (TN, mg/L), total phosphorus (TP, mg/L), total organic carbon (TOC, mg/L), chemical oxygen demand (COD, mg/L), pH, and chlorophyll-a concentration (Chl-a, mg/m3). Environmental factor data were obtained from publicly available datasets of Rural Agricultural Water Resource Information System (https://rawris.ekr.or.kr (accessed on 10 December 2024)), corresponding to the survey period of each lake (Table A2).

2.3. Classification of Functional Groups of Zooplankton

To distinguish species that directly interact with K. bostoniensis, potential competitor and predator species among the observed zooplankton were classified based on functional traits (Table 1, Table A3). The phylum Rotifera (referred to as rotifers), to which K. bostoniensis belongs, exhibits variations in feeding tendencies and consequent environmental responses depending on the structure of their mastax, specifically the trophi [24]. Species with malleate trophi, such as K. bostoniensis, feed by cutting and chewing their food. Rotifers can also be categorized by their habitat preferences into planktonic species (which live suspended in the water column) and littoral and periphytic species (which attach to substrates like aquatic plants or rocks) [25,26]. Accordingly, species which could compete with K. bostoniensis for food sources due to their malleate trophi and preference for planktonic habitat were classified as the potential competitor group. Meanwhile, the genus Asplanchna, which is a predatory rotifer, Asplanchna spp., and cyclopoid copepods, known as primary predators of rotifers [27,28,29], were classified as the potential predator group.

2.4. Statistical Methods

Differences in each environmental factors (environmental and biological factors) in the presence or absence of K. bostoniensis were analyzed using a permutation test (The R package “coin”) [30]. This method, a type of non-parametric statistical approach, mitigates issues related to sample size disparities and non-normality arising from differences in frequency between the presence and absence of K. bostoniensis. This analysis was conducted on nine environmental factors, zooplankton abundance by taxonomic group (rotifers, cladocerans, and copepods), and the abundance of potential competitor and predator groups.
Detrended Correspondence Analysis (DCA) was performed to analyze the response of K. bostoniensis to multivariate environmental factors and to compare these responses with those of potential competitor and predator groups (The R package “Vegan”) [31]. DCA was chosen to address non-linearities among variables caused by the absence data for K. bostoniensis (presence count: 0) [32].

3. Results

3.1. Seasonal Variation in the Abundance of Kellicottia bostoniensis, Zooplankton Factors, and Environmental Factors

Kellicottia bostoniensis was found in 7 out of 12 lakes, while it was absent in 5 lakes (Figure 1). For K. bostoniensis, the monthly average abundance was the lowest in May, with no occurrences recorded across all lakes.
Among the lakes where K. bostoniensis appeared, R-A-2, R-A-4, and R-A-8 exhibited very low abundance during most survey periods except for November, when the abundance of K. bostoniensis increased, showing similar seasonal appearance patterns (Figure 2(a1,a3,a6)). On the other hand, lakes R-A-3 and R-A-9 demonstrated the highest abundance in July, with lower occurrences during other survey periods (Figure 2(a2,a7)). Meanwhile, lakes R-A-6 and R-A-7 differed from other lakes by showing the highest abundance of K. bostoniensis in September, distinguishing seasonal appearance patterns by lake.
Potential competitors and predators also exhibited distinct seasonal occurrence patterns depending on the lake. For potential competitors, some lakes showed a pattern of increasing in July and then gradually decreasing (Figure 2(b2,b4,b6,b7)), while other lakes did not display a common seasonal pattern. In lakes R-A-2 and R-A-3, potential predators showed patterns similar to those of potential competitors’ patterns in many lakes, with high occurrences in July followed by gradual decreases (Figure 2(c1,c2)). Other patterns were further categorized into rapidly changing fluctuation patterns, with the lowest abundance in July (Figure 2(c3,c4)) and those with the highest abundance in July (Figure 2(c6,c7)). When summarized, the lakes where the fluctuation patterns of K. bostoniensis were similarly grouped differed from the lakes where the fluctuation patterns of potential competitors or predators were similarly grouped.
When examining the fluctuation patterns of environmental factors across the seven lakes where K. bostoniensis appeared (Figure 3), water temperature and DO exhibited similar patterns across lakes, indicating that they were the factors most influenced by seasonality. In contrast, other variables showed differing fluctuation trends depending on the lake. Similar to zooplankton factors, there were no factors in which the fluctuation trends of specific environmental factors grouped similarly across lakes that also grouped similarly for the fluctuation trends of K. bostoniensis.

3.2. Comparison of Differences in Environmental Factors Between Lakes with the Presence and Absence of Kellicottia bostoniensis

The results of the permutation test indicated that the environmental factors showing significant differences between lakes with and without the presence of K. bostoniensis were TOC and TP (Figure 4D,G). For TOC, lakes where K. bosotoniensis was present, TOC reached an average value of 3.94 mg/L, while in the ones where K. bosotoniensis was not present, it had an average value of 5.30 mg/L. Similarly, for TP, lakes with K. bostoniensis showed an average value of 0.03 mg/L, compared to 0.11 mg/L in lakes without its presence. This indicates that lakes where K. bostoniensis was absent had statistically significantly higher concentrations of TOC and TP compared to those where it was present (p < 0.01). In the case of R-A-5, a lake where K. longispina was observed among the non-occurrence lakes of K. bostoniensis, differences were observed in TOC, COD, TN, SS, and Chl-a when compared to lakes where K. bostoniensis was present (Figure 4D–F,H,I). TOC showed a significant difference, with 3.94 mg/L in lakes with K. bostoniensis and 2.00 mg/L in the lake with K. longispina (Figure 4D) (p < 0.01). COD also differed from 6.38 mg/L in lakes with K. bostoniensis and 3.15 mg/L in the lake with K. longispina (Figure 4E) (p < 0.001). Similarly, TN was 1.29 mg/L in lakes with K. bostoniensis and 0.83 mg/L in the lake with K. longispina (Figure 4F) (p < 0.05). SS was 5.38 mg/L in lakes with K. bostoniensis and 1.30 mg/L in the lake with K. longispina (Figure 4H) (p < 0.05). Chl-a showed a statistically significant difference between lakes with the presence of K. bostoniensis (13.95 mg/m3) and the lake with the presence of K. longispina (4.80 mg/m3) (Figure 4I) (p < 0.05).
Among the zooplankton abundances, there were no statistically significant differences between the lakes with the presence of K. bostoniensis and those without its presence (p > 0.05) (Figure 5). Additionally, there were no statistically significant differences in zooplankton abundances between the lakes with the presence of K. bostoniensis and those with the presence of K. longispina.

3.3. Comparison of the Responses of Kellicottia bostoniensis, and Its Potential Competitor/Predator Groups to Environmental Factors

In the DCA, COD, which was identified to exhibit multicollinearity, was excluded, and the analysis was conducted using the remaining eight variables. The analysis results showed that 7 out of the 8 environmental factors, excluding pH, had a significant influence on the formation of the first and second axes in the DCA dimensionality reduction (p < 0.05).
The first axis was primarily explained by DO and water temperature, while the second axis was explained by TOC, TN, TP, SS, and Chl-a, indicating that the explanatory contributions of the seven variables were generally distinct across the two axes (Figure 6). For the first axis, DO showed a positive relationship, whereas water temperature exhibited a negative relationship, indicating an inverse relationship between these two variables. On the second axis, TOC, TN, TP, SS, and Chl-a all displayed negative relationships, suggesting that the patterns among these variables were partially similar.
K. bostoniensis was positioned in the second quadrant of the DCA plot (Figure 6). Thus, its abundance showed a negative relationship with DO and a positive relationship with water temperature. On the other hand, K. bostoniensis exhibited negative relationships with TOC, SS, TP, TN, and Chl-a. These variables, which reflect the trophic state of the lakes, were found to negatively influence the abundance of K. bostoniensis.
When only species classified as potential competitors with K. bostoniensis, due to their similar feeding and habitat characteristics, were plotted on the DCA plane, most species exhibited response patterns to environmental factors that differed from those of K. bostoniensis (Figure 6A). Among the 13 potential competitor species observed, the response patterns to the first axis varied by species. However, all species showed similar responses to the second axis, with the exception of K. longispina. For DO and water temperature, each species displayed distinct response patterns. In contrast, the environmental factors representing the trophic state of the lakes (TOC, TN, TP, SS, and Chl-a) were found to exert a uniformly positive influence across species. Meanwhile, K. longispina, like K. bostoniensis, exhibited response patterns to environmental factors that were distinct from those of K. bostoniensis and its other potential competitors.
Cyclopoid copepods (Cyclopoidae) and rotifers of the genus Asplanchna, known predators of K. bostoniensis and other rotifers, were classified as potential predators and were plotted on the DCA plane (Figure 6B). Among the 12 predator species observed, the DCA plane showed that most species were primarily differentiated along the second axis. Excluding certain species, the majority exhibited a negative relationship with the second axis, indicating a positive influence from environmental factors representing the trophic state of the lakes. However, Acantocyclops sp., Cyclops sp., Mesocyclops sp., and Asplanchna sieboldii displayed negative relationships with these variables, differing from other potential predators. Among them, Mesocyclops sp. was positioned on the same plane as K. bostoniensis, showing a negative response to DO and a positive response to water temperature.
To compare the response patterns of K. bostoniensis and its potential competitor group to lake trophic states, TSI(TP) was calculated based on TP concentrations, and their abundances (log scale) were plotted (Figure 7). K. bostoniensis exhibited relatively high abundances in mesotrophic and lower-level eutrophic environments according to TSI(TP) but showed low abundances in oligotrophic and hypereutrophic environments. In contrast, the potential competitor group exhibited very high abundance in hypereutrophic environments. Despite showing substantial variability due to high residuals, the group displayed a pattern proportional to TSI(TP).

4. Discussion

The trophic state of a lake is a key variable determining the occurrence of K. bostoniensis [23,33,34,35,36,37]. In this study, TOC and TP were also analyzed to be closely associated with the presence/absence of K. bostoniensis (Figure 2). Furthermore, the multivariate analysis conducted to examine relationships with environmental factors revealed that the abundance of K. bostoniensis showed inverse correlation with variables representing the trophic state of the lakes (TOC, TN, TP, Chl-a) (Figure 6). However, when visualizing the univariate relationship between TSI(TP) and the abundance of K. bostoniensis, the species showed relatively low abundance or was absent in oligotrophic and hypereutrophic environments, while exhibiting high abundance in meso-eutrophic trophic conditions. Looking at past studies on the occurrence of K. bostoniensis along the TSI gradient, Macêdo, who studied lakes ranging from oligotrophic to eutrophic (TSI: 30–60), reported an inverse relationship between the calculated trophic state index (TSI) values of the lakes and the average abundance of K. bostoniensis [17]. However, Palazzo’s study, conducted in oligotrophic environments with Chl-a concentrations below 2 µg/L, also found an inverse relationship between lake trophic status and the abundance of K. bostoniensis [14]. When summarizing the contrasting findings from previous studies, it suggests that K. bostoniensis and trophic status do not exhibit a linear relationship but rather a bell-shaped relationship, with the highest abundance occurring within a specific range. Our study, which analyzed a broad range of trophic states from oligotrophic to hypereutrophic environments (TSI: 24.1–96.4), provided meaningful insights by identifying the specific range of trophic conditions suitable for the habitation of K. bostoniensis. Considering numerous studies that reported the occurrence of K. bostoniensis within the mesotrophic to eutrophic range [22,23,36,37], the trophic state suitable for its habitation and population growth is determined to be within the meso-eutrophic range based on the results of this analysis. This finding is significant, as it differentiates the response patterns of K. bostoniensis from those of its potential competitors and predators within the community.
The response of potential competitors, categorized based on trophic state and habitat characteristics, to environmental factors revealed a consistent positive relationship with trophic state variables (TOC, TN, TP, Chl-a) along the second axis of the DCA, excluding K. longispina (Figure 6A). Additionally, the total abundance of species within the potential competitor group was found to be proportional to the TSI values in the surveyed lakes. This indicates that the difference in preferred environments between potential competitors and K. bostoniensis is related to the trophic state of the lakes. Many species within the potential competitor group were from the genus Keratella and Brachionus, which are known to dominate in eutrophic lakes [38,39]. Additionally, species with malleate trophi, including most rotifers and species of the potential competitor group, tend to increase in abundance as eutrophication indicators rise [24,40,41]. Although their abundance trends did not exhibit a strong linear relationship (R2 = 0.21), most species showed high abundances in hypereutrophic lakes. This pattern of occurrence based on trophic state contrasts with the distribution of K. bostoniensis, which showed peak abundances in meso-eutrophic environments. These differences suggest that K. bostoniensis may prefer different habitat conditions compared to species within the potential competitor group, despite their functional similarities.
Rotifers exhibit different feeding characteristics depending on the structure of their trophi [42]. K. bostoniensis possesses malleate trophi, and species within this trophi group are known to gather food using their corona and then chew or cut it with their trophi for ingestion [43,44]. Representative genera with malleate trophi Kellicottia, Brachionus, and Keratella are all classified as polyphagous, known to commonly feed on fine detritus/organic aggregates, picoplankton, and nanoplankton [19,20,21]. However, their clearance rates for these food sources varied significantly. Bogdan’s study comparing the clearance rates of Keratella cochlearis and Kellicottia bostoniensis revealed that Keratella cochlearis exhibited higher clearance rates than Kellicottia bostoniensis across all tracer cell types [45]. Similarly, Starkweather’s research showed that K. cochlearis had higher clearance rates than K. bostoniensis for both live algae and heat-killed algae [46]. In another study by Starkweather, which extended the comparison to species within the genera Brachionus, it was found that species from the Brachionus genus also exhibited higher clearance rates than K. bostoniensis [47]. The physiological differences underlying the variation in feeding rates between K. bostoniensis and other species in the competitor group remain unknown. Factors such as differences in swimming speed [43] and the size range of ingestible food particles [20] may influence these differences. However, the specific reasons why K. bostoniensis exhibits uniquely lower feeding rates and efficiency compared to other species have yet to be elucidated. The low clearance rate of K. bostoniensis limits the maximum food concentration it can effectively utilize. Consequently, beyond a certain threshold food concentration (meso-eutrophic levels, according to the findings of this study), the population size of K. bostoniensis is expected to decrease.
In this study, rotifer predators from the genus Asplanchna and the order Cyclopoida were classified as the potential predator group for K. bostoniensis and analyzed (Figure 6B). However, most species within the predator group, similar to the potential competitor group, exhibited higher abundances in high-nutrient environments, demonstrating response patterns to environmental factors distinct from those of K. bostoniensis. Therefore, it is inferred that a tight predator–prey interaction between K. bostoniensis and its predators was not observed. K. bostoniensis has long spines, which it elongates further in environments with high predation pressure, serving as an effective defense mechanism against predators. This indicates that the population of K. bostoniensis is relatively independent of predation pressure compared to other rotifer species that serve as prey for the same predators. The results of this study suggest that there is no connection between the abundance of K. bostoniensis and the abundance of potential predators, and this disconnection can be inferred to result from its effective defensive adaptations. Some previous studies have shown that this independence from predation pressure can confer an advantage in apparent competition (a form of competition mediated by shared predators), resulting in a positive correlation between K. bostoniensis populations and predation pressure [22,23]. However, such a trend was not observed in this study. Meanwhile, an abundance of Asplanchna species and Cyclopoida is generally known to have a positive correlation with the trophic state of lakes [48,49,50]. This phenomenon occurs because, in addition to the general rotifers, including potential competitors, appearing in large numbers in high-nutrient lakes, alternative food sources such as phytoplankton are also present at high densities. In other words, the increase in prey density leads to a corresponding rise in the abundance of predators. Therefore, these predators’ abundance increases or decreases in response to the total biomass of prey, rather than specifically following changes in that of K. bostoniensis.
In conclusion, this study indicates that the population size of K. bostoniensis is more strongly influenced by environmental factors than by the population size of potential competitors or predators. As K. bostoniensis overlaps in functional groups with potential competitors, the species identified as competitors generally exhibited environmental responses characteristic of general rotifer species. However, this non-native species displayed preferences for different habitat conditions. The unique responses of K. bostoniensis to environmental factors have been noted in some previous studies, and our results highlight its distinct reaction to the trophic state of lakes—an indirect indicator of its feeding environment. Whether this difference in habitat preference arises from direct avoidance of competition when K. bostoniensis enters established native communities or from an evolutionary process of niche differentiation remains unclear from our study. To elucidate the mechanisms underlying the specificity of K. bostoniensis, further research is required that considers vertical distribution within habitats, accumulates data across diverse and broad environmental ranges, and integrates physiological and genetic approaches. What is clear, however, is that its distinctive responses reduce interactions with native communities, implying that the direct impact of K. bostoniensis on native communities after establishment may be minimal. At the same time, these traits suggest that the biological filter—a resistance mechanism of native communities—may have little effect in resisting the introduction of K. bostoniensis. Therefore, the findings of this study provide significant evidence for explaining the widespread (or ongoing) distribution of K. bostoniensis across continents.

Author Contributions

Conceptualization, G.-H.H., K.-H.C. and H.-J.O.; methodology, G.-H.H., K.-H.C. and H.-J.O.; validation, Y.C., D.-H.L., K.-H.C. and H.-J.O.; formal analysis, G.-H.H. and H.-J.O.; data curation,; G.-H.H., Y.C., J.-H.K., K.-H.C. and H.-J.O.; writing—original draft preparation, G.-H.H., K.-H.C. and H.-J.O.; writing—review, K.-H.C. and H.-J.O.; writing—editing, Y.C., D.-H.L., J.-H.K., H.-J.O. and K.-H.C.; visualization, G.-H.H. and K.-H.C.; supervision, K.-H.C. and H.-J.O.; project administration, J.-H.K. and K.-H.C.; funding acquisition, J.-H.K. and K.-H.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by a grant from Environmental Basic Research program funded by Committee for Management of Geum River Basin.

Data Availability Statement

The original contributions presented in this study are included in the article. Further inquiries can be directed to the corresponding author.

Conflicts of Interest

Author Jeong-Hui Kim was employed by the company EcoResearch Inc. The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Appendix A

Table A1. Physical properties of the studied lake.
Table A1. Physical properties of the studied lake.
Lake
No.
LatitudeLongitudeBasin AreaTotal Water StorageEffective StorageHeight of DikeReservoir LENGTH
(ha)(×103 m3)(×103 m3)(m)(m)
R-A-136.0656127.658351078475315.4187
R-A-236.2845127.307213752846282120.0187
R-A-336.3009128.012822005773567923.7276
R-A-436.6027127.231216205284486817.0226
R-A-536.4906127.843542268236816349.0170
R-A-636.8669127.3996847926,61826,37228.7401
R-A-736.8980127.571270614,32913,90733.0237
R-A-836.8685127.599836558791869027.2205
R-A-936.7408127.406033107590639026.7202
R-A-1036.2149128.020810002206220521.7123
R-A-1135.9930126.84801442542415.5200
R-A-1235.9504127.09111043423425.0150

Appendix B

Table A2. Environmental factors for the study lakes (mean ± standard deviation), (minimum–maximum).
Table A2. Environmental factors for the study lakes (mean ± standard deviation), (minimum–maximum).
Lake No.WTpHDOTOCCODTNTPSSChl-a
°C mg/Lmg/Lmg/Lmg/Lmg/m3mg/Lmg/m3
R-A-118.0 ± 5.07.9 ± 0.35.2 ± 4.45.3 ± 0.59.0 ± 1.71.5 ± 0.526.0 ± 8.44.3 ± 2.913.2 ± 10.4
11.4–23.57.5–8.11.2–10.94.9–5.97.6–11.21.0–2.219.0–36.00.9–7.01.9–27.0
R-A-217.4 ± 4.77.7 ± 0.58.9 ± 4.24.4 ± 1.07.5 ± 1.91.2 ± 0.252.5 ± 40.58.9 ± 6.120.4 ± 15.8
11.6–23.16.9–8.24.2–14.43.2–5.64.8–9.00.9–1.312.0–99.01.3–14.93.3–41.4
R-A-315.1 ± 6.17.2 ± 0.27.2 ± 3.33.3 ± 0.35.6 ± 1.20.9 ± 0.118.3 ± 6.92.4 ± 0.811.3 ± 8.5
8.1–22.36.9–7.44.4–11.42.9–3.64.6–7.20.8–1.011.0–27.01.7–3.43.4–22.6
R-A-420 ± 5.97.6 ± 1.09.7 ± 3.33.4 ± 0.46.4 ± 1.21.6 ± 0.446.5 ± 43.59.7 ± 8.513.4 ± 13.0
12.0–25.96.9–9.14.9–12.13.0–4.04.8–7.61.4–2.212.0–110.04.0–22.34.2–32.6
R-A-518.1 ± 8.17.3 ± 0.28.9 ± 2.82.0 ± 0.23.2 ± 0.30.8 ± 0.510.3 ± 7.11.3 ± 1.24.8 ± 2.4
7.3–25.07.0–7.56.9–13.01.7–2.32.8–3.60.3–1.44.0–20.00.5–3.12.3–8.1
R-A-619.0 ± 8.58.0 ± 0.99.4 ± 4.13.1 ± 0.55.6 ± 0.91.2 ± 0.320.0 ± 7.13.2 ± 0.814.1 ± 3.3
7.0–25.87.4–9.44.3–13.32.6–3.64.6–6.71.0–1.612.0–29.02.4–4.411.3–18.8
R-A-715.5 ± 8.28.0 ± 0.710.3 ± 3.23.7 ± 0.36.0 ± 0.71.1 ± 0.320.0 ± 3.63.2 ± 0.714.6 ± 1.7
5.9–24.27.5–9.16.3–13.23.4–4.15.4–7.00.6–1.317.0–25.02.7–4.213.1–17.1
R-A-815.6 ± 8.07.9 ± 0.58.9 ± 2.43.1 ± 0.45.1 ± 1.11.2 ± 0.319.0 ± 11.03.2 ± 2.78.6 ± 5.4
6.6–25.87.2–8.46.6–11.92.7–3.44.0–6.60.9–1.69.0–29.00.7–6.23.5–15.7
R-A-917.7 ± 6.78.4 ± 0.67.5 ± 4.74.6 ± 0.78.7 ± 2.41.8 ± 1.026.0 ± 21.37.3 ± 5.315.3 ± 10.5
9.7–25.28.0–9.20.7–11.33.7–5.25.8–10.60.7–3.09.0–56.01.0–12.45.5–28.4
R-A-1018.7 ± 8.07.8 ± 0.97.9 ± 5.25.3 ± 0.67.5 ± 0.61.1 ± 0.434.3 ± 13.73.8 ± 3.610.8 ± 2.0
9.2–26.77.0–9.01.1–134.6–5.87.0–8.20.8–1.720.0–47.01.0–9.18.5–12.9
R-A-1120.8 ± 5.77.6 ± 0.37.3 ± 1.66.2 ± 2.410.5 ± 4.41.7 ± 0.9122.0 ± 95.118.9 ± 12.220.3 ± 18.3
14.5–25.77.3–7.96.2–9.74.5–9.76.4–16.50.7–3.022.0–243.02.3–31.34.3–45.5
R-A-1220.7 ± 5.07.5 ± 0.37.6 ± 2.07.8 ± 1.513.2 ± 3.22.6 ± 1.5357.5 ± 209.610.5 ± 6.565.6 ± 51.3
15.1–26.07.1–7.85.4–9.86.5–9.910.0–16.10.8–4.5100.0–601.04.4–19.729.2–141

Appendix C

Table A3. Abundance of zooplankton in the study lakes categorized by taxonomic and functional groups (mean ± standard deviation), (minimum–maximum).
Table A3. Abundance of zooplankton in the study lakes categorized by taxonomic and functional groups (mean ± standard deviation), (minimum–maximum).
Lake No.K. bostoniensisK. longispinaCompetitorPredatorRotiferCladoceranCopepod
ind./Lind./Lind./Lind./Lind./Lind./Lind./L
R-A-10.0 ± 0.00.0 ± 0.09.1 ± 18.1117.1 ± 162.2158.9 ± 201.111.2 ± 18.18.3 ± 12.6
0.0–0.00.0–0.00.0–36.30.1–344.80.3–423.00.0–38.20.3–27.1
R-A-20.6 ± 1.10.0 ± 0.092.8 ± 83.711.0 ± 6.6183.5 ± 148.515.1 ± 11.423.4 ± 11.4
0.0–2.30.0–0.020.0–172.23.4–19.253.0–313.91.8–28.010.8–36.6
R-A-31.7 ± 3.30.0 ± 0.02.8 ± 4.87.2 ± 7.326.1 ± 26.48.5 ± 8.635.9 ± 36.2
0.0–6.60.0–0.00.0–10.02.9–18.23.3–55.62.0–21.215.7–90.1
R-A-433.8 ± 66.30.0 ± 0.010.8 ± 11.710.8 ± 10.3147.4 ± 13917.5 ± 9.957.7 ± 25.4
0.0–133.20.0–0.00.8–26.00.0–20.419.2–309.25.2–27.334.0–89.0
R-A-50.0 ± 0.07.7 ± 13.714.7 ± 12.82.3 ± 3.5150.6 ± 221.54.1 ± 4.232.1 ± 33.2
0.0–0.00.0–28.22.8–31.00.1–7.535.4–482.90.0–9.94.7–80.1
R-A-63.3 ± 5.30.0 ± 0.017.1 ± 31.28.8 ± 9.169.1 ± 74.415.0 ± 15.326.9 ± 31.9
0.0–11.20.0–0.00.4–63.82.0–21.613.3–178.40.1–35.91.6–73.4
R-A-71.2 ± 2.40.0 ± 0.09.8 ± 6.57.6 ± 8.276.6 ± 59.59.5 ± 9.130.8 ± 47.5
0.0–4.90.0–0.01.5–16.80.1–17.07.6–137.60.6–21.22.2–101.4
R-A-80.1 ± 0.30.0 ± 0.015.0 ± 20.112.1 ± 11.1192.6 ± 273.756.3 ± 97.154.8 ± 47.0
0.0–0.50.0–0.02.5–44.91.5–23.221.8–596.30.5–201.614.4–105.6
R-A-91.6 ± 2.30.0 ± 0.06.3 ± 10.71.0 ± 0.851.1 ± 88.45.4 ± 4.911.8 ± 11.9
0.0–4.90.0–0.00.6–22.40.0–2.01.4–183.70.6–12.20.0.–25.1
R-A-100.0 ± 0.00.0 ± 0.06.6 ± 10.09.6 ± 5.464.0 ± 73.88.7 ± 5.747.8 ± 31.7
0.0–0.00.0–0.00.6–21.44.5–17.29.4–173.12.9–14.014.7–75.3
R-A-110.0 ± 0.00.0 ± 0.0284.8 ± 191.449.0 ± 47.9611.2 ± 423.719.4 ± 25.8108.4 ± 73.6
0.0–0.00.0–0.010.6–434.20.3–113.737.8–943.50.7–57.319.5–197.4
R-A-120.0 ± 0.00.0 ± 0.0377.6 ± 6787.0 ± 9.1746.0 ± 1231.41.8 ± 2.860.7 ± 42.2
0.0–0.00.0–0.03.0–1393.60.4–19.87.6–2583.00.0–6.122.2–105.7

References

  1. Linders, T.E.W.; Schaffner, U.; Eschen, R.; Abebe, A.; Choge, S.K.; Nigatu, L.; Mbaabu, P.R.; Shiferaw, H.; Allan, E. Direct and indirect effects of invasive species: Biodiversity loss is a major mechanism by which an invasive tree affects ecosystem functioning. J. Ecol. 2019, 107, 2660–2672. [Google Scholar] [CrossRef]
  2. Winter, M.; Schweiger, O.; Klotz, S.; Nentwig, W.; Andriopoulos, P.; Arianoutsou, M.; Basnou, C.; Delipetrou, P.; Didžiulis, V.; Hejda, M.; et al. Plant extinctions and introductions lead to phylogenetic and taxonomic homogenization of the European flora. Proc. Natl. Acad. Sci. USA 2009, 106, 21721–21725. [Google Scholar] [CrossRef]
  3. DeRivera, C.E.; Ruiz, G.M.; Hines, A.H.; Jivoff, P. Biotic resistance to invasion: Native predator limits abundance and distribution of an introduced crab. Ecology 2005, 86, 3364–3376. [Google Scholar] [CrossRef]
  4. Padilla, D.; Williams, S. Beyond ballast water: Aquarium and ornamental trades as sources of invasive species in aquatic ecosystems. Front. Ecol. Environ. 2004; 2, 131–138. [Google Scholar]
  5. Jeschke, J.M.; Strayer, D.L. Usefulness of bioclimatic models for studying climate change and invasive species. Ann. N. Y. Acad. Sci. Ann. N. Y. Acad. Sci. 2008, 1134, 1–24. [Google Scholar] [CrossRef]
  6. Carlin, B. Die Planktonrotatorien des Motal: Taxonomie und Ökologie der Planktonrotatoren. Meddel. Lunds Univ. Limnol. Inst. 1943, 5, 1–260. [Google Scholar]
  7. Eloranta, P. Kellicottia bostoniensis (Rousselet), a planktonic rotifer species new to Finland. Ann. Zool. Fennici. 1988, 25, 249–252. [Google Scholar]
  8. Lopes, R.; Lansac-Toha, F.A.; Vale, R.; Serafim, M., Jr. Comunidade zooplanctônica do reservatório de Segredo. In Reservatório de Segredo: Bases Ecológicas para o Manejo; Eduem: Maringá; Maringa (Brazil) EDUEM: Maringa, Brazil, 1997; pp. 39–60. [Google Scholar]
  9. de Fátima Bomfim, F.; Mantovano, T.; Schwind, L.T.F.; Palazzo, F.; Bonecker, C.C.; Lansac-Tôha, F.A. Geographical spread of the invasive species Kellicottia longispina (Kellicott, 1879) and K. bostoniensis (Rousselet, 1908): A scientometric approach. Acta Sci. Biol. Sci. 2016, 38, 29–36. [Google Scholar] [CrossRef]
  10. Mantovano, T.; Diniz, L.P.; de Oliveira da Conceição, E.; Rosa, J.; Bonecker, C.C.; Bailly, D.; Ferreira, J.H.D.; Rangel, T.F.; Lansac-Toha, F.A. Ecological niche models predict the potential distribution of the exotic rotifer Kellicottia bostoniensis (Rousselet, 1908) across the globe. Hydrobiologia 2021, 848, 299–309. [Google Scholar] [CrossRef]
  11. Sudzuki, M.; Kawakita, M. A rare rotifer, Kellicottia bostoniensis from Japan. Zool. Sci. 1999, 16, 39. [Google Scholar]
  12. Yang, H.M.; Min, G.S. New record of Kellicottia bostoniensis and redescription of two freshwater rotifers from Korea (Rotifera: Monogononta). Anim. Syst. Evol. Divers. 2020, 36, 222–227. [Google Scholar]
  13. Hwang, S.J.; Kwun, S.K.; Yoon, C.G. Water quality and limnology of Korean reservoirs. Paddy Water Environ. 2003, 1, 43–52. [Google Scholar] [CrossRef]
  14. Palazzo, F.; Moi, D.A.; Mantovano, T.; Lansac-Tôha, F.A.; Bonecker, C.C. Assessment of the occurrence and abundance of an exotic zooplankton species (Kellicottia bostoniensis) across a neotropical wetland over 12 years. Limnology 2023, 24, 137–149. [Google Scholar] [CrossRef]
  15. Branco, C.W.C.; Santos-Cabral, L.C.; Kozlowsky-Suzuki, B.; Lopes, V.G.; Puga, A.L.; Macêdo, R.L. Persistence of the non-native Kellicottia bostoniensis (Rousselet, 1908) in a large tropical reservoir. Hydrobiologia 2024, 851, 3039–3060. [Google Scholar] [CrossRef]
  16. Zhdanova, S.M.; Lazareva, V.I.; Bayanov, N.G.; Lobunicheva, E.V.; Rodionova, N.V.; Shurganova, G.V.; Zolotareva, T.V.; Il’in, M.Y. Morphological variability of Kellicottia bostoniensis (Rousselet, 1908) (Rotifera: Brachionidae) in waterbodies of European Russia. Inland Water Biology. 2019, 12, 140–149. [Google Scholar] [CrossRef]
  17. Macêdo, R.L.; Franco, A.C.S.; Klippel, G.; Oliveira, E.F.; Silva, L.H.S.; dos Santos, L.N.; Branco, C.W. Small in size but rather pervasive: The spread of the North American rotifer Kellicottia bostoniensis (Rousselet, 1908) through Neotropical basins. BioInvasions Rec. 2020, 9, 287–302. [Google Scholar] [CrossRef]
  18. Gilbert, J.J. Observations on the susceptibility of some protists and rotifers to predation by Asplanchna girodi. In Rotatoria: Proceedings of the 2nd International Rotifer Symposium Held at Gent, September 17–21; Springer: Dordrecht, The Netherlands, 1979; pp. 87–91. [Google Scholar]
  19. Gilbert, J.J. Food niches of planktonic rotifers: Diversification and implications. Limnol. Oceanogr. 2022, 67, 2218–2251. [Google Scholar] [CrossRef]
  20. Sanders, R.W.; Porter, K.G.; Bennett, S.J.; DeBiase, A.E. Seasonal patterns of bacterivory by flagellates, ciliates, rotifers, and cladocerans in a freshwater planktonic community. Limnol. Oceanogr. 1989, 34, 673–687. [Google Scholar] [CrossRef]
  21. Erman, L.A. On the utilization of the reservoirs’ trophic resources by plankton rotifers. Byull. Mosk. Obshch. Ispyt. Prir. 1962, 67, 32–47. [Google Scholar]
  22. Karpowicz, M.; Zieliński, P.; Grabowska, M.; Ejsmont-Karabin, J.; Kozłowska, J.; Feniova, I. Effect of eutrophication and humification on nutrient cycles and transfer efficiency of matter in freshwater food webs. Hydrobiologia 2020, 847, 2521–2540. [Google Scholar] [CrossRef]
  23. Shurganova, G.V.; Zolotareva, T.V.; Kudrin, I.A.; Zhikharev, V.S.; Gavrilko, D.E.; Il’in, M.I. Abundance of related species, Kellicottia bostoniensis (Rousselet, 1908) and K. longispina (Kellicott, 1879) (Rotifera: Brachionidae), in the zooplankton communities of the Pustynskaya lake-river system (Nizhny Novgorod region). Russ. J. Biol. Invasions 2021, 12, 219–231. [Google Scholar] [CrossRef]
  24. Oh, H.J.; Jeong, H.G.; Nam, G.S.; Oda, Y.; Dai, W.; Lee, E.H.; Hwang, S.J.; Chang, K.H. Comparison of taxon-based and trophi-based response patterns of rotifer community to water quality: Applicability of the rotifer functional group as an indicator of water quality. Anim. Cells Syst. 2017, 21, 133–140. [Google Scholar] [CrossRef] [PubMed]
  25. Fontaneto, D.; De Smet, W.H. Rotifera. In Gastrotricha and Gnathifera; Schmidt-Rhaesa, A., Ed.; De Gruyter: Berlin, Germany, 2015; pp. 217–300. [Google Scholar]
  26. Sharma, S.S.S. Zooplankton diversity: Freshwater planktonic and semi-planktonic Rotifera. In Water Quality Assessment Biomonitoring and Zooplankton Diversity; North Eastern Hill University: Meghalaya, India, 2001; pp. 190–215. [Google Scholar]
  27. Karabin, A. The pressure of pelagic predators of the genus Mesocyclops (Copepoda, Crustacea) on small zooplankton. Ekol. Pol. 1978, 26, 241–257. [Google Scholar]
  28. Stemberger, R.S.; Evans, M.S. Rotifer seasonal succession and copepod predation in Lake Michigan. J. Great Lakes Res. 1984, 10, 417–428. [Google Scholar] [CrossRef]
  29. Gilbert, J.J.; Williamson, C.E. Predator-prey behavior and its effect on rotifer survival in associations of Mesocyclops edax, Asplanchna girodi, Polyarthra vulgaris, and Keratella cochlearis. Oecologia 1978, 37, 13–22. [Google Scholar] [CrossRef] [PubMed]
  30. Hothorn, T.; Hornik, K.; van de Wiel, M.A.; Zeileis, A. A Lego system for conditional inference. Am. Stat. 2006, 60, 257–263. [Google Scholar] [CrossRef]
  31. Oksanen, J.; Simpson, G.; Blanchet, F.; Kindt, R.; Legendre, P.; Minchin, P.; O’Hara, R.; Solymos, P.; Stevens, M.; Szoecs, E.; et al. vegan: Community Ecology Package; R package version 2.6-4. Available online: https://CRAN.R-project.org/package=vegan (accessed on 10 December 2024).
  32. Hill, M.O.; Gauch, H.G. Detrended correspondence analysis: An improved ordination technique. Vegetatio 1980, 42, 47–58. [Google Scholar] [CrossRef]
  33. Campbell, R.S. Vertical distribution of the plankton Rotifera in Douglas Lake, Michigan, with special reference to submerged depression individuality. Ecol. Monogr. 1941, 11, 2–19. [Google Scholar] [CrossRef]
  34. Arnemo, R.; Berzins, B.; Grönberg, B.; Mellgren, I. The dispersal in Swedish waters of Kellicottia bostoniensis (Rousselet) (Rotatoria). Oikos 1968, 19, 351–358. [Google Scholar] [CrossRef]
  35. Zhdanova, S.M.; Dobrynin, A.E. Kellicotia bostoniensis (Rousselet, 1908) (Rotifera: Brachionidae) in waterbodies of European Russia. Inland Water Biol. 2011, 4, 39–46. [Google Scholar] [CrossRef]
  36. Zolotareva, T.V.; Shurganova, G.V.; Kudrin, I.A.; Gavrilko, D.E.; Zhikharev, V.S.; Kolesnikov, A.A.; Il’in, M.Y. Morphological characteristics of Kellicottia bostoniensis (Rousselet, 1908) (Rotifera: Brachionidae) in water bodies of the Middle Volga basin. Povolzhskiy J. Ecol. 2021, 1, 16–34. [Google Scholar] [CrossRef]
  37. Podshivalina, V.N.; Semenova, A.S. Kellicottia bostoniensis (Rousselet, 1908) and K. longispina (Kellicott, 1879) (Rotifera: Brachionidae): Features of occurrence and distribution in the lakes of the Upper and Middle Volga Region. Russ. J. Biol. Invasions 2023, 14, 376–388. [Google Scholar] [CrossRef]
  38. García-Chicote, J.; Armengol, X.; Rojo, C. Zooplankton species as indicators of trophic state in reservoirs from Mediterranean river basins. Inland Waters 2019, 9, 113–123. [Google Scholar] [CrossRef]
  39. Ejsmont-Karabin, J. The usefulness of zooplankton as lake ecosystem indicators: Rotifer trophic state index. Pol. J. Ecol. 2012, 60, 339–350. [Google Scholar]
  40. Oh, H.J.; Chang, K.H.; Seo, D.I.; Nam, G.S.; Lee, E.H.; Jeong, H.G.; Yoon, J.D.; Oh, J.M. Zooplankton community as an indicator for environmental assessment of aquatic ecosystem: Application of rotifer functional groups for evaluating water quality in eutrophic reservoirs. J. Environ. Impact Assess. 2017, 26, 404–417. [Google Scholar]
  41. Qian, F.P.; Wen, X.L.; Xi, Y.L. Seasonal variations of rotifer communities in three climatic zones: Effects of environmental parameters on changes of functional groups. Limnology 2022, 23, 1–16. [Google Scholar] [CrossRef]
  42. Obertegger, U.; Smith, H.A.; Flaim, G.; Wallace, R.L. Using the guild ratio to characterize pelagic rotifer communities. Hydrobiologia 2011, 662, 157–162. [Google Scholar] [CrossRef]
  43. Koste, W.; Shiel, R.J. Rotifera from Australian inland waters. II. Epiphanidae and Brachionidae (Rotifera: Monogononta). Invert. Taxon. 1987, 7, 949–1021. [Google Scholar] [CrossRef]
  44. Lokko, K.; Virro, T.; Kotta, J. Seasonal variability in the structure and functional diversity of psammic rotifer communities: Role of environmental parameters. Hydrobiologia 2017, 796, 287–307. [Google Scholar] [CrossRef]
  45. Bogdan, K.G.; Gilbert, J.J.; Starkweather, P.L. In situ clearance rates of planktonic rotifers. In Rotatoria; Pejler, B., Ed.; Springer: Dordrecht, The Netherlands, 1980; pp. 73–77. [Google Scholar]
  46. Starkweather, P.L.; Bogdan, K.G. Detrital feeding in natural zooplankton communities: Discrimination between live and dead algal foods. Hydrobiologia 1980, 73, 83–85. [Google Scholar] [CrossRef]
  47. Starkweather, P.L. Aspects of the feeding behavior and trophic ecology of suspension-feeding rotifers. In Rotatoria; Pejler, B., Ed.; Springer: Dordrecht, The Netherlands, 1980; pp. 63–72. [Google Scholar]
  48. Stemberger, R.S.; Gilbert, J.J. Body size, food concentration, and population growth in planktonic rotifers. Ecology 1985, 66, 1151–1159. [Google Scholar] [CrossRef]
  49. Perbiche-Neves, G.; Saito, V.S.; Previattelli, D.; Da Rocha, C.E.; Nogueira, M.G. Cyclopoid copepods as bioindicators of eutrophication in reservoirs: Do patterns hold for large spatial extents? Ecol. Indic. 2016, 70, 340–347. [Google Scholar] [CrossRef]
  50. Chang, K.H.; Nishibe, Y.; Nakano, S.I. Feeding habits of omnivorous Asplanchna: Comparison of diet composition among Asplanchna herricki, A. priodonta, and A. girodi in pond ecosystems. J. Limnol. 2010, 69, 209–216. [Google Scholar] [CrossRef]
Figure 1. Geographic distribution of the lakes in this study. The green diamonds represent the lakes where Kellicottia bostoniensis appeared in our study, while the red diamonds indicate lakes where it did not appear.
Figure 1. Geographic distribution of the lakes in this study. The green diamonds represent the lakes where Kellicottia bostoniensis appeared in our study, while the red diamonds indicate lakes where it did not appear.
Water 17 00273 g001
Figure 2. The fluctuation trends of zooplankton factors (abundance of potential competitors (Competitors, b); abundance of potential predator (Predator, c); abundance of total rotifers (Rotifers, d); Abundance of total cladocerans (Cladocerans, e); abundance of total copepods (Copepods, f)), including Kellicottia bostoniensis, in the seven lakes where K. bostoniensis appeared are represented, with variables arranged in rows and lakes in columns (a). Abundance has been converted to a log scale. The numbers 17 in the subfigures refer to lakes R-A-2, R-A-3, R-A-4, R-A-6, R-A-7, R-A-8, and R-A-9, in that order.
Figure 2. The fluctuation trends of zooplankton factors (abundance of potential competitors (Competitors, b); abundance of potential predator (Predator, c); abundance of total rotifers (Rotifers, d); Abundance of total cladocerans (Cladocerans, e); abundance of total copepods (Copepods, f)), including Kellicottia bostoniensis, in the seven lakes where K. bostoniensis appeared are represented, with variables arranged in rows and lakes in columns (a). Abundance has been converted to a log scale. The numbers 17 in the subfigures refer to lakes R-A-2, R-A-3, R-A-4, R-A-6, R-A-7, R-A-8, and R-A-9, in that order.
Water 17 00273 g002
Figure 3. The seasonal dynamics of environmental factors in the seven lakes where Kellicottia bostoniensis appeared are represented, with variables arranged in rows and lakes in columns. From the top: water temperature (WT, °C; a), hydrogen ion concentration (pH; b), dissolved oxygen concentration (DO, mg/L; c), total organic carbon concentration (TOC, mg/L; d), chemical oxygen demand (COD, mg/L; e), total nitrogen concentration (TN, mg/L; f), total phosphorus concentration (TP, mg/L; g), suspended solids concentration (SS, mg/L; h), and chlorophyll a concentration (mg/mL; i). The numbers 17 in the subfigures refer to lakes R-A-2, R-A-3, R-A-4, R-A-6, R-A-7, R-A-8, and R-A-9, in that order.
Figure 3. The seasonal dynamics of environmental factors in the seven lakes where Kellicottia bostoniensis appeared are represented, with variables arranged in rows and lakes in columns. From the top: water temperature (WT, °C; a), hydrogen ion concentration (pH; b), dissolved oxygen concentration (DO, mg/L; c), total organic carbon concentration (TOC, mg/L; d), chemical oxygen demand (COD, mg/L; e), total nitrogen concentration (TN, mg/L; f), total phosphorus concentration (TP, mg/L; g), suspended solids concentration (SS, mg/L; h), and chlorophyll a concentration (mg/mL; i). The numbers 17 in the subfigures refer to lakes R-A-2, R-A-3, R-A-4, R-A-6, R-A-7, R-A-8, and R-A-9, in that order.
Water 17 00273 g003
Figure 4. Results of the permutation test for differences in environmental factors (water temperature (WT, °C; A), hydrogen ion concentration (pH; B), dissolved oxygen concentration (DO, mg/L; C), total organic carbon concentration (TOC, mg/L; D), chemical oxygen demand (COD, mg/L; E), total nitrogen concentration (TN, mg/L; F), total phosphorus concentration (TP, mg/L; G), suspended solids concentration (SS, mg/L; H), and chlorophyll a concentration, mg/mL, I) among the lakes with the presence of K. bostoniensis, the lakes without its presence, and the lakes with the presence of K. longispina (Permutation test: p < 0.05 *, p < 0.01 **, p < 0.001 ***).
Figure 4. Results of the permutation test for differences in environmental factors (water temperature (WT, °C; A), hydrogen ion concentration (pH; B), dissolved oxygen concentration (DO, mg/L; C), total organic carbon concentration (TOC, mg/L; D), chemical oxygen demand (COD, mg/L; E), total nitrogen concentration (TN, mg/L; F), total phosphorus concentration (TP, mg/L; G), suspended solids concentration (SS, mg/L; H), and chlorophyll a concentration, mg/mL, I) among the lakes with the presence of K. bostoniensis, the lakes without its presence, and the lakes with the presence of K. longispina (Permutation test: p < 0.05 *, p < 0.01 **, p < 0.001 ***).
Water 17 00273 g004
Figure 5. Results of the permutation test for differences in zooplankton factors (abundance of potential competitors (Competitors; D); abundance of potential predator (Predator; E); abundance of Total Rotifers (Rotifers; A); abundance of total cladocerans (Cladocerans; B); abundance of total copepods (Copepods; C)) among the lakes with the presence of K. bostoniensis, the lakes without its presence, and the lakes with the presence of K. longispina (permutation test: p < 0.05 *, p < 0.01 **, p < 0.001 ***).
Figure 5. Results of the permutation test for differences in zooplankton factors (abundance of potential competitors (Competitors; D); abundance of potential predator (Predator; E); abundance of Total Rotifers (Rotifers; A); abundance of total cladocerans (Cladocerans; B); abundance of total copepods (Copepods; C)) among the lakes with the presence of K. bostoniensis, the lakes without its presence, and the lakes with the presence of K. longispina (permutation test: p < 0.05 *, p < 0.01 **, p < 0.001 ***).
Water 17 00273 g005
Figure 6. Detrended Correspondence Analysis (DCA) of K. bostoniensis abundance, zooplankton community, and their responses to environmental factors. (A) Graph showing only the potential competitor group. (B) Graph showing only the potential predator group. Environmental factors with p < 0.05 are represented. Red circle indicate K. bostoniensis and, black circles indicate the target species that assumed to potentially have a biological interaction with K. bostoniensis ((A) potential competitor, (B) potential predator). Grey circles indicate remaining observed species, and plus signs indicate sampling sites.
Figure 6. Detrended Correspondence Analysis (DCA) of K. bostoniensis abundance, zooplankton community, and their responses to environmental factors. (A) Graph showing only the potential competitor group. (B) Graph showing only the potential predator group. Environmental factors with p < 0.05 are represented. Red circle indicate K. bostoniensis and, black circles indicate the target species that assumed to potentially have a biological interaction with K. bostoniensis ((A) potential competitor, (B) potential predator). Grey circles indicate remaining observed species, and plus signs indicate sampling sites.
Water 17 00273 g006
Figure 7. Responses of K. bostoniensis and the potential competitor group to TSI (TP). The dotted line represents the regression line for the potential competitor group in relation to TSI (TP). Oligotrophic: 20 < TSI < 40; Mesotrophic: 40 < TSI < 50; Eutrophic: 50 < TSI < 60; Hypereutrophic: 60 < TSI.
Figure 7. Responses of K. bostoniensis and the potential competitor group to TSI (TP). The dotted line represents the regression line for the potential competitor group in relation to TSI (TP). Oligotrophic: 20 < TSI < 40; Mesotrophic: 40 < TSI < 50; Eutrophic: 50 < TSI < 60; Hypereutrophic: 60 < TSI.
Water 17 00273 g007
Table 1. Classification of functional groups of rotifers that appeared in the studied lake (at genus level).
Table 1. Classification of functional groups of rotifers that appeared in the studied lake (at genus level).
Functional GroupHabitat Type
PlanktonicLittoral, Periphytic
trophiMalleate typeAnuraeopsis, Brachionus, Kellicottia, KeratellaEuchlanis, Lecane, Lepadella,
Mytilina, Platyias, Trichotria
Virgate typeAscomorpha, Cephalodella, Ploesoma,
Polyarthra, Synchaeta, Trichocerca
-
Incudate typeAsplachna-
Ramate typeRotariaPhilodina
Malleoramate typeConochilus, Filinia, Hexarthra, PompholyxTestudinella, Tetramastix
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hong, G.-H.; Choi, Y.; Lee, D.-H.; Kim, J.-H.; Chang, K.-H.; Oh, H.-J. The Response Patterns of Non-Native Rotifer Species Kellicottia bostoniensis (Rousselet, 1908) to Environmental Factors and Its Relationship with Potential Competitors and Predators. Water 2025, 17, 273. https://doi.org/10.3390/w17020273

AMA Style

Hong G-H, Choi Y, Lee D-H, Kim J-H, Chang K-H, Oh H-J. The Response Patterns of Non-Native Rotifer Species Kellicottia bostoniensis (Rousselet, 1908) to Environmental Factors and Its Relationship with Potential Competitors and Predators. Water. 2025; 17(2):273. https://doi.org/10.3390/w17020273

Chicago/Turabian Style

Hong, Geun-Hyeok, Yerim Choi, Dae-Hee Lee, Jeong-Hui Kim, Kwang-Hyeon Chang, and Hye-Ji Oh. 2025. "The Response Patterns of Non-Native Rotifer Species Kellicottia bostoniensis (Rousselet, 1908) to Environmental Factors and Its Relationship with Potential Competitors and Predators" Water 17, no. 2: 273. https://doi.org/10.3390/w17020273

APA Style

Hong, G.-H., Choi, Y., Lee, D.-H., Kim, J.-H., Chang, K.-H., & Oh, H.-J. (2025). The Response Patterns of Non-Native Rotifer Species Kellicottia bostoniensis (Rousselet, 1908) to Environmental Factors and Its Relationship with Potential Competitors and Predators. Water, 17(2), 273. https://doi.org/10.3390/w17020273

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop