Next Article in Journal
Silver Doped Zinc Stannate (Ag-ZnSnO3) for the Photocatalytic Degradation of Caffeine under UV Irradiation
Next Article in Special Issue
Additive Effects of Sediment and Nutrient on Leaf Litter Decomposition and Macroinvertebrates in Hyporheic Zone
Previous Article in Journal
Fuzzy or Non-Fuzzy? A Comparison between Fuzzy Logic-Based Vulnerability Mapping and DRASTIC Approach Using a Numerical Model. A Case Study from Qatar
Previous Article in Special Issue
Combined Effects of Summer Water Temperature and Current Velocity on the Distribution of a Cold-Water-Adapted Sculpin (Cottus nozawae)
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Relative Effects of Biotic and Abiotic Factors on the Recruitment of Freshwater Mussels (Margaritifera laevis)

1
Graduate School of Agriculture, Hokkaido University, N9W9 Sapporo, Hokkaido 060-8589, Japan
2
Forest Research Institute, Hokkaido Research Organization, Kosyunai, Bibai, Hokkaido 079-0198, Japan
3
Shiretoko Museum, Honmachi 49, Shari, Shari-gun, Hokkaido 099-4113, Japan
4
Department of Biology, University of North Carolina at Greensboro, Greensboro, NC 27412, USA
5
Aqua Restoration Research Center, Public Works Research Institute, Kawashima Kasada-machi, Kakamigahara, Gifu 501-6021, Japan
*
Author to whom correspondence should be addressed.
Water 2021, 13(9), 1289; https://doi.org/10.3390/w13091289
Submission received: 25 March 2021 / Revised: 27 April 2021 / Accepted: 30 April 2021 / Published: 3 May 2021
(This article belongs to the Special Issue Endangered Freshwater Ecosystems: Threats and Conservation Needs)

Abstract

:
Freshwater mussels, Unionoida, are endangered across the globe due to recruitment failure. In the present study, with general linear mixed models, we investigated the relative effects of biotic (host fish density) and abiotic (water depth, fine sediment, water temperature, and water quality) factors on the recruitment of Margaritifera laevis in 10 streams of Hokkaido, northern Japan. We additionally examined the factors regulating the density of the host fish Oncorhynchus masou masou with general linear models. The proportion of juvenile mussels had a unimodal relationship with the host density, which was the most influential factor among the others examined. The positive relationship between mussel recruitment and host density can be attributed to an increased host fish infection rate. The negative correlation between mussel recruitment and host density at high fish densities may be due to reduced larval growth on host fish that are in poor physical condition. We also found that host fish density was negatively affected by nutrient enrichment. Our results suggest that mitigating water quality degradation to recover host fish density should be prioritized to improve mussel recruitment. Although stock enhancement is effective for increasing the salmon population density, excess stocking can further disturb mussel recruitment.

1. Introduction

Biodiversity has rapidly declined during the sixth mass extinction [1]. Clarifying the causes of species decline is an essential first step for evidence-based conservation. Previous research has focused primarily on population size as an indicator of current population status e.g., [2]. However, population size does not necessarily reflect actual extinction risks. For example, recruitment can be impaired even in large populations of long-lived species [3]. Such apparently healthy populations (i.e., large populations) of long-lived species will go extinct if recruitment does not recover in the future. Therefore, the causes of recruitment failure must be carefully examined to conserve populations of long-lived species [4,5].
In freshwater ecosystems, unionid mussels (order Unionoida) are among the most long-lived species (often reaching 100 years of age), and they play significant roles in purifying the water and providing habitats for other freshwater organisms [6,7]. As with other long-lived species, many species of freshwater mussels have undergone substantial declines across the globe [8,9,10]. For example, approximately 70% of Japanese unionid mussels are listed on the national Red List as either endangered or threatened species [11,12]. While the acute impacts of physical habitat destruction are well-recognized as a major cause of mussel decline, less understood are the chronic stressors that impair the recruitment processes of unionid mussels. Therefore, an in-depth understanding of the factors limiting mussel recruitment is of critical importance for conserving this highly imperilled group of animals.
Biotic and abiotic factors may be involved in chronic recruitment failure in unionid mussels. Freshwater mussels have a complex life cycle in which larvae called glochidia are obligate parasites of the gills or fins of host fish. After the larval parasitic life stage (which typically lasts for several weeks to months), juvenile mussels shed from the host fish and burrow into the substrate [13,14,15]. Therefore, the recruitment of freshwater mussels is influenced by biotic (i.e., host fish) and abiotic factors during their parasitic and benthic life stages. The host fish density in the surrounding environment may regulate the host fish infection rate (see, e.g., [16,17,18,19]). Moreover, excessive sedimentation leads to recruitment failure in freshwater mussels through the impairment of benthic environments [20,21,22]. There is increasing awareness that a combination of biotic and abiotic factors may determine the recruitment of freshwater mussels [16,21,23]. An approach that incorporates both biotic and abiotic factors may illuminate the relative effects of these factors, providing useful insights into the conservation of unionid mussels. However, the relative effect is rarely tested because biotic and abiotic factors have often been considered separately.
Examining the factors that regulate host fish density is also critically important for the conservation of unionid mussels. This is because habitat degradation indirectly influences the parasitic stage, which is crucial to the initial growth of juvenile mussels, by reducing the host fish density [24,25]. A combined investigation of the determinants of mussel recruitment and host fish density can improve our understanding of the conservation of unionid mussels.
Margaritifera laevis is a Japanese freshwater mussel whose glochidia are obligate parasites of the gills of juvenile masu salmon, Oncorhynchus masou masou. This long-lived species has an estimated lifespan of 80 years and is distributed in cold streams in Sakhalin, Russia, and the Honshu and Hokkaido Islands of Japan [26]. The population of M. laevis has experienced nationwide decline over the past few decades, and M. laevis is now listed as an endangered species on the Red List of Japan [12]. Although northern populations of M. laevis in Hokkaido, Japan seem to remain intact, the disproportionate lack of juvenile mussels implies recent recruitment failures [9]. In the present study, we examined the relative effects of biotic (host fish density) and abiotic (sedimentation and water quality) factors on the recruitment of freshwater mussels. In doing so, we statistically evaluated the relative effects of these potential factors on the proportion of juvenile mussels in the population (i.e., a proxy for juvenile recruitment). Furthermore, we explored the abiotic factors limiting the density of O. masou masou hosting M. laevis larvae. We expected that biotic factors (i.e., host fish density) had the most influence on the recruitment of M. laevis because released glochidia from female mussels cannot begin the juvenile stage if they cannot encounter the proper host fish individuals.

2. Materials and Methods

2.1. Quantitative Sampling and Species Detection

The present study was conducted in 10 streams in Hokkaido, Japan (eight streams in northern Hokkaido, one stream in central Hokkaido, and another stream in eastern Hokkaido) (Figure 1). Hokkaido is a large island (83,450 km2) and one of the most famous agricultural prefectures in Japan. The proportion of farmland in the catchment area of each study stream was 0.25 ± 0.24 (measure ± SD) (Table S1). We sampled mussels in the summers of 2016 and 2017. We established a 2–10 km segment within each stream based on accessibility and mussel distribution. Within each segment, we established two or three sampling sites (width of the sites; 2.5–13.0 m) with 30–50 m reach lengths that included riffles and rapids. There were a total of 26 sampling sites. At each site, 15 quadrats were sampled with a Surber net (0.09 m2, 0.5 mm mesh size). Each quadrat was thoroughly investigated by excavating to a depth of 10 cm in the substrate and sieving the substrate immediately through a 2 mm mesh sieve. All collected mussels were photographed with a digital camera in the field, and the shell length and height were measured with image analysis software (ImageJ 1.50i) [27] in the laboratory. Margaritifera togakushiensis, which are taxonomically close to M. laevis, are found in Japan [28]. In our study, we distinguished M. laevis from M. togakushiensis by the method based on the ratio of the shell length to the shell height; this method has been shown to provide 80% accuracy for species identification in northern and central Hokkaido and 92% accuracy for species identification in eastern Hokkaido [29]. With this method, we eliminated the possibility of misidentifying M. laevis as M. togakushiensis to the extent possible.

2.2. Age Estimation

To assess the recruitment status of the study populations, we estimated the proportion of juvenile mussels (the number of juvenile mussels divided by the total number of mussels) at each study site. In this study, mussels under 10 years old, which are considered to be sexually immature, were treated as juvenile mussels [30]. Juvenile mussels were classified based on the following growth models. The age of mussels can be estimated by counting the number of growth rings engraved onto the shell ligament. However, this method requires the sacrifice of live mussels and is therefore inappropriate for use with endangered species [31]. Therefore, we modeled a shell length-age relationship to minimize the number of sacrificed individuals. We developed separate regional shell length–age relationships for northern Hokkaido, central Hokkaido, and eastern Hokkaido given the geographical distance between the study streams (Figure 1). We collected 150 mussels from the study streams (50 mussels in northern Hokkaido, 28 mussels in central Hokkaido, and 72 mussels in eastern Hokkaido) and measured their shell length with image analysis software (ImageJ 1.50i) [27]. We estimated the age of the mussels by counting the number of growth rings under a dissecting microscope. Since the top of the ligament was often damaged, it was often not possible to read the exact age. Therefore, the growth rings of the damaged portion were estimated using the length of the ligaments and the number of growth rings of undamaged individuals [32]. We modeled the shell length–age relationships in the three regions with the following four nonlinear growth models (see, e.g., [31,33,34]):
von   Bertalanffy   growth   curve :   L t = L ( 1 e k ( t t 0 ) )
Gompertz   function :   L t = L e ae kt
Hyperbolic   saturation   function :   L t = L k ( t t 0 ) 1 + k ( t t 0 )
Logistic   function :   L t = L 1 + e k ( t a ) ,
where L t is the shell length (mm) of a mussel at age t, L is the theoretical maximum length (mm), k is a growth coefficient, t 0 is the theoretical age at zero length, and a is a constant. It must be noted that the classical von Bertalanffy growth model used in previous studies with freshwater mussels is modified here, following Urban [34], and is not the “generalized” von Bertalanffy growth model. The performance of each model was assessed on the basis of their sum of squared residuals (RSS) values, and the model with the minimum RSS value for each region was applied.

2.3. Abiotic Factors

We measured five abiotic factors: water depth, proportion of fine sediment, summer maximum water temperature, total nitrogen concentration (TN), and total phosphorous concentration (TP). The water depth and proportion of fine sediment were measured at each study site, whereas the summer maximum water temperature, TN, and TP were measured at a single location within each study segment.
We measured the water depth (cm) at the corners and center of each quadrat. The mean water depth was calculated as the average of all the measurements at each study site. We took three sediment samples (right bank side, left bank side, centerline of the site; 300 mL) from the substrate with a trowel and sieved them with 2 mm and 0.042 mm mesh sieves to identify sediments fine enough for a mussel to burrow in. The sieved sediments that contained <2 mm and >0.042 mm particles were brought back to the laboratory in a cooler box with ice packs and stored in a freezer. In the laboratory, we sieved these sediments with a 0.25 mm mesh sieve (i.e., fine enough to impede interstitial water flow) to classify the sediments into fine sediments (<0.25 mm) that went through the sieve and coarse sediments (>0.25 mm) that remained on the sieve. After drying the sediments at 60 °C for 24 h, we weighed and calculated the proportion of fine sediment (fine sediment/fine + coarse sediments). The mean proportion of fine sediment was calculated as the average for the three samples at each study site. The maximum summer water temperature (in °C) was measured using data loggers (CO-U20L-04, Onset Computer Corporation, U.S.) located in each study segment. We defined the maximum summer water temperature as the maximum water temperature during July and August of 2018 based on hourly records of water temperature. To measure TN and TP (10−2 g l−L), we collected surface water samples from each study segment in autumn 2017 and filtered them immediately through precombusted glass-microfiber filters (Whatman GF/F, GE Healthcare UK Ltd., Buckinghamshire, UK) using syringes. The filtered water samples were quickly brought back to the laboratory in a cooler box with ice packs and stored in a freezer at −15 °C until further analyses. The TN and TP of each sample were analyzed by the colorimetric method using a continuous flow injection analyzer (AACS-4, BL-Tech Co. Ltd., Osaka, Japan) after NaOH-K2S2O8 digestion. The index of nutrient enrichment was calculated by principal component analysis of TN and TP because these values showed a strong correlation with each other (r > 0.7).

2.4. Biotic Factors

We assembled abundance data for the host fish, O. masou masou, across the 10 study streams. Only 1-year-old and younger fish were observed in the streams. In one or two survey reaches (survey reach length; 40.6 ± 23.2 m (measure ± SD)) within 2 km of each mussel survey study segment, one-pass host fish sampling was conducted with a combination of cast-net sampling and electrofishing from 2005 to 2016. We used the cast-net method as an efficient method to collect swimming fishes in riffles, which are difficult to collect with electrofishing. In the analysis, we used data from five or six non-consecutive years during the sampling period that were one or two years apart (Table S2). We calculated the mean host fish density during the sampling period and used that value in the following statistical analyses. The data were provided by the Hokkaido Regional Development Bureau and the Salmon and Freshwater Fisheries Research Institute.

2.5. Statistical Analysis

First, we developed general linear mixed models (GLMMs; n = 26) with a binomial error distribution via the R package “lme4” [35] to examine the relative effects of biotic and abiotic factors on the recruitment success of M. laevis. We used the proportion of juvenile mussels as a response variable and the biotic (logarithmically transformed host fish density) and abiotic factors (water depth, proportion of fine sediment, maximum summer water temperature) as explanatory variables. Since the host fish density and the index of nutrient enrichment showed a strong correlation (r < −0.7), we did not add these variables to the model simultaneously to account for multicollinearity. Thus, we used the host fish density as the explanatory variable after a preliminary analysis that compared the Akaike information criterion (AIC) values of the four GLMMs developed with each variable (Table S3). All explanatory variables were normalized (i.e., scaled to a mean of 0 and a standard deviation of 1). Stream identity was included as a random effect to account for the effect of unmodeled factors attributable to differences among the streams. To identify plausible models, we selected the best model with the lowest AIC via the R package “MuMIn” [36]. We performed model averaging when multiple best models (ΔAIC < 2) were found. We interpreted the significance of the parameters in the best model or the averaged model on the basis of a 95% confidence interval that did not include 0. We evaluated the relative effect of each factor based on estimates of the standardized partial regression coefficients of the explanatory variables. Because the standard partial regression coefficient does not depend on the unit of each variable, it was possible to compare the degree to which the explanatory variable affected the objective variable.
Second, we used general linear models (GLMs; n = 10) with a Gaussian error distribution to examine the abiotic factors affecting the host fish density. We used the host fish density as a response variable and the four abiotic factors (water depth, proportion of fine sediment, summer maximum water temperature and index of nutrient enrichment) as explanatory variables. The same procedures used in the analyses of the recruitment success of M. laevis were adopted for the model selection. The model selection was based on the AIC for a small sample situation (AICc). All analyses were performed in R software v. 3.6.3 [37].

3. Results

The best-fitting models for estimating the shell length–age relationship were the Gompertz function in northern and central Hokkaido and the hyperbolic saturation function in eastern Hokkaido (Figure S1, Table S4). The estimated shell length at 10 years of age was 49.76 mm (in northern Hokkaido), 17.00 mm (in central Hokkaido), and 28.56 mm (in eastern Hokkaido).
The number of individuals and the proportion of juvenile mussels varied among streams (Figure 2). Among the 2279 mussels collected (shell length: 8.5–148.7 mm), 388 mussels were estimated as juveniles based on the shell length-age relationships. The proportion of juvenile mussels was 0.15 ± 0.17 (measure ± SD), and no juveniles were found at five sites in the three streams. The mean water depth and the proportion of fine sediment of the studied sites were 29.8 ± 10.2 cm and 0.18 ± 0.11, respectively, and the maximum summer water temperature, TN, TP, and host fish density of the studied segment were 21.5 ± 2.79 °C, 0.41 ± 0.34 10−2 g L−1, 0.02 ± 0.01 10−2 g L−1, and 0.24 ± 0.19 ind. m−2, respectively (Table S1). The first principal component explained variations in the observed TN and TP values, which had a strong positive correlation with both variables (r > 0.9).
Host fish density was the most influential factor of all the potential influencing ones. The model selection process revealed that the proportion of juvenile mussels was affected only by the host fish density (Table 1, Marginal R2: 0.246, Conditional R2: 0.250). The host fish density had a unimodal relationship with the proportion of juvenile mussels (Figure 3a). The other tested factors had no relationship with the proportion of juvenile mussels (Figure 3b–d). The host fish density was affected only by the index of nutrient enrichment (Table S5, R2: 0.540), and this relationship was negative (Figure 3e).

4. Discussion

In the present study, we examined the relative effects of biotic and abiotic factors on the recruitment success of endangered M. laevis and established that host fish density was the most influential factor. In addition, we clarified the factors regulating the host fish density, which is closely related to the mussels’ life cycle. Freshwater mussels are among the most threatened groups of animals in freshwater ecosystems [38,39]. Improving recruitment will be the key for the conservation of this taxon. Our findings provide critical information for the prioritization of conservation measures for target freshwater mussel species.
The results revealed that the host fish density, which showed a unimodal relationship with the proportion of M. laevis juveniles, was the most important factor influencing the recruitment of M. laevis. At low to moderate host fish densities, the proportion of juveniles increased with host fish density. Such a positive relationship was also observed in previous studies (see, e.g., [17]) and may be attributable to the high infection rate at moderate host densities [16]. However, the proportion of juvenile mussels decreased from moderate to the highest host fish densities. One possible reason for this negative correlation is decreased larval growth in host fish that are in poor physical condition. Among stream-dwelling Salmo species, density-dependent competition has been well-documented. Previous studies showed that increasing fish density caused a decline in individual growth rates or body mass (see, e.g., [40,41]). Hasegawa et al. [42] clarified that intraspecific competition reduces the stomach contents and growth rate of our target species, O. masou masou. Decreased larval growth in M. margaritifera was also confirmed when the host fish could not grow sufficiently [43]. Hence, we hypothesize that poor host fish conditions were induced by high host fish density, which finally deteriorated mussel recruitment by reducing larval growth. Elucidating the details of this mechanism will require further research.
Abiotic factors had no direct impact on mussel recruitment. However, we should note that physicochemical environmental factors may indirectly affect recruitment by decreasing the host fish density (see, e.g., [24]). In this study, nutrient enrichment was found to have a negative effect on the host fish O. masou masou. Similar patterns were observed in Ishiyama et al. [44], in which the population of O. masou declined with increasing nutrient enrichment in the Hokkaido stream. In that study, using data on Ephemeroptera, Plecoptera, and Trichoptera (EPT) abundance, the authors suggested that the primary cause of the salmon population decline was the decline in the abundance of preferred food resources.
Understanding the distribution and age composition of populations is a key issue for the conservation of M. laevis. However, there have been very few studies that quantitatively investigated juveniles of M. laevis in Japan (but see [9,45]). Our quantitative study has a meaningful outcome and provides the following management implications for improving the health of M. laevis populations. We found that conserving host fish density is the most important measure for the recruitment of this species. In the study streams, water quality improvements would help to maintain a high host fish density because the index of nutrient enrichment negatively affected the host fish density. In many regions of Hokkaido, catchments have been converted to farmlands, which have decreased freshwater biodiversity (see, e.g., [46,47]). A positive relationship between the index of nutrient enrichment and the proportion of farmland in the basin was shown in our study streams (Figure S2). Therefore, excessive agricultural land use should be avoided to improve water quality and conserve host fish populations. Restoring riparian forests is an effective measure for improving water quality and promoting fish conservation because forests around streams filter nutrients supplied from catchments (see, e.g., [48,49]). On the other hand, the negative effect of excessive increases in host fish density must also be considered. The host fish, O. masou masou, is one of the key fishery resources in northern Japan, but excessive O. masou stocking programs may increase the risk of M. laevis recruitment impairment.
Few conservation measures targeting M. laevis have been implemented, and more efforts to promote the recruitment success of this species, as well as ex-post-evaluations of these efforts are needed [11]. Our findings regarding high-priority measures will contribute to the future management of the recruitment success of this species. However, none of the potential factors were examined in the present study. For example, the flow regime can greatly affect the recruitment success of M. laevis, because higher flows may cause bedload movement that is detrimental to juvenile mussels [50]. Further studies are needed to gain a deeper understanding of conservation measures. Moreover, we should note that it will take a long time (i.e., several decades) to accurately evaluate conservation outcomes for this species because the mortality rate of benthic juvenile mussels is high and because the shells of young juveniles are too small to observe for the first several years after they start to inhabit streambeds [30]. We believe that elucidating the relative effects of biotic and abiotic factors on mussel recruitment could be effective for understanding the ecology of and conservation priorities for other endangered freshwater mussels.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/w13091289/s1, Figure S1: Relationships between shell length and age in three regions, Figure S2: Relationship between the index of nutrient enrichment and proportion of farmland, Table S1: Abiotic and biotic factors used in the statistical analysis, Table S2: Host fish density in the study streams (A–J) from 2005 to 2016, Table S3: Results of preliminary analysis (GLMMs) for selecting explanatory variables from the host fish density and the index of nutrient enrichment, Table S4: Results of the estimations of four nonlinear growth models in the three regions, Table S5: Results from the GLMs examining the abiotic factors affecting the host fish density in ascending order of AICc.

Author Contributions

Conceptualization, K.K., N.I., K.M., A.T., M.S. and F.N.; methodology, K.K., N.I., K.M., A.T., M.S. and F.N.; software, K.K.; validation, K.K. and N.I.; formal analysis, K.K., N.I. and A.T.; investigation, K.K., N.I. and K.M.; resources, K.K. and N.I.; data curation, K.K. and N.I.; writing—original draft preparation, K.K. and N.I.; writing—review and editing, K.K., N.I., K.M., A.T., M.S. and F.N.; visualization, K.K.; supervision, K.K.; project administration, K.K., N.I. and F.N.; funding acquisition, K.K., N.I., K.M. and F.N. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Pro Natura Foundation Japan’s 27th Pro Natura Fund, JSPS Research Fellow Grant Number JP18J12458, JSPS KAKENHI Grant Number 18K18221, the Takara Harmonist Fund, and the Zoshinkai Fund in 2017 for Protection of Endangered Animal and Environment Research and Technology Development Fund (S15 Predicting and Assessing Natural Capital and Ecosystem Services [PANCES]) of the Ministry of the Environment, Japan.

Institutional Review Board Statement

Mussel sampling was undertaken with permission from the Fisheries Adjustment Rules of Hokkaido Prefecture.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data that support the findings of this study are available from the Supplementary Materials or corresponding author upon reasonable request.

Acknowledgments

We are grateful to H. Urabe of the Hokkaido Research Organization for sharing data. We thank the members of the Forest Ecosystem Management Laboratory and the members of the Watershed Conservation & Management Laboratory of Hokkaido University for their help in conducting the field and laboratory work. We also thank the members of the Teshio Experimental Forest, Field Science Center for the Northern Biosphere, Hokkaido University, for their support of our field survey.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Barnosky, A.D.; Matzke, N.; Tomiya, S.; Wogan, G.O.U.; Swartz, B.; Quental, T.B.; Marshall, C.; McGuire, J.L.; Lindsey, E.L.; Maguire, K.C.; et al. Has the Earth’s sixth mass extinction already arrived? Nature 2011, 471, 51–57. [Google Scholar] [CrossRef]
  2. O’Grady, J.J.; Reed, D.H.; Brook, B.W.; Frankham, R. What are the best correlates of predicted extinction risk? Biol. Conserv. 2004, 118, 513–520. [Google Scholar] [CrossRef]
  3. Hylander, K.; Ehrlén, J. The mechanisms causing extinction debts. Trends Ecol. Evol. 2013, 28, 341–346. [Google Scholar] [CrossRef] [PubMed]
  4. Myers, R.A.; Bowen, K.G.; Barrowman, N.J. Maximum reproductive rate of fish at low population sizes. Can. J. Fish. Aquat. Sci. 1999, 56, 2404–2419. [Google Scholar] [CrossRef]
  5. Aguilar, R.; Ashworth, L.; Galetto, L.; Aizen, M.A. Plant reproductive susceptibility to habitat fragmentation: Review and synthesis through a meta-analysis. Ecol. Lett. 2006, 9, 968–980. [Google Scholar] [CrossRef] [PubMed]
  6. Haag, W.R.; Rypel, A.L. Growth and longevity in freshwater mussels: Evolutionary and conservation implications. Biol. Rev. 2011, 86, 225–247. [Google Scholar] [CrossRef] [PubMed]
  7. Vaughn, C.C.; Hoellein, T.J. Bivalve impacts in freshwater and marine ecosystems. Annu. Rev. Ecol. Evol. Syst. 2018, 49, 183–208. [Google Scholar] [CrossRef] [Green Version]
  8. Hastie, L.C.; Young, M.R.; Boon, P.J.; Cosgrove, P.J.; Henninger, B. Sizes, densities and age structures of Scottish Margaritifera margaritifera (L.) populations. Aquat. Conserv. Mar. Freshw. Ecosyst. 2000, 10, 229–247. [Google Scholar] [CrossRef]
  9. Akiyama, Y. Factors causing extinction of a freshwater pearl mussel Margaritifera laevis in Japan (Bivalvia: Unionoida). Grad. Sch. Environ. Earth Sci. 2007. [Google Scholar] [CrossRef]
  10. Geist, J. Strategies for the conservation of endangered freshwater pearl mussels (Margaritifera margaritifera L.): A synthesis of conservation genetics and ecology. Hydrobiologia 2010, 644, 69–88. [Google Scholar] [CrossRef] [Green Version]
  11. Negishi, J.N.; Kayaba, Y.; Tsukahara, K.; Miwa, Y. Unionid mussels as imperiled indicator organisms: Habitat degradation processes and restoration approaches. Ecol. Civ. Eng. 2008, 11, 195–211, (In Japanese with English summary). [Google Scholar] [CrossRef]
  12. Ministry of the Environment Japan. Red Lists 2020. Available online: https://www.env.go.jp/press/files/jp/114457.pdf (accessed on 2 May 2021). (In Japanese).
  13. Anthony, J.L.; Kesler, D.H.; Downing, W.L.; Downing, J.A. Length-specific growth rates in freshwater mussels (Bivalvia: Unionidae): Extreme longevity or generalized growth cessation? Freshw. Biol. 2001, 46, 1349–1359. [Google Scholar] [CrossRef] [Green Version]
  14. Wächtler, K.; Dreher-Mansur, M.C.; Richter, T. Larval Types and Early Postlarval Biology in Naiads (Unionoida). In Ecology and Evolution of the Freshwater Mussels Unionoida; Bauer, G., Wächtler, K., Eds.; Springer: Berlin/Heidelberg, Germany, 2001; pp. 93–125. [Google Scholar] [CrossRef]
  15. Geist, J.; Auerswald, K. Physicochemical stream bed characteristics and recruitment of the freshwater pearl mussel (Margaritifera margaritifera). Freshw. Biol. 2007, 52, 2299–2316. [Google Scholar] [CrossRef]
  16. Österling, M.E.; Greenberg, L.A.; Arvidsson, B.L. Relationship of biotic and abiotic factors to recruitment patterns in Margaritifera margaritifera. Biol. Conserv. 2008, 141, 1365–1370. [Google Scholar] [CrossRef]
  17. Arvidsson, B.L.; Karlsson, J.; Österling, M.E. Recruitment of the threatened mussel margaritifera margaritifera in relation to mussel population size, mussel density and host density. Aquat. Conserv. Mar. Freshw. Ecosyst. 2012, 22, 526–532. [Google Scholar] [CrossRef]
  18. Modesto, V.; Ilarri, M.; Souza, A.T.; Lopes-Lima, M.; Douda, K.; Clavero, M.; Sousa, R. Fish and mussels: Importance of fish for freshwater mussel conservation. Fish Fish. 2018, 19, 244–259. [Google Scholar] [CrossRef]
  19. Terui, A.; Miyazaki, Y. A “parasite-tag” approach reveals long-distance dispersal of the riverine mussel Margaritifera laevis by its host fish. Hydrobiologia 2015, 760, 189–196. [Google Scholar] [CrossRef]
  20. Bauer, G. Threats to the freshwater pearl mussel Margaritifera margaritifera L. in Central Europe. Biol. Conserv. 1988, 45, 239–253. [Google Scholar] [CrossRef]
  21. Österling, M.E.; Arvidsson, B.L.; Greenberg, L.A. Habitat degradation and the decline of the threatened mussel Margaritifera margaritifera: Influence of turbidity and sedimentation on the mussel and its host. J. Appl. Ecol. 2010, 47, 759–768. [Google Scholar] [CrossRef]
  22. Varandas, S.; Lopes-Lima, M.; Teixeira, A.; Hinzmann, M.; Reis, J.; Cortes, R.; Machado, J.; Sousa, R. Ecology of southern European pearl mussels (Margaritifera margaritifera): First record of two new populations on the rivers Terva and Beça (Portugal). Aquat. Conserv. Mar. Freshw. Ecosyst. 2013, 23, 374–389. [Google Scholar] [CrossRef]
  23. Strayer, D.L.; Malcom, H.M. Causes of recruitment failure in freshwater mussel populations in southeastern New York. Ecol. Appl. 2012, 22, 1780–1790. [Google Scholar] [CrossRef] [PubMed]
  24. Österling, M.E.; Högberg, J.O. The impact of land use on the mussel Margaritifera margaritifera and its host fish Salmo trutta. Hydrobiologia 2014, 735, 213–220. [Google Scholar] [CrossRef]
  25. Sousa, R.; Amorim, Â.; Sobral, C.; Froufe, E.; Varandas, S.; Teixeira, A.; Lopes-Lima, M. Ecological Status of a Margaritifera margaritifera (Linnaeus, 1758) Population at the Southern Edge of its Distribution (River Paiva, Portugal). Environ. Manag. 2013, 52, 1230–1238. [Google Scholar] [CrossRef] [PubMed]
  26. Kobayashi, O.; Kondo, T. Age Determination of the Freshwater Pearl Mussel Margaritifera laevis (Bivalvia: Margaritiferidae) in the Chubu-Nougu River, Nagano Prefecture. Venus 2008, 67, 61–71, (In Japanese with English summary). [Google Scholar] [CrossRef]
  27. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef]
  28. Kondo, T.; Kobayashi, O. Revision of the Genus Margaritifera (Bivalvia: Margaritiferidae) of Japan, with Description of a New Species. Venus 2005, 64, 135–140. [Google Scholar] [CrossRef]
  29. Miura, K.; Ishiyama, N.; Kawajiri, K.; Atsumi, K.; Tachibana, M.; Nagasaka, Y.; Machida, Y.; Yiyang, G.; Negishi, J.N.; Koizumi, I.; et al. Simple nonlethal identification criteria for two endangered freshwater pearl mussels, Margaritifera laevis and Margaritifera togakushiensis, in Hokkaido, northern Japan. Ecol. Res. 2019, 34, 667–677. [Google Scholar] [CrossRef]
  30. Kondo, T. Monograph of Unionoida in Japan (Mollusca: Bivalvia); Special Publication of the Malacological Society of Japan: Tokyo, Japan, 2008. [Google Scholar]
  31. Hastie, L.C.; Young, M.R.; Boon, P.J. Growth characteristics of freshwater pearl mussels, Margaritifera margaritifera (L.). Freshw. Biol. 2000, 43, 243–256. [Google Scholar] [CrossRef]
  32. Nakamura, K.; Cucala, L.; Mestre, A.; Mesquita-Joanes, F.; Elbaile, E.; Salinas, C.; Muñoz-Yanguas, M.Á. Modelling growth in the critically endangered freshwater mussel Margaritifera auricularia (Spengler, 1793) in the Ebro basin. Hydrobiologia 2018, 810, 375–391. [Google Scholar] [CrossRef]
  33. San Miguel, E.; Monserrat, S.; Fernández, C.; Amaro, R.; Hermida, M.; Ondina, P.; Altaba, C.R. Growth models and longevity of freshwater pearl mussels (Margaritifera margaritifera) in Spain. Can. J. Zool. 2004, 82, 1370–1379. [Google Scholar] [CrossRef]
  34. Urban, H.J. Modeling growth of different developmental stages in bivalves. Mar. Ecol. Prog. Ser. 2002, 238, 109–114. [Google Scholar] [CrossRef] [Green Version]
  35. Bates, D.; Maechler, M.; Bolker, B.; Walker, S.; Christensen, R.H.B.; Singmann, H. lme4: Linear mixed-effects models using Eigen and S; R Package Version 1.4. 2014. Available online: https://www.researchgate.net/publication/279236477_Package_Lme4_Linear_Mixed-Effects_Models_Using_Eigen_and_S4 (accessed on 2 May 2021).
  36. Barton, K. MuMIn: Multi-Model Inference. R Package Version 1. 0. 0. 2009. Available online: https://cran.r-project.org/web/packages/MuMIn/index.html (accessed on 2 May 2021).
  37. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical 527 Computing: Vienna, Austria, 2020; Available online: https://www.R-project.org/ (accessed on 2 May 2021).
  38. Ricciardi, A.; Rasmussen, J.B. Extinction Rates of North American Freshwater Fauna. Conserv. Biol. 1999, 13, 1220–1222. [Google Scholar] [CrossRef]
  39. Haag, W.R.; Williams, J.D. Biodiversity on the brink: An assessment of conservation strategies for North American freshwater mussels. Hydrobiologia 2014, 735, 45–60. [Google Scholar] [CrossRef]
  40. Jenkins, T.M.; Diehl, S.; Kratz, K.W.; Cooper, S.D. Effects of population density on individual growth of brown trout in streams. Ecology 1999, 80, 941–956. [Google Scholar] [CrossRef]
  41. Grant, J.W.A.; Imre, I. Patterns of density-dependent growth in juvenile stream-dwelling salmonids. J. Fish Biol. 2005, 67, 100–110. [Google Scholar] [CrossRef]
  42. Hasegawa, K.; Morita, K.; Ohkuma, K.; Ohnuki, T.; Okamoto, Y. Effects of hatchery chum salmon fry on density-dependent intra- and interspecific competition between wild chum and masu salmon fry. Can. J. Fish. Aquat. Sci. 2014, 71, 1475–1482. [Google Scholar] [CrossRef]
  43. Österling, M.E.; Larsen, B.M. Impact of origin and condition of host fish (Salmo trutta) on parasitic larvae of Margaritifera margaritifera. Aquat. Conserv. Mar. Freshw. Ecosyst. 2013, 23, 564–570. [Google Scholar] [CrossRef]
  44. Ishiyama, N.; Miura, K.; Inoue, T.; Sueyoshi, M.; Nakamura, F. Geology-dependent impacts of forest conversion on stream fish diversity. Conserv. Biol. 2020, cobi.13655. [Google Scholar] [CrossRef]
  45. Terui, A.; Miyazaki, Y.; Yoshioka, A.; Kaifu, K.; Matsuzaki, S.I.S.; Washitani, I. Asymmetric dispersal structures a riverine metapopulation of the freshwater pearl mussel Margaritifera laevis. Ecol. Evol. 2014, 4, 3004–3014. [Google Scholar] [CrossRef] [PubMed]
  46. Ishiyama, N.; Akasaka, T.; Nakamura, F. Mobility-dependent response of aquatic animal species richness to a wetland network in an agricultural landscape. Aquat. Sci. 2014, 76, 437–449. [Google Scholar] [CrossRef]
  47. Sueyoshi, M.; Ishiyama, N.; Nakamura, F. β-diversity decline of aquatic insects at the microhabitat scale associated with agricultural land use. Landsc. Ecol. Eng. 2016, 12, 187–196. [Google Scholar] [CrossRef]
  48. Nakamura, F.; Yamada, H. Effects of pasture development on the ecological functions of riparian forests in Hokkaido in northern Japan. Ecol. Eng. 2005, 24, 539–550. [Google Scholar] [CrossRef]
  49. Valle Junior, R.F.; Varandas, S.G.P.; Pacheco, F.A.L.; Pereira, V.R.; Santos, C.F.; Cortes, R.M.V.; Sanches Fernandes, L.F. Impacts of land use conflicts on riverine ecosystems. Land Use Policy 2015, 43, 48–62. [Google Scholar] [CrossRef] [Green Version]
  50. Vaughn, C.C.; Taylor, C.M. Impoundments and the Decline of Freshwater Mussels: A Case Study of an Extinction Gradient. Conserv. Biol. 1999, 13, 912–920. [Google Scholar] [CrossRef]
Figure 1. Locations of the study streams where quantitative sampling was performed (gray circles) and the sampling areas where mussels were collected for age estimation (circles) in Hokkaido. The name and specific locations of the study streams are withheld to protect the habitat of M. laevis.
Figure 1. Locations of the study streams where quantitative sampling was performed (gray circles) and the sampling areas where mussels were collected for age estimation (circles) in Hokkaido. The name and specific locations of the study streams are withheld to protect the habitat of M. laevis.
Water 13 01289 g001
Figure 2. Shell length distribution in each study stream. Streams (AH) are in northern Hokkaido, stream (I) is in central Hokkaido, and stream (J) is in eastern Hokkaido. The name and specific locations of the study streams are withheld to protect the habitat of M. laevis. The percentage next to the stream name is the percentage of juvenile mussels in each population.
Figure 2. Shell length distribution in each study stream. Streams (AH) are in northern Hokkaido, stream (I) is in central Hokkaido, and stream (J) is in eastern Hokkaido. The name and specific locations of the study streams are withheld to protect the habitat of M. laevis. The percentage next to the stream name is the percentage of juvenile mussels in each population.
Water 13 01289 g002
Figure 3. Relationships between the proportions of juvenile mussels and investigated factors: host fish density, water depth, proportion of fine sediment and summer maximum water temperature (ad), and the relationship between host fish density and index of nutrient enrichment (e). The points represent the observed values. The solid line represents the predicted values (Table 1 and Table S5).
Figure 3. Relationships between the proportions of juvenile mussels and investigated factors: host fish density, water depth, proportion of fine sediment and summer maximum water temperature (ad), and the relationship between host fish density and index of nutrient enrichment (e). The points represent the observed values. The solid line represents the predicted values (Table 1 and Table S5).
Water 13 01289 g003
Table 1. Results of the averaged models examining the relative effects of factors affecting the recruitment success of M. laevis. Significant results according to the significance test (α = 0.05) are indicated with an asterisk.
Table 1. Results of the averaged models examining the relative effects of factors affecting the recruitment success of M. laevis. Significant results according to the significance test (α = 0.05) are indicated with an asterisk.
Explanatory VariableCoefficient95% Confidence
Interval
Significancep-Value
Log (Host fish density)−5.23−8.74-−1.73*0.003
Log (Host fish density)2−5.31−8.78-−1.83*0.003
Water depth−0.04−0.22-0.15 0.919
(Water depth)2−0.05−0.22-0.13 0.857
Fine sediment−0.22−0.53-0.10 0.462
(Fine sediment)2−0.19−0.44-0.06 0.462
Maximum water temperature−1.35−7.12-4.42 0.732
(Maximum water temperature)20.22−7.18-7.62 0.971
Intercept−2.32−3.06-−1.57*<0.001
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kawajiri, K.; Ishiyama, N.; Miura, K.; Terui, A.; Sueyoshi, M.; Nakamura, F. The Relative Effects of Biotic and Abiotic Factors on the Recruitment of Freshwater Mussels (Margaritifera laevis). Water 2021, 13, 1289. https://doi.org/10.3390/w13091289

AMA Style

Kawajiri K, Ishiyama N, Miura K, Terui A, Sueyoshi M, Nakamura F. The Relative Effects of Biotic and Abiotic Factors on the Recruitment of Freshwater Mussels (Margaritifera laevis). Water. 2021; 13(9):1289. https://doi.org/10.3390/w13091289

Chicago/Turabian Style

Kawajiri, Keita, Nobuo Ishiyama, Kazuki Miura, Akira Terui, Masanao Sueyoshi, and Futoshi Nakamura. 2021. "The Relative Effects of Biotic and Abiotic Factors on the Recruitment of Freshwater Mussels (Margaritifera laevis)" Water 13, no. 9: 1289. https://doi.org/10.3390/w13091289

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