The Sympatric Coexistence Mechanism: A Case Study of Two Penahia Species in the Beibu Gulf, South China Sea

Simple Summary In this study, the feeding habits, trophic niches, and spatial niches of Pennahia pawak and Pennahia anea were investigated with the aim of exploring the coexistence and competition mechanism of these two species in the Beibu Gulf. The results indicate that both P. pawak and P. anea exhibit feeding shifts and differentiation in their trophic and spatial niches. Abstract The study of trophic relationships among closely related species plays an important role in deepening our understanding of the resource utilization characteristics, differentiation patterns, and population dynamics of co-occurring species in the same habitat. This research uses two congeneric fish species, Pennahia pawak and Pennahia anea, as examples. Based on a stomach content analysis and a carbon–nitrogen stable isotope analysis, a comparative analysis of their feeding habits and trophic niches is conducted. Additionally, a spatial niche analysis is employed to explore the coexistence and competitive mechanisms between these two closely related fish species. The results show that specialized feeding habits mitigate intraspecific competition as the population densities increase. The carbon and nitrogen stable isotope analysis reveals variations in the feeding habits and trophic levels with body length, indicating adaptive shifts in prey selection. Despite similar food resources, niche differentiation arises due to differences in dominant prey, facilitating coexistence. Differences in spatial niche further contribute to niche separation and coexistence. In resource-limited environments, species such as Pennahia utilize trophic and spatial niche differentiation to collectively exploit resources and achieve coexistence, with implications for fishery management favoring Pennahia resource occupancy capabilities.


Introduction
Competition is defined as a rivalry where two or more parties strive for limited resources.Competition for food and space between different species is relatively common in communities, and the result is almost always a victory for one side and the exclusion of the other.Two species with similar food habits compete for any limited resource, but it is common for species to compete for food resources [1,2].Food competition is one of the important manifestations of food relationships, including intraspecific and interspecific food competition.Intraspecific food competition can be broadly defined as the effort made by two or more individuals of the same species to access a crucial and limited Animals 2024, 14, 849 2 of 17 food resource.In the short term, competing individuals may adopt different behavioral strategies and feeding patterns or establish dominance hierarchies and feeding areas; in the long term, intraspecific food competition can drive ecological role turnover [3].Studying the similarity in the compositions of prey organisms consumed by different organisms (overlapping feeding niches) can help us understand food competition between species; a higher the degree of overlap in food between fish species leads to competition for food; on the contrary, if there is little food overlap between fish species, this means that the greater the differences in prey, the less competition between species.When prey organisms in the environment are limited, fish with similar niches will have severe food competition, which may lead to the turnover of dominant species and changes in the community structure [4].The study of interspecific food relationships among fishes of the same genus in the same sea area can reflect the status of fish utilization of food resources and the existence of competition in the feeding process, and this is the basis for studying the trophic dynamics of food webs [5].Many scholars, both domestically and internationally, have conducted research on fish feeding relationships [6][7][8][9][10], but most of them have focused on relationships between different species or non-congeneric species.The interrelationships among closely related species are particularly unique and significant aspects of species relationships [11].Due to the advancements in theoretical research and scientific technology, there is an increasing focus on the study of food relationships among closely related species [12][13][14].It is anticipated that this attention will lead to a deeper understanding of the mechanisms underlying species coexistence.
Feeding relationship research tools mainly include the traditional stomach content analysis and stable isotope analysis [15,16].A stomach content analysis is a traditional method that has been used for a long time in marine feeding ecology, and its advantage is that it can accurately and intuitively reflect the composition and taxonomic characteristics of the prey consumed by fish [17].However, a traditional stomach content analysis is timeconsuming and labor-intensive, and it has the drawbacks of only being able to characterize transient feeding information as well as not being able to determine the fragmented prey organisms being digested, whereas stable isotopes of carbon and nitrogen in tissues, which can record feeding information on medium to long time scales, have been proven to be the best alternatives and complements to stomach content analysis [18], and they have been widely used in the reconstruction of feeding habits of aquatic organisms.However, the stable isotope technique also has some limitations in that it is unable to visually obtain information on ingested food types and sizes, so a stable isotope analysis and stomach content analysis were combined as complementary techniques to more comprehensively and accurately reflect the food relationships of organisms in the study of feeding ecology.
Pennahia pawak and Pennahia anea have relatively important ecological niches as relatively stable resource fish species in the Beibu Gulf [19].Studies on P. pawak and P. anea have mainly focused on fishery biology [20,21], growth and mortality [22], reproductive biology [23], population structural characteristics and spatial-temporal distribution [24,25], and feeding ecology [26], whereas studies on the interspecies feeding relationship between P. pawak and P. anea have not been reported.These two co-distributed species are not only closely related, but also similar in size, shape, and ecology [19].So, there may be potential feeding competition between them.In this study, using a stomach content analysis and a carbon-nitrogen stable isotope analysis, we compared and analyzed the differences in the feeding habits and trophic niches of P. pawak and P. anea.We included a spatial niche analysis to elucidate the feeding relationship between these two closely related species, aiming to explore the mechanisms of coexistence of sympatric species in the Beibu Gulf.

Sample Collection and Processing
Samples for trophic niche analysis were collected in 2020 from Guangdong Jianghong and Beihai fishing harbors, with 1506 and 1684 samples of P. pawak and P. anea, respectively, which were randomly sampled in four seasons (   season were frozen and sent back to the laboratory for the next step of processing.Further classification and identification of the collected samples will be conducted to determine the correct research species (Figure 1).After thawing the samples in the laboratory, biological indicators such as total length, body length, and weight were measured according to Specifications for oceanographic survey [27], and stomach content samples were retained for freezing and preservation.

Sample Collection and Processing
Samples for trophic niche analysis were collected in 2020 from Guangdong Jianghong and Beihai fishing harbors, with 1506 and 1684 samples of P. pawak and P. anea, respectively, which were randomly sampled in four seasons (Table 1).The samples collected in each season were frozen and sent back to the laboratory for the next step of processing.Further classification and identification of the collected samples will be conducted to determine the correct research species (Figure 1).After thawing the samples in the laboratory, biological indicators such as total length, body length, and weight were measured according to Specifications for oceanographic survey [27], and stomach content samples were retained for freezing and preservation.

Stomach Contents Analysis
After thawing the stomach content samples, the food mass was removed and placed in a clean Petri dish, and the excess water on the surface of the food mass was carefully blotted out with filter paper.Subsequently, efforts were continued to identify the various species of prey within the food mass, aiming to classify these prey organisms down to the finest taxonomic order possible.The food mass was weighed (to the nearest 0.01 g), and the data were recorded.
In this study, the repletion index (RI) and vacuity coefficient (VC) were used to determine the feeding intensity of fish [28].The formulas are as follows: where Wf is the total weight, f, of the fish, and ∑  is the total weight, f, of all the food in the stomach of the fish and the vacuity coefficient (VC).The importance of prey organisms was evaluated using the Index of Relative Importance percentage (IRI%) [29] with the following formula:

Stomach Contents Analysis
After thawing the stomach content samples, the food mass was removed and placed in a clean Petri dish, and the excess water on the surface of the food mass was carefully blotted out with filter paper.Subsequently, efforts were continued to identify the various species of prey within the food mass, aiming to classify these prey organisms down to the finest taxonomic order possible.The food mass was weighed (to the nearest 0.01 g), and the data were recorded.
In this study, the repletion index (RI) and vacuity coefficient (VC) were used to determine the feeding intensity of fish [28].The formulas are as follows: VC = no. of empty stomachs total no. of stomachs (2) where W f is the total weight, f, of the fish, and ∑ W i f is the total weight, f, of all the food in the stomach of the fish and the vacuity coefficient (VC).The importance of prey organisms was evaluated using the Index of Relative Importance percentage (IRI%) [29] with the following formula: where N is the number of a given prey as a percentage of the total number of prey; W is the weight of a given prey as a percentage of the total weight of prey; and F is the percentage of frequency of occurrence of each prey.The Shannon-Wiener diversity index (H ′ ) was used to evaluate the width of the trophic niche of the two fish species.The overlap of trophic niches between the two species [30,31] was also calculated by using Pianka's coefficient of ecosystem overlap (O ij ) to evaluate the intensity of competition between the two species.
H ′ i is the width of the trophic niche for fish i, and P ik is the percentage of the number of prey k in the food of fish i(N).
where s is the total number of prey species consumed by the two fish species, and P ik and P jk represent the number percentage (N) of prey k in the prey composition of fishes i and j, respectively.The value of O ij ranges from 0 to 1.The larger the value, the higher the food similarity, that is, the fiercer the food competition.O ij > 0.3 means that the overlap is effective, and O ij > 0.6 has significant overlap.Individual specialization (IS) was used to evaluate the intensity of intraspecific competition between the two species [32], and the Whin-individual component (WIC) to the width of its trophic niche (H ′ ) was calculated.The range was 0-1, and the smaller the ratio, the lower the individual overlap and the higher individual specialization [33].The related calculation formula is as follows: where P i is the percentage of prey consumed by an individual i to the population of prey in that population, and P ik is the percentage of prey consumed by an individual i to the population of prey in that population.

Carbon and Nitrogen Stable Isotope Analysis
A rectangular cut was made below the dorsal fin and above the lateral line of the fish, the fish epidermis was peeled off from the inner side of the skin, and an appropriate amount of muscle sample was taken from the notch in a 2 mL centrifuge tube and cryopreserved for subsequent processing.
The muscle sample was placed in the sample tray of the freeze dryer (Christ, Osterode, Germany, Alpha1-4/2-4LD Plus), and it was freeze-dried for 48 h at −48 • C until the muscle maintained a constant weight.The sample was taken out, and 2 small steel balls were put into each centrifuge tube containing the sample.Then, it was put in the homogenizer (BIOSPEC MiniBeadbeater-16, Biospec, Bangor, PA, USA) and ground for 1 min.The ground muscle powder was embedded and sent to the Isotope Laboratory of the School of Marine Meteorology, Guangdong Ocean University, where it was measured using an EA Isolink Elemental Analyzer (Thermo Fisher scientific, Waltham, MA, USA) and a 253 Plus Isotope Mass Spectrometer (Thermo Fisher scientific, Waltham, MA, USA).
The stable isotope ratios are expressed according to the internationally accepted δvalue [34], which is calculated using the following formula: Animals 2024, 14, 849 5 of 17 where δX denotes the carbon stable isotope ratio δ 13 C or nitrogen stable isotope ratio δ 15 N, R sample denotes the ratio of carbon or nitrogen stable isotope in the sample, and R standard denotes the ratio of carbon or nitrogen stable isotope in the standard.The carbon isotope ratio is 13 C/ 12 C, and the nitrogen stable isotope ratio is 15 N/ 14 N.
The trophic level formula is as follows [35]: where TL denotes the trophic level of the fish; δ 15 N sample denotes the nitrogen stable isotope signature of the sample; δ 15 N 0 is the nitrogen stable isotope signature of the baseline organism; δ 15 N c denotes the nitrogen isotope enrichment (3.4‰); and TL b denotes the trophic level of the baseline organism.According to our team's pre-monitoring study, Amusium pleuronectes are qualified as baseline organisms because they have stable feeding habits, longevity, and ease of sampling, feeding mainly on plankton and organic detritus [36].Therefore, we chose them as the baseline organisms for calculating the trophic level in this study.Their nitrogen isotope value was determined to be 8.89‰, and the trophic level was set at 2. The degree of overlap between SEAc was used to quantify the area of trophic niche overlap between P. pawak and P. anea [35], and the overlap of the two species was expressed as the proportion of overlap in the maximum likelihood estimate of SEAc, which represents the point estimate of trophic niche overlap between the two species, with an overlap close to 0 indicating that the two ellipsoids are separated, and a value close to 1 indicating that the two ellipsoids are completely overlapped.According to Schoener's D index and according to the results of other stable isotope niche area overlap studies [37], B ij is the niche overlap index, ranging from 0 to 1, with larger values indicating higher overlap, and a B ij value greater than 0.3 is regarded as a meaningful overlap, and a value greater than 0.6 is regarded as a significant overlap [38].

Spatial Niche
The summer close season is from 12:00 a.m.May 1 to 12:00 a.m.August 16.The data used for the spatial niche analysis came from surveys of 26 stations in the Beibu Gulf before (April) and after (August) the period of fishing moratorium.The survey vessel had a power of 441 KW, and the net was a bottom trawl with a 20 m wide net opening, a maximum mesh of 5 cm, and a sac mesh of 2 cm.One hour of trawling was conducted at each station.Samples were collected and preserved on ice and brought back to the laboratory for classification and biological determination.The number and weight of each species were recorded and converted to biomass (kg/km 2 ) and abundance (ind/km 2 ) per unit area.Sample sampling and analysis were performed according to Specifications for oceanographic survey [27]. ( where D is the resource density (kg•km −2 ); C is the catch rate (kg•h −1 ); E is the escape rate taking the empirical value (0.5) [24]; and a is the hourly swept area (km 2 ) of the survey vessel, where the swept width is taken as 1/2 of the length of the float line, and the towing speed is taken as the average towing speed of 3.0 kn.
(2) Spatial niche width and overlap Spatial niche width and overlap were estimated based on catch, where spatial niche width was expressed using the Shannon-Wiener index (H ′ ) [39], and the spatial niche overlap index (Q ij ) was expressed using the Pianka index [30] with the following formula: Animals 2024, 14, 849 where H ′ denotes the spatial ecological niche width, and P ij denotes the mass percentage of fish i in the total biomass of the jth station.Q ij is the ecological niche overlap index, which ranges from 0 to 1, with larger values indicating a higher degree of overlap; specifically, a Q ij value greater than 0.3 is regarded as a meaningful overlap, and a value greater than 0.6 is regarded as a significant overlap [38].

Data Processing
Data were processed, analyzed, and graphed using Excel 2016 and SPSS 25.0, and Pearson analysis was used to test for linear correlation between body length and trophic level for two Pennahia species.Trophic niche of δ 13 C and δ 15 N of the two Pennahia species was determined using the SIAR [40] data package in R 3.6.3.And mean centrifugal distance (CD), carbon range (CR), nitrogen range (NR), mean nearest neighbor distance (MNND), standard deviation of nearest neighbor distance (SDNND), standard ellipse corrected area (SEAc), the area of the convex hull (TA), and other ecological indicators were also used todescribe the trophic niche.[35,38].The distribution of resources was mapped using ArcGIS 10.8 software.

Differences in Feeding Habits within Species
The intraspecific feeding habits of P. pawak varied significantly between seasons.In the spring and winter, the stomach contents of P. pawak had the highest proportions of Macrura (77.99% and 63.78%, respectively), and in the summer and autumn, they had the highest proportions of Pisces (47.37% and 34.73%, respectively); Gastropoda were also present only in the spring (Figure 2).The prey overlap coefficients among the length groups of P. pawak were all greater than 0.6 (except for 171-191 mm), and there was significant food overlap.The prey overlap coefficients between the 71-80 mm length group and the 91-110 mm length group, the 111-130 mm length group, and between the 91-110 mm length group and the 111-130 mm length group were close to 1, with an almost complete overlap of prey species (Table 2).Meanwhile, the individual specialization index of the groups of different lengths ranged from 0.06 to 0.22, with an overall increasing trend with the increase in length (Figure 3).3.1.2.Differences in Feeding Habits between Species P. pawak and P. anea have similar feeding habits and some feeding preferences.The following are some of the most common species in the genus (Table 3).Although both species of Pennahia genus take Pisces and Macrura as prey, the proportion of Pisces in the prey of P. anea is obviously higher than that of P. pawak; there are also some differences in specific species of the two species of Pennahia genus, among which P. anea are dominated by Alpheus, Bregmaceros, and Stolephorus, whereas for P. pawak, in addition to Alpheus and The intraspecific feeding habits of P. anea differed markedly between seasons.The highest proportions of Macrura in the stomach contents of P. anea were found in the spring and winter, and the highest proportions of Pisces were found in summer and autumn; in addition, gastropods were present only in the spring (Figure 2).The P. anea prey overlap coefficients were greater than 0.6 between the length groups, except between the 171-190 mm length group and the 111-130 mm length group, the 131-151 mm length group, and the 151-170 mm length group, which were less than 0.6, and there was significant food overlap (Table 2).Among them, the highest prey overlap coefficient was found between the 71-90 mm length group and the 171-190 mm length group, followed by the 131-150 mm length group and 151-170 mm length group.Meanwhile, the fluctuation range of the individual specialization index of different length groups was 0.08-0.26,which showed a fluctuating upward trend with the increase in length (Figure 3).

Differences in Feeding Habits between Species
P. pawak and P. anea have similar feeding habits and some feeding preferences.The following are some of the most common species in the genus (Table 3).Although both species of Pennahia genus take Pisces and Macrura as prey, the proportion of Pisces in the prey of P. anea is obviously higher than that of P. pawak; there are also some differences in specific species of the two species of Pennahia genus, among which P. anea are dominated by Alpheus, Bregmaceros, and Stolephorus, whereas for P. pawak, in addition to Alpheus and Bregmaceros, Metapenaeopsis barbata and Alpheus bisincisus also accounted for certain proportions.The trophic niche widths for P. pawak and P. anea were 2.21 and 2.29, respectively, and the prey overlap coefficient was 0.56, reaching the level of meaningful overlap.

Isotopic Characteristics and Trophic Levels
The δ 13 C and δ 15 N isotopic characteristics of P. pawak and P. anea are different to some extent.The δ 13 C of P. pawak ranged from −18.08 to −15.22‰, with a mean value of −16.57‰, and the δ 15 N ranged from 14.68 to 16.40‰, with a mean value of 15.60‰ (Figure 4).After the Pearson analysis, the length of P. pawak had a significant positive correlation with δ 15 N (Pearson r = 0.65, p < 0.01), but no significant correlation with δ 13 C (p > 0.05); the δ 13 C of P. anea ranged from −18.12 to −15.10‰, with a mean value of −16.58‰, and δ 15 N ranged from 13.65 to 16.63‰, with a mean value of 15.18‰ (Figure 4).After the Pearson analysis, the length of P. anea had a significant positive correlation with δ 15 N (Pearson r = 0.44, p < 0.01), while there was no significant correlation with δ 13 C (p > 0.05).The δ 13 C ranges of P. pawak and P. anea were similar, with that of P. anea being slightly higher, while P. pawak had a larger δ 15 N range than P. anea.
the increase in the length, with a decreasing trend from the 131-150 mm length gro The trophic level of P. anea gradually increased with the body length (Figure 5).trophic span of P. pawak is smaller than that of P. anea.With the increase in length gro the trophic level of P. pawak increased continuously, while that of P. anea decreased w it increased to a certain extent (Figure 5).The one-way ANOVA showed that the trophic level differences between P. pawak and P. anea were significant (p < 0.01).The trophic levels of P. pawak ranged from 3.70 to 4.21 based on the carbon and nitrogen stable isotope ratios, with an average of 3.97 and a trophic span of 0.51.The trophic levels of P. anea ranged from 3.30 to 4.28, with an average of 3.85 and a trophic span of 0.7.The trophic level of P. pawak showed a rising trend with the increase in the length, with a decreasing trend from the 131-150 mm length group.The trophic level of P. anea gradually increased with the body length (Figure 5).The trophic span of P. pawak is smaller than that of P. anea.With the increase in length group, the trophic level of P. pawak increased continuously, while that of P. anea decreased when it increased to a certain extent (Figure 5).

Trophic Niche
From the stable isotope distribution of carbon and nitrogen (Figure 6), it can be seen that the trophic niches of P. pawak and P. anea overlapped to some extent, but the degree of overlap was only 0.29, which was not significant.Based on the carbon and nitrogen stable isotope analysis, we calculated the trophic niche indexes of P. pawak and P. anea.The result showed that (Table 4) the basal food source, food chain length, level of trophic diversity, overall density, range of trophic niches, space of core niches, and total niche area of P. anea were larger than those of P. pawak.

Trophic Niche
From the stable isotope distribution of carbon and nitrogen (Figure 6), it can be seen that the trophic niches of P. pawak and P. anea overlapped to some extent, but the degree of overlap was only 0.29, which was not significant.Based on the carbon and nitrogen stable isotope analysis, we calculated the trophic niche indexes of P. pawak and P. anea.The result showed that (Table 4) the basal food source, food chain length, level of trophic diversity, overall density, range of trophic niches, space of core niches, and total niche area of P. anea were larger than those of P. pawak.

Spatial Niche
The resources of the Beibu Gulf of P. pawak and P. anea have obvious spatial distribution characteristics, and both of them are mainly distributed in the inner waters of the Beibu Gulf, while the P. pawak is also distributed in the middle and the mouth of the Gulf.P. pawak is distributed in the area of 108.46°~109.50°E, 18.00°~21.33°N, while P. anea is distributed in the area of 108.51°~109.50°E, 20.00°~21.33°N (Figure 7).The resource densities of the two species of Pennahia genus range from 0.08 to 109.52 kg/m 2 and from 0.28 to 29.65 kg/m 2 , respectively, with the average resource density of P. pawak (21.19 kg/m 2 ) being greater than that of P. anea (8.74 kg/m 2 ).The spatial niche widths of P. pawak and P.

Spatial Niche
The resources of the Beibu Gulf of P. pawak and P. anea have obvious spatial distribution characteristics, and both of them are mainly distributed in the inner waters of the Beibu Gulf, while the P. pawak is also distributed in the middle and the mouth of the Gulf.P. pawak is distributed in the area of 108.46 • ~109.50 • E, 18.00 • ~21.33 • N, while P. anea is distributed in the area of 108.51 • ~109.50 • E, 20.00 • ~21.33 • N (Figure 7).The resource densities of the two species of Pennahia genus range from 0.08 to 109.52 kg/m 2 and from 0.28 to 29.65 kg/m 2 , respectively, with the average resource density of P. pawak (21.19 kg/m 2 ) being greater than that of P. anea (8.74 kg/m 2 ).The spatial niche widths of P. pawak and P. anea differ before and after the summer close season, with P. anea having a greater niche width than P. pawak before the summer close season, while P. pawak had a greater spatial niche width than P. anea after the summer close season (Table 5).The spatial niche overlap index of the two Pennahia species was 0.20, which was present but not significant.Neither of the two Pennahia species reached a meaningful level of overlap before or after the summer close season(Table 5). of the two Pennahia species reached a meaningful level of overlap before or after the summer close season(Table 5).

Discussion
This study revealed significant food overlap among individuals of different length groups in both P. pawak and P. anea.The optimal foraging theory suggests that individual feeding differences depend on the phenotype (size, sex, or developmental stage) and prey availability of the individual [41].More precisely, although individuals may consume a wide variety of prey, they adopt more specialized feeding habits to avoid intraspecific competition, depending on the state of the individual and the bait resources available in the habitat.At low population densities, individuals compete for the same dominant prey, but as the population densities increase, this prey resource becomes scarce, and individuals begin to feed on other, different prey, suggesting that increasing population densities lead to increased intraspecific competition, which, in turn, increases the degree of feeding specialization [42].Xia et al. [43] used a stomach content analysis and stable isotope analysis to study the feeding habits of Megalobrama terminalis, and the results showed that the decrease in intraspecific competition was mainly due to individual feeding specialization.However, in recent decades, the fishery resources in the Beibu Gulf have had a declining status due to over-exploitation by the fishery [44,45].There is a significant overlap in the length group of the two species, indicating the potential presence of intraspecific competition.The degree of feeding specialization in both P. pawak and P. anea at the body length stage of 71-190 showed a decrease with the increase in the body length, indicating that there was a possible feeding shift at this stage.

Discussion
This study revealed significant food overlap among individuals of different length groups in both P. pawak and P. anea.The optimal foraging theory suggests that individual feeding differences depend on the phenotype (size, sex, or developmental stage) and prey availability of the individual [41].More precisely, although individuals may consume a wide variety of prey, they adopt more specialized feeding habits to avoid intraspecific competition, depending on the state of the individual and the bait resources available in the habitat.At low population densities, individuals compete for the same dominant prey, but as the population densities increase, this prey resource becomes scarce, and individuals begin to feed on other, different prey, suggesting that increasing population densities lead to increased intraspecific competition, which, in turn, increases the degree of feeding specialization [42].Xia et al. [43] used a stomach content analysis and stable isotope analysis to study the feeding habits of Megalobrama terminalis, and the results showed that the decrease in intraspecific competition was mainly due to individual feeding specialization.However, in recent decades, the fishery resources in the Beibu Gulf have had a declining status due to over-exploitation by the fishery [44,45].There is a significant overlap in the length group of the two species, indicating the potential presence of intraspecific competition.The degree of feeding specialization in both P. pawak and P. anea at the body length stage of 71-190 showed a decrease with the increase in the body length, indicating that there was a possible feeding shift at this stage.
The δ 13 C value is less variable during food chain transmission and can indicate its food source, while the δ 15 N value is relatively enriched in the organism and is generally used to determine the trophic level of the study object [35].It has been shown that the carbon and nitrogen stable isotope compositions of fish are not constant at different growth stages [43].In this study, the carbon stable isotope values of the two species did not vary much with the body length, which may be related to the small enrichment of carbon stable isotopes in the living body on the one hand [46], and on the other hand, although there was a shift in feeding, the main feeding taxa, especially the final trophic sources (elements), were relatively stable.In addition, the nitrogen stable isotopes and trophic levels showed overall increasing trends with the body length, which may be related to the increased feeding capacity of individuals.This pattern of change has been observed in many fish species [47,48], with the growth of individuals, the feeding and digestive organs of fish continuing to improve, the swimming ability continuing to increase, the range of prey selection becoming wider, and the prey of fish shifting from a low trophic level to a higher trophic level.Differences in prey resources across various sea regions can lead to variations in the trophic levels of fish.In a 2009 study on the predominant fish species in the Leizhou Bay area, Lu et al. observed that P. pawak exhibited a trophic level of 2.8 [49].It is evident that the trophic level of P. pawak in Leizhou Bay is lower compared to that in the Beibu Gulf.Furthermore, this discrepancy may be associated with the size range of the collected samples.This study encompassed a diverse size range in the collected samples, potentially contributing to an overall higher trophic level for P. pawak.
The trophic niche width represents the strength of an organism's utilization of habitat and resources, as well as its competitive ability [50,51].In situations of limited resource supply, species with broader niches may exhibit stronger competitive abilities [52][53][54].Strong interspecific competition causes species to reduce feeding on the same prey, thereby mitigating interactions with other species [55].Niche overlap represents the frequency of encounters between two species on the same spatial resource, which can also reflect the potential competitive relationship between species [56].In this study, the niche overlap indexes (O ij = 0.56 and B ij = 0.29) for both P. pawak and P. anea did not reach a significant level.Although niche overlap indicates potential interspecific competition, its intensity also depends on consumer and resource abundance [57].In the present study, it was found that although Macrura organisms were the main prey for two Pennahia species, P. anea predominantly consumed Alpheus, while P. pawak consumed a percentage of Metapenaeopsis barbata and Alpheus bisincisus in addition to Alpheus.This suggests that P. pawak and P. anea, when they have the same prey base, show different feeding preferences by enhancing their feeding on different types of prey.This is consistent with the Evolutionarily Stable Strategy, in which closely related species inhabiting the same sea area will coexist by regulating their own feeding habits to attenuate food competition between species, thus maximizing equilibrium [13].Based on the carbon and nitrogen stable isotope analysis, the difference in δ 13 C variation between P. pawak and P. anea was not significant, but there was no significant niche overlap between the two species, and the total niche width (TA) and the core niche (SEAc) of P. anea were larger than those of P. pawak.Combined with the analysis of food composition, although the food resources of the two species are generally similar, a differentiation in their trophic niches arises due to differences in their dominant feeding prey.
The width of the spatial niche reflects the spatial distribution of the species and its ability to utilize spatial resources [58].The spatial niche width of P. pawak increased while that of P. anea decreased after the summer close season.P. pawak spawns in the April-August period [25], and the summer close season (May 1) facilitates the protection of the spawning population and juveniles, thus increasing its resource density.P. anea spawns in the March-June period [23], two months earlier than the summer close season.This is not enough to protect the spawning population and juveniles is insufficient, so its resource density is in a decreasing trend.Differences in the width of spatial niche is one of the important conditions for species to realize coexistence [58], and this change has a certain mitigating effect on the pressure of spatial resource competition between the two.According to the theory of resource competition [59], in an environment with limited resources, complete competitors cannot coexist, indicating niche differentiation along a certain resource dimension.The overlap between P. pawak and P. anea did not reach a significant level.The differences in their spatiotemporal distribution resulted in a low degree of overlap, reflecting a distinct niche differentiation between the two species.
The classical ecological theory suggests that the niche occupied by a species is limited by a variety of ecological factors, among which biological factors include its own population size, feeding and being fed on, and food competition, and that an increase in population size leads to increased competition [60].When similar organisms with similar diets coexist in a resource-limited environment, dominant species eat whatever they want, while weaker competitors may be forced to eat different items.In marine ecosystems, marine organisms inhabiting the same sea area mainly reduce interspecific competition by realizing the differentiation of trophic, spatial, and temporal niches, among which the differentiation of trophic and spatial ecological niches is particularly important [57].As demonstrated in the present study, in a resource-limited environment, coexisting species of Pennahia may potentially utilize the resources of the Beibu Gulf through the differentiation of their trophic and spatial ecological niches, and they may collectively utilize resources in the Beibu Gulf and achieve coexistence.In addition, the trophic and spatial niches of P. pawak are more advantageous than those of P. anea, suggesting a potentially stronger resource occupancy capability of P. pawak.This has important reference value for fishery management.

Conclusions
This study explores the feeding habits and ecological niches of P. pawak and P. anea in the Beibu Gulf.It highlights how individual feeding differences relate to prey availability, leading to specialized feeding habits to mitigate intraspecific competition.As population densities increase, competition intensifies, prompting individuals to shift their dominant prey.A stable isotope analysis revealed changes in trophic levels with the body length and differences in prey selection across growth stages.Despite having similarities in food resources, the two species exhibit distinct trophic niches due to differences in dominant prey.The spatial niche width also varies between the species, influenced by spawning seasons and resource availability.Niche differentiation aids in reducing interspecific competition, facilitating coexistence.These findings underscore the importance of trophic and spatial niche differentiation in marine ecosystems for species coexistence and have implications for fishery management.Funding: This research was funded by the National Natural Science Foundation of China (U20A2087) and the start-up funds of Guangdong Ocean University (no.060302022104).
Institutional Review Board Statement: Not applicable.All specimens were collected from commercial fishery catches, and data were obtained from dead fishes that were not endangered or protected.The experimental protocols concerning the handling of fishes were in accordance with the requirements of the Institutional Animal Care and Use Committee of Guangdong Ocean University, China.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data presented in this study are available on request Konglan Luo (email: luokonlan@gmail.com).The data are not publicly available due to [the need to reserve the data for future research use].

Figure 2 .
Figure 2. Diet composition by season for Pennahia pawak (A) and Pennahia anea (B) in Beibu Gulf.Circles are in order from inside out: spring, summer, fall, and winter.Eight kinds of prey are displayed with different colors, respectively.

Figure 2 .
Figure 2. Diet composition by season for Pennahia pawak (A) and Pennahia anea (B) in Beibu Gulf.Circles are in order from inside out: spring, summer, fall, and winter.Eight kinds of prey are displayed with different colors, respectively.

Figure 2 .
Figure 2. Diet composition by season for Pennahia pawak (A) and Pennahia anea (B) in Beibu Gulf.Circles are in order from inside out: spring, summer, fall, and winter.Eight kinds of prey are displayed with different colors, respectively.

Animals 2024 , 18 Figure 6 .Table 4 .
Figure 6.The trophic niches of two Pennahia species in the Beibu Gulf.The red circles represent P. pawak, and the blue triangles represent P. anea; the ellipses represent the core niche area, and the convex polygons represent the total niche area.Table 4. The trophic niche indicators of two Pennahia species in the Beibu Gulf.The CD (centrifugal distance), CR (carbon range), NR (nitrogen range), MNND (mean nearest neighbor distance), SDNND (standard deviation of nearest neighbor distance), SEAc (standard ellipse corrected area), TA (the area of the convex hull)values are shown to provide insight into the trophic niches.An overlap index of > 0.3 means that the overlap is effective, and an overlap index of > 0.6 represents significant overlap.

Figure 6 .
Figure 6.The trophic niches of two Pennahia species in the Beibu Gulf.The red circles represent P. pawak, and the blue triangles represent P. anea; the ellipses represent the core niche area, and the convex polygons represent the total niche area.

Table 4 .
The trophic niche indicators of two Pennahia species in the Beibu Gulf.The CD (centrifugal distance), CR (carbon range), NR (nitrogen range), MNND (mean nearest neighbor distance), SDNND (standard deviation of nearest neighbor distance), SEAc (standard ellipse corrected area), TA (the area of the convex hull)values are shown to provide insight into the trophic niches.An overlap index of >0.3 means that the overlap is effective, and an overlap index of >0.6 represents significant overlap.

Table
).The samples collected in each

Table 1 .
Sampling information of two Pennahia species in Beibu Gulf in 2020.

Table 1 .
Sampling information of two Pennahia species in Beibu Gulf in 2020.

Table 2 .
The diet overlap coefficient of Pennahia pawak and Pennahia anea in each length group.Values >0.3 mean that the overlap is effective, and values > 0.6 represent significant overlap.

Table 2 .
The diet overlap coefficient of Pennahia pawak and Pennahia anea in each length group.Values > 0.3 mean that the overlap is effective, and values > 0.6 represent significant overlap. P.

Table 3 .
The diet compositions of the Pennahia pawak and Pennahia anea captured in the Beibu Gulf.The F% (frequency of occurrence), W% (percentage in biomass), N% (percentage of number), IRI (index of relative importance), and IRI% (index of relative importance expressed as a percentage) values are shown to provide insight into each prey item."+" indicates that the ratio was smaller than 0.01%, and "-" indicates absence.

Table 5 .
Spatial niche width and overlap index of two Pennahia species in Beibu Gulf; SCS indicates summer close season.Overlap index of >0.3 means overlap is effective, and overlap index of >0.6 represents significant overlap.

Table 5 .
Spatial niche width and overlap index of two Pennahia species in Beibu Gulf; SCS indicates summer close season.Overlap index of > 0.3 means overlap is effective, and overlap index of > 0.6 represents significant overlap.