Seasonal feeding plasticity can facilitate coexistence of dominant omnivores in Neotropical streams

Coexistence of ecomorphologically similar species in diverse Neotropical ecosystems has been a focus of long-term debate among ecologists and evolutionary biologists. Such coexistence can be promoted by trophic plasticity and seasonal changes in omnivorous feeding. We combined stomach content and stable isotope analyses to determine how seasonal variation in resource availability influences the consumption and assimilation of resources by two syntopic fish species, Psalidodon aff. gymnodontus and P. bifasciatus, in the Lower Iguaçu basin. We also tested the impact of seasonality on trophic niche breadth and diet overlap of these two dominant omnivores. Seasonal changes in resource availability strongly influenced the consumption and assimilation of resources by the two fish species. Both species exhibited high levels of omnivory, characterized by high diversity of allochthonous resources in the wet season. Terrestrial invertebrates were the main component of diet during this season. However, in the dry season, both species reduced their isotopic niches, indicating diet specialization. High diet overlap was observed in both seasons, but the isotopic niche overlap was smaller in the dry season. Substantial reduction in the isotopic niche of P. bifascistus and a shift toward aquatic invertebrates can facilitate coexistence during this season of resource shortage. Feeding plasticity allows omnivorous fish to adjust their trophic niches according to seasonality, promoting the exploitation of different resources during periods of greater resource diversity. This seasonal variation could be an important mechanism that contributes to the resource partitioning and coexistence of dominant omnivores in Neotropical streams.


Introduction
The ecological mechanisms that enable the coexistence of similar species in diverse Neotropical communities are not fully understood, but such coexistence can be facilitated by high trophic plasticity, which can be favored by subtle morphological differences, and changes in omnivorous feeding during periods of food shortage (Chase and Leibold 2003). Omnivory is a widespread feeding strategy where consumers feed at multiple trophic levels (Pimm and Lawton 1978;Pimm 1982;Pimm et al. 1991;Kratina et al. 2012). Although it generally refers to feeding on both animals and plants, omnivory also includes cases where an intermediate consumer and its predator compete for shared resources (Polis et al. 1989;Holt and Polis 1997;Ingram et al. 2012). However, the effects of temporal changes in the diet of sympatric omnivores that utilize the same resources remain unclear, especially in Neotropical streams. Variation in diet breath could be influenced by seasonal changes in precipitation, temperature and nutrient concentrations, which all govern the availability and diversity of food resources (Lowe-McConnell 1999;Abujanra et al. 2009;Ruiz-Cooley et al. 2017). A shift in trophic niche away from possible competitors during periods of low availability of resources may facilitate coexistence of omnivorous species, and reduce the risk of local extinctions Pereira et al. 2010;Wooton 2017).
Omnivorous fish are common consumers in freshwater ecosystems (Barros et al. 2017;Garcia et al. 2016;Costa-Pereira et al. 2017;Silva et al. 2017), but ecologists still seek to understand how seasonality influences their trophic niches, and, consequently, the diversity of ecomorphologically similar species. Competition theory postulates that intraspecific competition during seasons with a shortage of food is dampened by the narrowing of trophic niches and a reduction in diet overlap among species (Pianka 1974;Schoner 1974). On the other hand, omnivores can expand their diet in seasons with high resource availability, resulting in increased overlap of trophic niches among species. However, empirical findings indicate that not all syntopic, congeneric species can expand their diet breath in periods of greater diversity and abundance of resources (Correa and Winemiller 2014;Silva et al. 2017;Neves et al. 2018a, b). This may be related to the traits and food preferences of species, and thus be of further assistance to understanding species coexistence (Chesson 2000;Walker et al. 2013;Alves et al. 2020).
Omnivorous fish in Neotropical streams experience a season with high resource abundance, followed by a season of resources shortages, corresponding to the wet and dry seasons respectively (Pujarra et al. 2017;Quirino et al. 2017). Continuous precipitation during the wet season increases water flow, hindering the establishment of the benthic invertebrate community (Pujarra et al. 2017) and increasing inputs of terrestrial resources (Lisboa et al. 2015;Tonin et al. 2017). The allochthonous subsidies are substantially reduced during the dry season, when aquatic consumers rely predominantly on autochthonous resources (Soares et al. 2013;Correa and Winemiller 2014). It has been hypothesized that these strong seasonal differences in resource availability and composition are responsible for the dominance of omnivorous fish species in the tropics (González-Bergonzoni et al. 2012;Silva et al. 2017). Omnivorous feeding in the tropics may be further promoted by a high abundance of primary producers and high environmental temperatures (González-Bergonzoni et al. 2012).
Omnivory in Neotropical streams is mainly represented by small fish species from the family Characidae. Among these, Psalidodon and Astyanax are highly diverse genera with a wide geographical distribution (Eschmeyer et al. 2016;Rossini et al. 2016;Bonato et al. 2017;Silva et al. 2017;Terán et al. 2020). Although species of these genera are considered to be omnivores (Baumgartner et al. 2012), responses in their feeding and resource assimilation to seasonal changes are poorly understood. Recent works have suggested that seasonal variation in precipitation alters the availability of resources in floodplains, lakes and rivers, which, in turn, promotes changes in trophic niches of characid fish (Quirino et al. 2015;Costa-Pereira et al. 2017;Silva et al. 2017;Neves et al. 2018a, b). Thus, a greater mechanistic understanding of how seasonal changes in resource use affect the trophic ecology of these fish could provide new insights about the coexistence of diverse communities of Neotropical consumers.
In the Iguaçu river basin (Brazil), there are 12 species of the genera Astyanax and Psalidodon (Haluch and Abilhoa 2005;Abilhoa and Duboc 2007;Pavanelli and Oliveira 2009;Baumgartner et al. 2012;Ingenito and Duboc 2014); with Psalidodon bifasciatus (Garavello and Sampaio 2010) and Psalidodon gymnodontus (Eigenmann 1911) among the most common and dominant species (Baumgartner et al. 2012;Neves et al. 2015;Larentis et al. 2016;Delariva et al. 2018;Baldasso et al. 2019). These species were recently relocated to the genus Psalidodon by Terán et al. (2020), previously being recognized as Astyanax gymnodontus and Astyanax bifasciatus. This basin has a high number of endemic tetra fish species but the overall diversity of fish fauna is low in comparison with other basins (Zawadzki et al. 1999;Baumgartner et al. 2012). This region is characterized by distinct dry and wet seasons, and represents a suitable model system to test hypotheses about the trophic ecology and coexistence of congeneric omnivorous fish.
Herein, we investigated seasonal variation in the trophic ecology of the dominant syntopic omnivores, P. aff. gymnodontus and P. bifasciatus, in a stream of the Lower Iguaçu river basin. We combined stomach content and stable isotope analyses (Nielsen et al. 2018) to test how seasonal changes in resource availability and diversity alter: (1) the consumption and assimilation of resources; (2) the trophic position; (3) the trophic niche breadths and niche overlap of both species. We hypothesized a shift in resource consumption and assimilation from mainly aquatic resources in the dry season, to mainly terrestrial resources during the wet season. Based on the predictions of competition theory (Pianka 1974;Schoner 1974), we expected wider trophic niches, higher niche overlap and a greater degree of omnivory in the wet season due to the high input of abundant and diverse allochthonous resources. In contrast, we expected reduced trophic niches in the dry season characterized by reduced resource diversity, and dominance of autochthonous resources. We also hypothesized that both species would feed at higher trophic positions during the dry season, due to a diet comprising mainly animal resources. Finally, we expected both species to share food resources during the period of greatest resource diversity (wet season) and to focus their diet on preferred resources during the period of resource scarcity (dry season). Such information would improve our understanding of the mechanisms that facilitate the coexistence of ecomorphologically similar omnivores in diverse Neotropical ecosystems.

Study area
The Lower Iguaçu river basin is an important hotspot for biodiversity (Abell et al. 2008). In this basin, the Iguaçu National Park protects one of the few remaining areas of Atlantic Rain Forest. About 70% of the fish fauna is considered endemic, comprising small and medium sized species, mainly belonging to the genera Astyanax and Psalidodon (Baumgartner et al. 2012;Terán et al. 2020). Despite its importance, the Lower Iguaçu river basin is threatened by intensive soybean and corn agriculture.
We studied seasonal differences in the trophic ecology of two dominant syntopic species of the genus Psalidodon in the Arquimedes stream (25°9 0 10.25''S 53°16 0 41.86''W). This third-order headwater stream (* 700 m altitude) is approximately 4 m wide and between 0.20 and 0.8 m deep. About 66% of the stream basin area is covered with native vegetation, and surface geology is composed mainly of sand, pebbles and gravel (Neves et al. 2018a, b). Basal resources fueling the stream food webs are dominated by native terrestrial vegetation and aquatic periphyton associated with stones. The climate of this region is humid subtropical (Köppen 1936;Baumgartner et al. 2012), and characterized by cold and dry winters, followed by hot and humid summers (Baumgartner et al. 2012; Fig. S1). Whereas increased precipitation in the wet season leads to abrupt changes in stream water flow, prolonged periods of flooding are not commonly observed like in floodplains. During this period there is also reduced availability of autochthonous resources, because the high flows associated with precipitation hinder the establishment of the benthic communities in these headwater streams (Pujarra et al. 2017).

Field sampling
We performed two field sampling campaigns during the middle of the dry (July) and the wet (December) seasons of 2017 (Fig. S1). We collected fish from 50 m stream reaches using depletion electrofishing with three passes. We highlight that in both seasons, sampling was standardized with the same effort to capture. During this sampling period, P. aff. gymnodontus and P. bifasciatus were the only species of Characidae and the most abundant fish species detected (68% of the total number of individuals). After identification, all specimens were measured and weighed. Before processing for analysis of stomach contents, all specimens were anesthetized in eugenol, fixed in 10% formaldehyde and preserved in 70% alcohol. For stable isotope analysis, adult specimens of each species were immediately euthanized via spinal section and stored on ice for further processing in the laboratory. During the two fish sampling campaigns, we also collected aquatic and terrestrial basal resources and invertebrate prey in the same 50 m stream reaches where fish were collected. Basal resources included terrestrial leaves and seeds from live plants in the riparian vegetation, aquatic biofilm and sedimentary organic matter (SOM). We stored vegetation resources in plastic bags. We sampled biofilm by scraping and washing rocks with distilled water. We also manually sampled SOM from sediment deposits in the three random sections of the streams by gently removing sediment deposits and stored in plastic bottles. We stored all samples of SOM and biofilm in plastic bottles. We sampled aquatic invertebrates (Trichoptera, Hemiptera, Coleoptera, Megaloptera, Odonata, and Aeglidae) manually, using a kick net method. We also used entomological nets to sample terrestrial invertebrates (Coleoptera, Odonata, Hymenoptera, Orthoptera, Diptera, and Oligochaeta) in riparian vegetation around the edges of the stream. We stored both aquatic and terrestrial invertebrate samples in Eppendorf tubes. We froze all fish samples, basal resources and prey individuals at -20°C until processing.
To determine the seasonal differences in the availability of aquatic macroinvertebrate prey, in the two fish sampling campaigns, we randomly collected Surber samples (0.04 m 2 and 200 lm of mesh size, five minutes' duration) from three different substrates across the stream reaches (50 m) concomitantly with the fish sampling (in July and December 2017). We rinsed the samples with running water in a set of sieves of different mesh sizes (2.0, 1.0 and 0.2 mm) and preserved the macroinvertebrate samples in 70% ethanol for later enumeration and identification under a stereoscopic microscope (Opton TIM-2B WF10X). We estimated the potential prey availability as the number of macroinvertebrate individuals found per square meter of stream (density). For this, we divided the number of collected individuals by the Surber sampler area (0.04 m 2 ).

Laboratory procedures
We inspected the gonads of all sampled fish according to the methodology proposed by Vazzoler (1996) and we selected only adult specimens with mature gonads for further analyses. For stomach content analysis, we used a total of 138 fish specimens of similar standard lengths (SL) for both species. In the dry season, we analyzed 60 specimens of P. aff. gymnodontus (mean and standard deviation SL: 8.73 ± 1.24 cm), and 32 specimens of P. bifasciatus (SL: 7.80 ± 1.24 cm). In the wet season, we used 36 specimens of P. aff. gymnodontus (SL: 7.96 ± 1.33 cm), and 15 specimens of P. bifasciatus (SL: 6.56 ± 1.33 cm). We removed fish stomachs and used optical and stereoscopic microscopes (Opton TIM-2B WF10X) to identify the stomach contents to the lowest taxonomic resolution possible, using specific literature for algae (Bicudo and Bicudo 1970) and invertebrates (Mugnai et al. 2010). We quantified food items using the volumetric method (Hyslop 1980), with graduated test tubes and a glass counting plates (Hellawell and Abel 1971). Detritus in the stomachs was considered as sedimentary organic matter (SOM) in different stages of decomposition, and mineral particles.
For stable isotope analysis, we analyzed a total of 40 fish specimens. The stomach contents of these specimens were also analyzed. In each season, we selected ten adult specimens of P. aff. gymnodontus (Dry:10.54 ± 0.69 cm; Wet: 8.99 ± 1.46 cm) and ten adult specimens of P. bifasciatus (Dry: 8.27 ± 1.34 cm; Wet: 7.50 ± 1.03 cm). According to the Cucherousset and Villéger (2015), this sample number is enough to get a robust estimate of the isotopic niches of the species. We dissected samples of dorsal muscle tissue, washed them with distilled water, then lyophilized and homogenized them using a mortar and pestle. We also washed all basal resource and prey items with distilled water. We then lyophilized and homogenized aquatic and terrestrial invertebrates (separated by taxonomic groups), and basal resources. The small insects, including Ephemeroptera, Chironomidae, Coleoptera and Hymnoptera, were macerated as whole individuals. For shrimps and crabs (Aegla sp.), we extracted, lyophilized and macerated muscle tissues. We then stored all samples in 2 mL Eppendorf tubes, before weighing them into tin capsules (1.6 ± 0.2 mg of dry animal tissue and 3.6 ± 4.2 mg for the basal resources samples). We analyzed the samples for nitrogen ( 15 N/ 14 N) and carbon ( 13 C/ 12 C) stable isotope ratios at the Center for Nuclear Energy in Agriculture (CENA) at the University of São Paulo, Brazil. To determine the stable isotope ratios, we used a mass spectrometer system in the continuous-flow (CF-IRMS) mode with a Carlo Erba elemental analyzer (CHN 1110), coupled to a Delta Plus mass spectrometer (Thermo Scientific). We presented all values using the delta notation, which represents the deviation of stable isotope ratios ( 13 C: 12 C and 15 N: 14 N) from universal standards: PDB limestone for carbon and atmospheric nitrogen for nitrogen. We did not correct fish muscle d 13 C for lipids because C:N ratios were below 3.5, indicating negligible lipid content in the samples (Hoffman et al. 2015).

Stomach content analysis
In order to test if diet differed in composition between the two species and two seasons, we applied a two-way PERMANOVA to a matrix of food items of individual fish (Table 1) using the Bray-Curtis index with 9999 permutations (Anderson 2001). To compare seasonal differences in intra-and inter-specific dietary niche breadth, we applied the permutational analysis of multivariate dispersions (PERMDISP, Anderson, 2006), where niche breadth was measured through the dispersion of the diet in space. A priori, we calculated the distance of a group from the centroid through principal coordinate analysis (PCoA), using the dissimilarity measure of Bray-Curtis and the volume matrix of food items (Table 1). The two groups represented the species and the seasons. To test the null hypothesis that niche breadth did not differ among the groups, we calculated F-statistics to compare the mean distance of each sample with the median of the group. We determined the statistical significance from 9999 permutations of the residuals of least squares (Anderson 2006). We then performed post hoc pairwise comparisons using Tukey's Honest Significant Difference test. All statistical analyses were performed in R, version 3.5.2 (R Core Team 2019), using the vegan package (Oksanen et al. 2019).
To determine the dietary overlap between species in each season, we calculated the food overlap index (Pianka 1973) based on the volume matrix of food items (Table 1) using EcoSimR package version 0.1.0 (Gotelli et al. 2015). Following Grossman (1986), we defined overlap as low (0-0.39), intermediate (0.4-0.6) or high (0.61-1). We then used the null model with a RA3 algorithm to evaluate the The asterisk indicates a contribution less than 0.1% significance of Pianka's index (Winemiller and Pianka 1990). We performed 9999 Monte Carlo randomizations that represented a null model, that we then compared with the data (mean niche overlap values for all group pairs). Whereas mean overlap that is significantly lower than expected by chance indicates food partitioning, values higher than those expected by chance indicates food sharing (Albrecht and Gotelli 2001;Bonato et al. 2018). The observed overlap was compared with the null model at a significance level a = 0.05 to evaluate if the observed pattern differed from the randomly simulated overlap values.

Stable isotope analysis
Turnover rate of animal tissues is an important consideration when interpreting temporal variation in dietary assimilation (Kaymak et al. 2015;Garcia et al. 2016;Oliveira et al. 2017 , we assumed that the isotopic ratios of adult fish collected in July reflected the assimilation of resources consumed since April (Fig. S1) and, therefore, represented the diet during the wet season. Similarly, the isotopic ratios of adult fish collected in December reflected the assimilation of food resources consumed since July (Fig. S1), representing the diet during the dry season. We tested whether carbon (d 13 C) and nitrogen (d 15 N) stable isotope ratios differed between species, and between wet and dry seasons using two-way Analysis of Variance (ANOVA) in the car package (Fox and Weisberg 2019). Also, we used a one-way ANOVA to test for differences in density of aquatic macroinvertebrates between seasons. We assured that the assumptions of normality (Shapiro-Wilk test) and homoscedasticity (Levene's test) were met.
To estimate the trophic position of both fish species in each season, we applied Bayesian trophic positions measured in tRophicposition package (Quezada-Romegialli et al. 2019), considering 20,000 iterations and five Markov chains. As the model implementation considered a maximum of two different baselines, we pooled (1) biofilm, SOM, and terrestrial vegetation into a basal resource baseline and (2) aquatic and terrestrial invertebrates into a prey baseline. We then calculated pairwise comparisons among species and seasons to determine any differences in estimated trophic position.
Fractionation is a relevant point to estimate the trophic position and resource assimilation. Here, for fish, we used Trophic Discrimination Factors (TDF) of 1.3 ± 0.3% for C, and 2.9 ± 0.32% for N (McCutchan et al. 2003). This enrichment is appropriate for muscle tissues of omnivorous fish that consume mixtures of plant and animal diet (McCutchan et al. 2003). We also considered specific TDF values for plant resources (Bastos et al. 2017). Here, we followed the method of Neres-Lima et al. (2016) and doubled the mean discrimination factor and the variability estimate (SD) by propagation of error (H (2.SD 2 ), yielding the values 2.6 ± 0.42% for C, and 5.8 ± 0.45% for N.
To estimate the degree of omnivory and isotopic niche overlap between the two species in each season, we calculated the isotopic metrics proposed by Cucherousset and Villéger (2015). In order to quantify isotopic niche size for the two species, we used convex hull area to determine isotopic similarity and isotopic nestedness (Cucherousset and Villéger 2015). Isotopic similarity (ISim) is the ratio between the isotopic overlap of two species and the total volume that both species occupy in the stable isotope space (Villéger et al. 2011). However, this index can be influenced by differences in the size of the convex hull of the two groups of organisms. Therefore, isotopic nestedness (Ines) was used as a complementary measure, representing the ratio between the volume of the intersection and the minimal volume filled by a species. Ines ranges from 0 with no isotopic overlap to 1 when the group with the lowest isotopic richness fills a subset of the isotopic space filled by the group with the highest isotopic richness (Cucherousset and Villéger 2015). We inferred the degree of omnivory based on the variation in nitrogen (d 15 N) stable isotope ratios calculated as ISim and Ines, where large variation in d 15 N values indicates assimilation of resources from different trophic levels. This approach quantitatively estimates several aspects of isotopic diversity and provides abundance-weighted and unitless indices, which can be easily compared. These metrics proposed by Cucherousset and Villéger (2015) are mathematically independent of the number of organisms used in the calculation and they can account for the differences in population density or biomass.
To determine the relative contributions of different diet sources assimilated by P. aff. gymnodontus and P. bifasciatus in each season, we applied Bayesian stable isotope mixing models (Moore and Semmens 2008;Parnell et al. 2010) from the MixSIAR package (Stock and Semmens 2016a). We compared the relative contribution of four isotopically homogeneous diet categories (Table S2, Fig. S4): aquatic invertebrates (Odonata nymph, Hemiptera nymph, Trichoptera larvae, Megaloptera larvae, Coleoptera larvae, Aeglidae adult), terrestrial invertebrates (Hymenoptera, Orthoptera, Araneae, Oligochaeta, Coleoptera adult, Diptera adult, Odonata adult), aquatic biofilm and SOM, and terrestrial plants (Pteridophita, Spermatophyta). To prevent erroneously attributing food items that are not consumed by the fish, we used the stomach content data as initial information. Because isotopic composition of terrestrial plants was very different from the composition of fish consumers (Fig. S7), we performed the mixing model analyses with and without terrestrial plants and our final model were composed of aquatic and terrestrial invertebrates and SOM. Both analyses provided qualitatively similar results and showed little contribution of terrestrial plants to the diet of fish consumers (Table S5 and S6). Therefore, we present the analysis without terrestrial plants below. We fit the model with a Markov chain Monte Carlo sampling with number of chains = 3; chain length = 100,000; burn in = 50,000; thin = 50, model 4 (Resid*Process) error structure (Stock and Semmens 2016b), and with ''uninformative/generalist'' prior in MixSIAR. A hierarchical Bayesian mixing model computes the effects of covariates (factors: species and seasons) when estimating food assimilation patterns by consumers (Semmens et al. 2009). Here, we used this approach to evaluate potential changes in the assimilation of basal food sources considering species and seasons (wet and dry) as fixed factors, with species nested with seasons. Together with the source isotope values, we included elemental carbon (%C) and nitrogen (%N) composition data as concentration dependence because the C and N composition of plant and animal resources can differ (Phillips et al. 2014).
We examined model convergence using diagnostic tests (Gelmin-Rubin, Heidelberger-Welch and Geweke) and trace plots.
To verify that the isotopic variability of consumers, after accounting for fractionation, fell within the variability of their assimilated resources (Phillips et al. 2014), we performed isotope mixing polygon simulations using the sp (Pebesma and Bivand 2005;Bivand et al. 2013) and splancs (Rowlingson and Diggle 2017) packages. This method quantitatively determines the boundaries of possible source values in the d 13 C-d 15 N biplots that can explain the isotope mixture (Smith et al. 2013). When individual consumers are located outside the 95% mixing polygon region, they could not be confidently explained by the available food sources (Phillips et al. 2014). We determined that it was not necessary to remove any individual from the final isotope mixing models (Fig. S6).

Results
Stomach content analysis revealed that the main food items in the diets of both fish species were leaves, seeds, insects and detritus (comprising 99.3% of all diet, Fig. 1, Table 1). As predicted, seasonal changes in resource availability and diversity altered the consumption and assimilation of resources, the trophic niche breadths and the niche overlap of both species. There were significant differences in the diet of the two species (Pseud-F = 2.24, p = 0.029), and between the dry and wet season (Pseud-F = 2.77, p = 0.003), with a significant interaction between species and season (Pseud-F = 2.27, p = 0.003). In the dry season, the stomach content of P. aff. gymnodontus consisted mostly of leaves, seeds (56.7%) and terrestrial insects (17%), whereas the stomach content of P. bifasciatus consisted mostly of leaves, seeds (40.4%) and detritus (39.7%). In the wet season, the stomach content of P. aff. gymnodontus consisted mostly of terrestrial insects (39.9%), followed by leaves and seeds (28.7%), whereas the stomach content of P. bifasciatus had a greater proportion of leaves and seeds (57.3%), followed by aquatic (20.5%) and terrestrial (19.5%) insects ( Fig. 1; Table 1).
The two fish species were significantly different in their trophic niche breadth in different seasons (ANOVA, F 3,126 = 4.74, p = 0.004; Fig. 2; Table 2). Psalidodon bifasciatus had a larger diet breadth during the dry season (PERMDISP, distance of centroid, 0.57), whereas P. aff. gymnodontus had a larger diet breath during the wet season (0.60). There was relatively high diet overlap in both dry (Pianka's index = 0.74) and wet (Pianka's index = 0.81) seasons. This high overlap in both seasons resulted from sharing leaves, seeds, detritus, aquatic insects, adult Coleoptera and Hymenoptera as resources ( Fig. S3; Table S1).
There was no difference in d 13 C between species (ANOVA: F 1,36 = 2.28, p = 0.14), but there was a significant difference between the wet and dry seasons (ANOVA: F 1,36 = 8.25, p = 0.007), indicating that both fish species utilized different carbon sources in each season. There was no interaction effect between species and season (ANOVA: F 1,36 = 0.002, p = 0.96). For both species, the d 13 C had slightly Fig. 1 The relative contribution (%) of main food resources consumed by P. aff. gymnodontus (Pag) and P. bifasciatus (Pbi) during the dry and wet seasons. Proportions of resources in diet were estimated by volume from stomach content analysis Fig. 2 Inter-seasonal differences in breadth of diet (measured as distance to centroid) of P. aff. gymnodontus (Pag) and P. bifasciatus (Pbi), assessed by permutational analysis of multivariate dispersions (PERMDISP). Box lower and upper endpoints represent the 25th and 75th quartiles, respectively and the horizontal bar inside each box represents median diet breadth (for p-values, see Table 2). Different letters indicate significant differences (Post-hoc test)  higher values during the dry season (-24.32 to -20.57%) than the wet season (-25.39 to -22.76%: Fig. S5). There was no difference in mean d 15 N between species (ANOVA: F 1,36 = 0.05, p = 0.81) or season (ANOVA: F 1,36 = 0.53, p = 0.47) and no significant interaction between species and season (ANOVA: F 1,36 = 1.59, p = 0.21). Contrary to our expectation, there was no significant difference in mean trophic position of both species between the dry season (P. aff. gymnodontus: 3.34; P. bifasciatus: 3.29) and the wet season (P. aff. gymnodontus: 2.55; P. bifasciatus: 2.75; Fig. 3; Table 2). Both species contracted their isotopic niches during the dry season (cf. wet season), although this effect was more pronounced in P. bifasciatus (Fig. 4). There was higher isotopic similarity (0.36; Fig. 4; Table S4) and isotopic nestedness (0.77) between the two species in the wet season compared with the dry season (ISim = 0.17; Ines = 0.70), suggesting that both species expanded their isotopic niches and shared more similar resources in wet season. Major variation in ISim and Ines indicated a higher degree of omnivory during the wet season than in the dry season.
In agreement with our hypothesis, we observed a shift in resource consumption and assimilation from mainly aquatic resources in the dry season, to mainly terrestrial resources during the wet season. The relative contribution of the main resource types to the diet of both fish species, inferred from the stable isotope mixing model, differed between the dry and wet season ( Fig. 5; Table S6). In the dry season, P. bifasciatus assimilated mainly aquatic (79 ± 10%, mean ± SD) and terrestrial (14 ± 8%) invertebrates, while P. aff. gymnodontus assimilated aquatic (67 ± 13%) and terrestrial (24 ± 12%) invertebrates. In addition to a larger contribution of aquatic invertebrates to the diet of both species, there was a higher density of aquatic macroinvertebrates in the dry season than in wet season (ANOVA: F 1,4 = 12.16, p = 0.025; Fig. S2), suggesting increased availability. In the wet season, when availability of aquatic macroinvertebrates was lower (Fig. S2), both fish species increased consumption and assimilation of terrestrial invertebrates. Psalidodon aff. gymnodontus assimilated the highest proportion of terrestrial invertebrates (68 ± 19% of their total diet), together with a smaller proportion of aquatic invertebrates (17 ± 9%) and SOM (15 ± 17%). Similarly, P. bifasciatus assimilated a great proportion of terrestrial invertebrates (59 ± 19% of their total diet), a low proportion of aquatic invertebrates (27 ± 12%), and a similar proportion of SOM (14 ± 16%) as P. aff. gymnodontus.

Discussion
Seasonal changes strongly influenced the consumption and assimilation of resources by the two dominant fish species, altering the diet breadth and isotopic niche of both species, and the extent of isotopic overlap between them. In dry season, P. bifasciatus and P. aff. gymnodontus assimilated high proportions of aquatic invertebrates. In this season, the high availability of aquatic macroinvertebrates was associated with a narrowing of isotopic niches of both species and reduced isotopic niche overlap (ISim = 0.17; Ines = 0.70), indicating greater resource specialization. During the wet season, when Neotropical aquatic ecosystems are known to have an increased availability of terrestrial resources (Correa and Winemiller 2014;Lisboa et al. 2015;Tonin et al. 2017), both species assimilated high proportions of terrestrial invertebrates. In this season, both fish species amplified their isotopic niches, indicating a greater level of omnivory, and increased their isotopic niche overlap Although the stomach content analysis indicated high resource overlap in both seasons, the low overlap of isotopic niche may suggest species-specific differences in assimilation of consumed resources. These findings demonstrate that seasonal changes in resource availability and composition lead to partially differential resource use by two syntopic omnivores with otherwise widely overlapping trophic niches.
Although both species rely largely on terrestrial invertebrates during the wet season, our findings agree with previous studies from streams (Baldasso et al. 2019), rivers (Pini et al. 2019) and reservoirs (Delariva et al. 2013;Mise et al. 2013), suggesting higher consumption of terrestrial invertebrates by P. gymnodontus. Leaves and seeds were consumed in considerable proportions by both species (Fig. 1), but their assimilation into the fish tissues was very low ( Fig. S7; Table S5). Similar high consumption but low assimilation of terrestrial plants has been observed in other tetra fish species (Bastos et al. 2017;Bonato et al. 2018). At least five factors could explain this apparent discrepancy between diet composition inferred from stomach contents and food assimilation inferred from stable isotope analysis. First, individual food items differ in their digestibility and nutritional quality (protein content; Bowen et al. 1995), and the two methods focus on trophic processes that operate on different time scales (Hyslop 1980;Bastos et al. 2017;Nielsen et al. 2018). Such differences in assimilation efficiency may be related to a lack of specialization in the digestive tract to facilitates the breakdown, digestion and assimilation of plant material (German et al. 2009;Pelster et al. 2015), and to the higher nutritional quality of animal prey (Gerking 1994). Second, the high intake of plant material may result from an opportunistic feeding strategy (Winemiller 1989), where consumers feed on the most available resources (in our case terrestrial plants) during shortages of other resources, even when these resources are not assimilated. Third, some studies suggest that plant material can be accidentally ingested when fish attempt to capture terrestrial insects on the water surface (Bastos et al. 2017;Bonato et al. 2018). However, the high proportion of ingested plant material suggests that accidental ingestion plays a minor role in this case. Finally, ingested plant material may include a biofilm, which is the main resource assimilated into fish tissues. The fifth reason for why fish may ingest but not assimilate plants is that this material is used to change the intestinal microbiota and to improve the intestinal transit and digestion time for more nutritious food (Baldo et al. 2015;Serra et al. 2019). The relative contribution of these factors remains a promising venue for future research. Isotopic overlap between P. aff. gymnodontus (Pag) and P. bifasciatus (Pbi) during the dry and wet seasons, illustrated in d 13 C and d 15 N isotopic space. Isotopic overlap metrics were measured using the isotopic richness of the two species and the volume of isotopic space they shared. Isotopic overlap on each stable isotope axis is illustrated by the overlap of the colored segments symbolizing range of values for each species In agreement with optimal foraging theory (Stephens and Krebs 1986; Schoener 1971), seasonal variation in resource supply may allow exploitation of preferred diet items during periods of high abundance. The greater feeding opportunities in the wet season have been associated with an expansion of trophic niches (Costa- , and here, with an increase in omnivory. During the wet season, P. aff. gymnodontus assimilated mainly terrestrial invertebrates, whereas P. bifasciatus assimilated terrestrial invertebrates, aquatic invertebrates and SOM. Similarly, high contributions of terrestrial invertebrates to the diets of aquatic consumers during periods of greatest rainfall (wet season) have been recorded by studies in temperate environments (Perkins et al. 2008). Both species contracted their isotopic niches during the dry season (cf. wet season), although this effect was more pronounced in P. bifasciatus. This supports the hypothesis that consumers have more specialized diets when their preferred resources are more abundant (Stephens and Krebs 1986;Perry and Pianka 1997). Moreover, our findings agree with observations from Neotropical ecosystems that not all fish species expand their niches in the rainy season or contract them during the dry season (Correa and Winemiller 2014;Silva et al. 2017).

Fig. 5
Density plots depicting the relative isotopic contribution of food resources to P. aff. gymnodontus (Pag) and P. bifasciatus (Pbi) during the dry and wet seasons. The relative diet contributions were estimated using a MixSIAR Bayesian mixing model. SOM denotes sedimentary organic matter Variation in d 15 N across the time reflects differences in the degree of omnivory (Post and Takimoto 2007), the relative assimilation of nutrients (Hopkins and Ferguson 2012), or the protein content of consumed resources (Kelly and Martinez del Rio 2010). Here, the absence of significant differences in d 15 N resulted in similar trophic positions of the two species in both seasons. This indicates that, although the species differ in their carbon sources, both species maintained an omnivorous diet across seasons. Moreover, the similarity in the trophic position of species is likely related to the fact that the fish mainly assimilated animal resources in both seasons. Therefore, seasonal changes in the diet of both species alter their trophic niche breadth more than their trophic positions, which may have fundamental implications for how these species interact and compete for food resources.
Multiple mechanisms are likely to combine to facilitate coexistence in diverse tropical ecosystems (Vergnon et al. 2009). Stabilizing mechanisms act to reduce niche overlap through resource partitioning and spatio-temporal variation in environmental conditions (Chesson 2000). Equalizing mechanisms result from trade-offs that reduce fitness differences between potential competitors (Tilman 1990). Any two species must differ in both their requirements and impact on resource availability in order to locally coexist (Chase and Leibold 2003). In this study, the coexistence of two Psalidodon species was facilitated by seasonal variations in availability of resources and trophic plasticity of these omnivores, allowing them to partition resources. Both intrinsic (feeding preferences, species interactions) and extrinsic (environmental heterogeneity) changes over time and space determine the dynamics and coexistence of species (Chase and Leibold 2003;Hart et al. 2016). The interactions among species reflect trade-offs between their relative competitive abilities and their responses to these environmental variations (Wootton 1998). Similar isotopic composition of food sources in both seasons would indicate that the isotopic niche of consumers was reduced by preferential prey consumption (Table S2). Then, our findings suggest that the high abundance of aquatic insects in dry season ( Fig. S2; Pujarra et al. 2017), might lead both species to specialize in these resources and increase the competitive pressure between them. Specially for P. bifasciatus, a shift toward feeding on aquatic invertebrates and a contraction of trophic niche by may be the mechanism that alleviates this pressure and facilitates coexistence during this season of resource shortage.
In conclusion, seasonal variation in environmental conditions influences trophic niche breadth and resource overlap of syntopic omnivorous fish species. These changes are germane for the detailed understanding of coexistence of phylogenetically related species in diverse Neotropical ecosystems. Trophic plasticity allows them to exploit alternative resources during periods of greater availability. To fully understand this trophic plasticity, future research should investigate how the food preferences of each species alleviate interspecific competition for shared resources. Then, we could infer if such shifts in resource supply can favor better competitors for aquatic resources during the dry season, but better competitors for terrestrial resources during the wet season. Considering future climate forecasts for South America (Chou et al., 2014), which indicate a drier and warmer climate (Chou et al. 2014;FIOCRUZ 2016), these findings highlight the need for a better mechanistic understanding of coexistence of species in changing freshwater ecosystems. otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.