1. Introduction
The effectiveness of stock enhancement programs is significantly influenced by the post-release performance of cultured fish. While large-scale stock enhancements are commonly conducted worldwide, successful cases have been limited [
1]. One of the reasons for this failure is the insufficient understanding of the biology and ecology of the stocks targeted for enhancement [
2]. Additionally, studies on the performance and behavior of fish stocks have been scarce and more challenging compared to research on their general biological characteristics such as growth, food habits, and reproduction. Previous reports on the post-release performance of fish have primarily focused on adult or large-sized fish species, such as
Acipenser sturio,
Argyrosomus japonicus,
Anguilla japonica,
Etelis coruscans, and
Myxocyprinus asiaticus [
3,
4,
5,
6,
7]. In China, numerous stock enhancement programs have been implemented in recent decades, with their scale expanding each year. However, most of these programs in China have targeted small-sized juvenile cultured fish, highlighting the need for further research on the post-release performance of small-sized cultured fish.
There is currently no recognized standard for what indicators to use and how long the research period should be in order to accurately evaluate the performance of released fishes. There are various indicators utilized to assess the post-release performance of cultured fishes, although they lack a unified approach. Previous research has employed indicators such as growth [
8,
9], survival [
10], dispersal or behavior [
5,
6], and physiology [
11] to evaluate post-release performance. In addition, conducting a long-term recapture study is challenging and labor-intensive. Typically, the initial stage of release is of utmost importance in determining the viability of cultured fish, as the majority of mortalities occur within a short timeframe following their release [
12]. Furthermore, monitoring and assessing their initial performance is considerably more feasible.
Mark–release–recapture is the fundamental method used to study the post-release performance of artificially released populations. Common fish tagging methods include listing, visible implant fluorescence (VIE), coded wire tags (CWT), and acoustic telemetry. The selection of appropriate tagging methods should be based on factors such as fish species, fish size, tagging cost, and characteristics of the study area. Due to its advantages of individual recognition and no need to recapture, acoustic telemetry was widely used in studies about the post-release movement of fishes, such as applications on
Acipenser sturio [
3],
Argyrosomus japonicus [
4],
Anguilla japonica [
5],
Etelis coruscans [
6], and
Myxocyprinus asiaticus [
7]. However, the acoustic telemetry was expensive and had large tags, which made it unsuitable for small-sized juvenile fish. Alternatively, VIE and CWT were recommended to the mark–release–recapture studies of small-sized fish as they were low cost and did minimal harm. For example, Leber et al. [
13] reported on the recruitment patterns of juvenile, cultured Pacific threadfin (
Polydactylus sexfilis) from the Polynemidae family. The study released the fish along sandy marine shores in Hawaii using CWT. Hervas et al. [
9] also revealed the growth, survival, and distribution of released
Atractoscion nobilis using the CWT method. The initial post-release movement of cultured
Cyprinus pellegrini juveniles was evaluated using the VIE method [
14].
Cyprinus chilia is endemic to the plateau lakes in Yunnan Province, China [
15].
Cyprinus chilia is omnivorous and mainly feeds on large benthic invertebrates. Under natural conditions,
C. chilia can grow up to 4 kg and is an important economic fish [
15]. However, the natural resource of
C. chilia has declined rapidly over the last decades due to overfishing, deterioration of the water environment, habitat loss, invasion of alien species, and other factors. At present, it is difficult to find wild populations of
C. chilia in its naturally distributed lakes. To protect this endangered species, artificial breeding [
16] and stock enhancement of
C. chilia have been successfully implemented. However, there is a lack of information regarding the post-release performance and evaluation of these stock enhancement efforts. Also, a suitable tagging method for
C. chilia juveniles has not been evaluated yet. Therefore, the management and planning of
C. chilia stock enhancement programs remain unclear.
VIE was one of the most popular tagging methods used to mark released fish in China, but there are few reports about the effects of VIE tagging on the assessment of stock enhancement. Therefore, the objectives of this study are (1) to assess the effectiveness of VIE as a short-term tagging method for C. chilia juveniles and (2) to investigate the initial post-release performance of cultured C. chilia juveniles by monitoring their movement, spatial distribution, gut fullness, and gut microbiota. Correspondingly, two hypotheses are proposed: (1) the VIE should be suitable for the short-term tagging of C. chilia juveniles, and (2) cultured C. chilia juveniles should go through an obvious transitional period in the new environment. The findings of this study will contribute to a better understanding of the post-release performance of cultured Cyprinus fishes and will be valuable for the management and evaluation of stock enhancement programs related to C. chilia.
2. Materials and Methods
2.1. Fish
The fish were artificially hatched C. chilia juveniles produced at the Raofu fish farm, Tonghai County. Before the mark–release–recapture trial, fish were about 4 months old and had a mean total length (L) of 10.8 ± 1.2 cm (mean ± standard deviation, n = 46) and mean body weight (W) of 17.54 ± 5.19 g.
2.2. Mark–Release–Recapture Trial
A mark–release–recapture trial was conducted to examine the movement, spatial distribution, gut fullness, and gut microbiota of C. chilia juveniles after release. VIE tags were used in the study as its low cost and minimal harm to small-sized fishes.
A total of 20,000 C. chilia juveniles were selected for the mark–release–recapture trial. All of the selected fish were tagged with VIE tags to help the identification of fish samples in the subsequent recapture study. The tagging work was carried out on 25–26 July 2022. Prior to the tagging operation, fish were anesthetized using 30 mg/L of MS-222 for 1~2 min. VIE tagging procedures were as follows: red VIEs were injected into the epidermis of the head skin using hand-pressurized syringes (Northwest Marine Technology, Anacortes, WA, USA). The presence of the tags was visually confirmed. Following the tagging process, the fish were transferred to clean water to facilitate their recovery (recovery rate 100%).
The designated water area for release is the Qilu Lake in Yunnan Province, China (
Figure 1), as it is one of the main distribution lakes of
C. chilia. The Qilu Lake has a maximum length of 10.4 km, maximum width of 4.4 km, water area of 38.86 km
2, a maximum water depth of 540 cm, and an average water depth of 403 cm. There were 11 species of fish historically distributed in Qilu Lake, including
C. chilia,
Cyprinus pellegrini,
Cyprinus ilishaestomus,
Cyprinus yunnanensis,
Anabarilius qiluensis,
Carassius auratus,
Misgurnus anguillicaudatus,
Oryzias latipes,
Channa argus,
Monopterus albus,
Silurus grahami [
17]. There has been no record of
C. chilia in Qilu Lake since 2000s. According to our investigation of the entire Qilu Lake on 12–14 July 2022, the present fish species include
Cyprinus carpio,
Carassius auratus,
Abbottina rivularis,
Hemicculter Leuciclus,
Pseudorasbora parva,
Cultrichthys erythropterus,
Pelteobagrus fulvidraco,
Misgurnus anguillicaudatus,
Paramisgurnus dabryanus,
Rhinogobius giurinus,
Hyporhamphus intermedius. The release site was situated on the south shore of Qilu Lake, with a water depth of 220 cm, transparency of 37 cm, water temperature of 27.4 °C, dissolved oxygen level of 10.46 mg/L, a pH of 8.83, and conductivity of 891 μs/cm. On 27 July 2022, the tagged fish were packed into double-layered nylon fish bags (80 cm × 40 cm) with a density of 150~200 fish per bag. The volume ratio of oxygen to water in the bag was approximately 3:1. The fish bags were then transported (approximately 0.5 h) to the release site for release. The entire release process was conducted manually. The fish bags were initially submerged in the water at the lakeshore to balance the temperature for 0.5 h. Afterward, the fish bags were untied, and the tagged fish were carefully released into the lake.
The recapture work was divided into two stages, including a 7-day continuous recapture and a long-term recapture.
The 7-day continuous recapture was conducted immediately following the release of fish. The lake was divided into key recapture area and regular recapture area. The key recapture area was around the release site (
Figure 2a). A total of 13 fixed recapture sites were set up in the key recapture area, including 1 site (S1) at the release site, 4 sites (E1, E2, E3, E4) on the eastern lakeshore, 4 sites (W1, W2, W3, W4) on the western lakeshore, and 4 sites (N1, N2, N3, N4) in the northern lake area (
Figure 2a). This created a 4-layer arc recapture monitoring network, with the release point as the origin. The distances from the recapture sites on the 1st, 2nd, 3rd, and 4th layers to the release site were set as 150, 300, 500, and 800 m, respectively. However, due to limitations in terrain, interference from wind and waves on the positioning ship, and other factors, there was some variation between the actual distance between adjacent recapture sites and the preset values (
Table 1 and
Table 2). The regular recapture area included lake areas outside of key recapture area (
Figure 2b). A total of 12 fixed recapture sites were set up in the regular recapture area, including sites D1, D2, D3, D4, D5, D6, D7, D8, D9, D10, D11, and D12 (
Figure 2b).
The long term recapture was conducted one month and three months after release. The recapture sites for the long term recapture was the same as the recapture sites in the regular recapture area of the 7-day continuous recapture, plus the release sites (S1). For each long-term recapture, 3-day continuous recapture was conducted.
Recapture gears used were trap nets with dimensions of 10 m in length, 45 cm in width, 33 cm in height, and a mesh size of 7 mm. The maximum volume of the net was 1.485 m3. A density of one net per site was set, and the nets were deployed 2 h prior to fish release. After the fish were released, the catch in the nets was collected every 24 h, and the nets were then repositioned during the recapture periods. It should be noted that the net at S1 had a water entry time of only 6 h on the fifth day after release, as it was affected by an incorrect removal of the net.
2.3. Processing of the Recaptured Fish Samples
All recaptured fish samples were identified by species. For C. chilia juveniles, the tag presence of fish was determined through visual examination. The recaptured C. chilia juveniles were counted separately based on the date and recapture site. When fewer than 30 released fish samples were recaptured at a site per day, the total length and weight of all fish samples were measured. If more than 30 released fish were recaptured at a site per day, a random sample of 30 fish was measured for total length (accurate to 0.1 cm) and weight (accurate to 0.1 g).
According to the number of fish samples, 3–20 fish were randomly selected for dissection every day to assess the feeding intensity of the released fish by examining their gut fullness levels [
18]. The level of gut fullness ranged from 0 to 5, representing feed intake from low to high (
Table 3). Additionally, freshly dissected guts of three fish samples were collected every day for analysis of gut microbiota structure. According to the post-release days and sampling order, the samples were named as D0–1, D0–2, D0–3, D1–1, D1–2, D1–3, D2–1, D2–2, D2–3, D3–1, D3–2, D3–3, D4–1, D4–2, D4–3, D5–1, D5–2, D5–3, D6–1, D6–2, D6–3, D7–1, D7–2, and D7–3, respectively. After sampling, any remaining surviving fish were transferred to indoor breeding tanks for temporary cultivation.
2.4. Gut Microbiota Structure Analysis
Total microbiota genomic DNA was extracted from approximately 0.3 g of the freshly dissected intestine of each fish using the FastDNA spin kit for soil (MP, Solon, OH, USA) following the manufacturer’s protocols. The 515F (5′-GTGYCAGCMGCCGCGGTA-3′) and 907R (5′-CCGTCAATTCMTTTRAGT-3′) primers [
19,
20] were used to amplify the bacterial 16S rRNA gene V4–V5 fragments. PCR integration and protocols were carried out as follows: 94 °C for 3 min followed by 30 cycles of 94 °C for 40 s, 56 °C for 60 s, 72 °C for 60 s, and a final extension at 72 °C for 10 min until the reaction was halted by the user. The amplicons were purified and sequenced using the Illumina Miseq platform at Guangdong Meilikang Bio-Science Ltd. (Foshan, China).
Raw reads were merged using FLASH 1.2.8 and subsequently processed using QIIME 1.9.0, as described [
21]. Briefly, all the merged sequences were assigned to each sample based on their barcode sequences, and the trimmed barcodes and primer sequences were removed using QIIME 1.9.0 software. Low-quality and chimeric sequences were removed using QIIME 1.9.0 and UCHIME 4.2.40, respectively. Subsequently, the remaining high-quality sequences were clustered into operational taxonomic units (OTUs) with 97% identity using UPARSE 7.1 [
22]. The taxonomy of each OTU was assigned using the RDP classifier [
23] in the gg_13_8 database. Alpha diversity indices were calculated using QIIME version 1.9.0. The Illumina sequencing raw data obtained from this study were deposited in the NCBI Sequence Read Archive with accession No. PRJNA999371.
2.5. Data Analysis
The condition factor, recapture rate, tag retention rate, fish distribution density, water depth selectivity, and daily movement velocity were calculated.
Condition factor:
where
CF was condition factor,
W was body weight (g), and
L was total length (cm).
Recapture rate:
where
C was recapture rate,
n was the number of recaptured releasing fish, and
N was the number of releasing fish.
Tag retention rate:
where
R was the tag retention rate,
np was the number of recaptured
C. chilia with VIE tag, and
n was the total number of recaptured
C. chilia.
Fish density:
where
Di was the density (ind./m
3) of fish at recapture site
i,
Ni was the number of recaptured releasing fish per day at recapture site
i, and
Vi was the fully extended volume of the trap net used at recapture site
i.
To understand any preference of released fish for water depth, a water depth selectivity index was calculated. The recapture sites were divided into 4 water depth groups based on their respective water depths (≤250 cm, 251–350 cm, 351–450 cm, ≥451 cm). The water depth selectivity index was calculated using the following function:
where
Sj was the selectivity of the releasing fish for the water depth group
j,
Nj was the mean number of recaptured fish at the recapture sites from the water depth group
j,
m was the number of water depth groups. The range of water depth selectivity values is 0~1 (0 means complete avoidance, 1 means total selectivity).
The daily movement velocity of the released fish in different directions (east, west, and north) in the first 7 days after release was calculated to assess the movement pattern of releasing fish using the following function:
where
V was mean daily movement velocity of the releasing fish towards a certain direction,
Li was the linear distance from the release site (S1) to site
i in a certain direction,
Nij was the number of recaptured releasing fish at recapture site
i on day
j, and
Tj was the number of days post release.
The α-diversity indices of the gut microbiota were determined using the richness index, Shannon index, Simpson index, and Chao1 index of observed OTUs in each sample.
The difference in individual size and α-diversity indexes of the gut microbiota among groups was compared using one-way ANOVA. A correlation heatmap of dominant microbial genera was analyzed using the corrplot R package. Beta diversity was calculated through unconstrained principal coordinate analysis (PCoA) based on weighted UniFrac distance to display the differences in gut microbial community structure among different sampling days. A significant difference was set at a value of p < 0.05.
4. Discussion
Marking methods suitable for juvenile fish are crucial to improving the efficiency of behavioral studies on released fish and evaluating the effects of numerous ongoing stock enhancement programs. Previous studies on the post-release behaviors of fish have primarily relied on acoustic telemetry [
4,
5,
6,
7,
24] and satellite archival tags [
25]. These methods enable the monitoring of individual fish movement and habitat selection without the need for physical capture. However, these marking methods had large tags and were therefore unsuitable for smaller-sized juvenile fish. Furthermore, the cost associated with acoustic telemetry and satellite archival tags makes it impractical to tag a large number of fish. It is worth noting that electronic tags may also affect the short-term behavior of fish [
26]. Consequently, marking methods such as VIE and CWT are preferred for mark–release–recapture studies of large quantities of small-sized fishes due to their affordability, minimal damage, ease of operation, and ease of identification. The findings of the present study demonstrate high VIE tag retention and minimal damage in all recaptured fish samples, indicating that VIE is a suitable marking method for short-term mark–release–recapture studies of
C. chilia.
Many studies have confirmed that the hatchery environment can affect the post-release behavior of fishes [
27]. On the first day after release, the majority of individuals remained near the release point, which is referred to as the “residency period” [
4]. The localized behavior exhibited by the cultured
C. chilia juveniles on the first day after release might be due to unfamiliarity in a new environment, or it could be an extension of the settlement behavior developed in the artificial rearing environment. Similarly, hatchery-reared Mulloway (
Argyrosomus japonicus) [
4] and
Micropterus cataractae [
28] also showed a period of minimal movement when released into new environments with varying durations. From the second day after release, the
C. chilia juveniles began to increase their movement and exploration of the environment. This may reflect the tradeoff between avoiding predation and resource utilization in the new environment. As time passes, factors such as hunger may motivate
C. chilia juveniles to leave the release point and spread to more unfamiliar waters in order to increase their opportunity to obtain resources.
The released fish showed significant habitat selectivity. As an important parameter of fish habitat selectivity, water depth preference were reported in many fish species. For example,
Acanthopagrus schlegelii prefers water depths of 8–10 m, while
Lutjanus argentimaculatus prefers a water depth of 1–3.5 m [
29]. The movements of land-locked Atlantic salmon (
Salmo salar) in a large lake were nearshore (<2 km from shore) from spring to summer at ~20 m bathymetric depths [
30]. The present study found the released fish mainly remained in shallow waters near the release point. Selectivity of water depth showed that released fish preferred water depths of up to 300 cm. Furthermore, the movement of released fish towards the eastern lakeshore was greater in scale than that towards the western lakeshore. It is speculated that the difference in complexity of the lake shorelines in the east and west may be related. The shoreline near the western side of the lake was relatively straight, while the shoreline near the eastern side of the lake was winding and complex. This complex shoreline provides more diverse habitat conditions, which can be beneficial for fish by providing shelter and feeding opportunities. Other species, like
Acanthopagrus schlegelii, also prefer to stay in areas where food is more abundant and easy to hide in the initial stages after release [
31].
The post-release movement velocity of
C. chilia juveniles fluctuated over the 7-day trial. The movement of
C. chilia juveniles decreased over the first three days after release and remained low for the following four days. This suggests that the majority of movement by
C. chilia juveniles occurs within the first three days after release. The dispersal or movement of different fish often varies after release. For example, the horizontal movement velocity of released deep-water longtail red snapper (
Etelis coruscans) was 2.2 km/d [
6], while the movement velocity of artificially bred Japanese eels (
Anguilla japonica) after release into the sea can reach 2.31 km/h [
5]. Lake trout (
Salvelinus namaycush) released into Lake Ontario had an average swimming speed of 1.64 km/h over the first day [
32]. Cyprinus fishes are generally sedentary. For example, the mean daily movements of grass carp (
Ctenopharyngodon idella) in Lake Erie ranged from <0.01 to 2.49 km/d, and only 25% of studied grass carp had mean daily movements greater than 0.88 km/d [
33]. Therefore, the movement of
C. chilia juveniles in our study was also relatively slow. The maximum recorded movement velocity in the present study was 813 m/d, and the overall movement gradually decreased over time, which is similar to that of some other fish species. For example, the movement velocity of mature Chinese sucker (
Myxocyprinus asiaticus) released into the Yangtze River also showed a declining trend over time [
7]. The gradual adaptation to the new environment and the acquisition and retention of an ideal habitat space may be reasons for the declining daily movement distance of
C. chilia juveniles in this study.
The levels of gut fullness in released
C. chilia juveniles showed a trend of initially decreasing and then increasing, although no statistical differences were found among recapture days. Limited sample sizes may explain these results. However, the mean levels of gut fullness in the released fish decreased from day 1 to day 4 after release, and then increased from day 4 to day 7 after release. This suggests that the released fish may begin feeding on the 5th day after release. Based on changes in the movements of released fish, it may be inferred that the first four days post-release are a transitional period for
C. chilia juveniles in adapting to the new environment. During this period, they need to find suitable habitats and adapt to consuming natural prey. By the fifth day after release, released fish have entered an adaptation period during which they have relatively stable habitats and feeding conditions. Similarly, the gut fullness of the hatchery-reared honmasu salmon
Oncorhynchus rhodurus ×
masou parr, released a week prior, was found to be comparable to that of wild parr [
11].
This study indicates that the structure of the gut microbial community in
C. chilia juveniles undergoes significant changes after they were released. Our study provided the first report on the gut microbiota of
C. chilia. From the perspective of the composition of gut microbial phyla, the variations within 7 days after release were not significant. However, there is a fluctuation in the dominant genus composition of gut microbes in
C. chilia juveniles after release. The dominant genera of gut microbiota on the first day after release mainly included Firmicutes other, Microcystis PCC-7914,
Cetobacterium, PeM15_unclassified genus, and Planktothrix NIVA-CYA 15. The members of Firmicutes and
Cetobacterium perform digestion functions. For example, PeM15_unclassified genus is a member of the Actinobacteria phylum, which has been verified as an important group of PAOs (phosphorus-accumulating organisms) in enhanced biological phosphorus removal systems. Additionally, some members of this genus may contain nitrite reductase genes that are involved in denitrification [
34]. Microcystis PCC-7914 and Planktothrix NIVA-CYA 15 are indicators of water body eutrophication, which can be used to assess the eutrophication levels in ponds used for fish rearing. Starting from the second day after release, the dominant genera continued to change, indicating that the
C. chilia juveniles were undergoing a transitional period of adaptation to the new environment. Generally, the dominant genera after release were
Aeromonas and
Romboutsia.
Aeromonas can secrete enzymes related to pathogenicity and environmental adaptation, such as hemolytic enterotoxin, lipase, protease, and amylase [
35]. These enzymes improve the digestibility of food and play an important role in the digestion process.
Romboutsia can ferment and metabolize macromolecular carbohydrates that are difficult for a host fish to digest into short-chain fatty acids, such as butyric acid. This process helps reduce intestinal pH, improve the host’s immune regulation ability, and maintain the balance of intestinal microecology [
36]. The increase in these two genera may be an adaptive strategy of
C. chilia juveniles to cope with the increased risk of food scarcity and stress infection after release.
Only the 7-day continuous recapture collected the released fish successfully, and the reasons for the failure to obtain samples over a longer period (1 month and 3 months after release) are unknown. We speculate that the dispersal of C. chilia after release may lead to a continuous decrease in density, thereby reducing the probability of catching released fish, or perhaps they experienced high post-release mortality. However, the present results may be sufficient to reveal the initial post-release performance of cultured C. chilia juveniles.