Biodiversity Data Journal 9: e60745 CO) doi: 10.3897/BDJ.9.e60745 open access Research Article Aquatic macroinvertebrates in Uruguayan rice agroecosystem Leticia Baot, Sebastian Martinez$, Monica Cadenazzi!, Monica Urrutia, Lucia Seijas?, Enrique Castiglioni* + Crop Protection Department, University of Republic, Montevideo, Uruguay, Montevideo, Uruguay § Instituto Nacional de Investigacion Agropecuaria (INIA), Programa Nacional de Investigacion en Produccién de Arroz, Laboratorio de Patologia Vegetal, Estacién Experimental INIA Treinta y Tres, Ruta 8 km 281, Treinta y Tres, C.P. 33000, Treinta y Tres, Uruguay | Department of Statistical Biometrics and Computing, University of Republic, Paysandu, Uruguay q University of Republic, Montevideo, Uruguay # Eastern Region University Centre, University of Republic, Rocha, Uruguay Corresponding author: Leticia Bao (baoleticia@gmail.com) Academic editor: Ben Price Received: 13 Nov 2020 | Accepted: 18 Feb 2021 | Published: 04 Mar 2021 Citation: Bao L, Martinez S, Cadenazzi M, Urrutia M, Seijas L, Castiglioni E (2021) Aquatic macroinvertebrates in Uruguayan rice agroecosystem. Biodiversity Data Journal 9: e60745. https://doi.org/10.3897/BDJ.9.e60745 Abstract This work is a first approach to the knowledge of insects and other aquatic macroinvertebrates of rice agroecosystems from eastern Uruguay. The composition of the groups collected may represent an approximation to the knowledge of the quality of water sources associated with Uruguayan rice production. Sampling of aquatic macroinvertebrates was carried out during the grain-filling stage in crops without insecticide use, in three localities of Treinta y Tres Department. In each crop, macroinvertebrates were collected with a Surber-type network at the inlet and outlet of water to and from the paddy field and a neighbouring control area. Differences in morphospecies composition were found according to the location and source of water. Insecta was the most represented class in macroinvertebrate samplings (41.5%). Diptera (59.9%), Hemiptera (16.3%) and Ephemeroptera (14.0%) were the most abundant orders within insects. The Richness and Shannon Diversity Indices were higher than those recorded for similar studies in Costa Rica, Italy and Australia. © Bao L et al. This is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. 2 Bao L etal Keywords biodiversity, bio-indicators, water quality Introduction According to the international Ramsar convention, rice fields are considered as artificial wetlands. Around 100 Ramsar sites in the world include rice areas and perform important ecological functions which sustain remarkable biodiversity (Rizo-Patron et al. 2013). Ecosystemic services provided by rice crops include nutrient recycling, soil generation and biodiversity, with contributions that vary depending on the region (Ramsar Convention Wetlands 2018). Natural water sources harbour more than 6% of all insect species in the world and a number of those can colonise artificial water ponds formed for rice cultivation or paddy fields (Dijkstra et al. 2014). Freshwater animal diversity is represented by approximately 126,000 species, with some underestimation of the true diversity for several important groups, which represent 10% of the total animal species actually recognised (Balian et al. 2008, Dijkstra et al. 2014). Insects represent 60% of all those known animal species reported from freshwater habitats and account for around 100,000 species, followed by vertebrates (14.5%), crustaceans (10%), arachnids (5%) and molluscs (4%) (Balian et al. 2008, Dijkstra et al. 2014). With the reduction in natural wetlands, the occurrence of rice agroecosystems over wide areas of eastern Uruguay could substitute for different groups of animals for the loss of habitat, particularly for aquatic organisms. Therefore, in this context, rice paddies can be evaluated as habitats that can act as a viable refuge for organisms that inhabit surrounding natural habitats (Lupi et al. 2013, Wakhid et al. 2020). The rice crop in Uruguay is produced under irrigation, so the semi-aquatic temporary environment resulting from tilling to pre-harvest could be favourable for insects or other groups that have one or more of their life stages in aquatic environments. Changes in the abundance of aquatic organisms, including insects, have been observed in relation to different rice growth phenological stages (Molozzi et al. 2007). Although marked seasonal variations in abundances of aquatic organisms have been reported, most abundance patterns in different phenological stages of rice crop have been reported as similar amongst different growing seasons (Che Salmah et al. 2017). The short period of flooding in rice fields can affect survival of some species with longer life cycles while favouring some groups (Lupi et al. 2013). When rice is flooded, the water used for irrigation is subjected to rapid chemical and physical changes, such as thermal and nutrients fluctuations. Advancing the growth stages, height and foliar area of rice plants produces thermal fluctuations that affect the creation of a microclimate in shaded areas (Lupi et al. 2013). Additionally, the abundance and diversity of aquatic organisms in rice farms can be modified under variations in food sources, habitat availability, soil drainage and management practices, mostly by machinery use and agrochemical application (Che Aquatic macroinvertebrates in Uruguayan rice agroecosystem 3 Salmah et al. 2017). Agrochemical application is one of the main management practices affecting diversity and abundance of aquatic organisms, not only by the direct toxic effects, but also by changing the physico-chemical conditions of water (Rizo-Patron et al. 2013, Che Salmah et al. 2017). For example, Rizo-Patron et al. (2013) found a greater number of families and species of aquatic invertebrates in organic rice (i.e. no agrochemical use) than in conventional rice cultivation. Thus, it would also be interesting to know the impact on the conservation of the diversity of aquatic organisms in rice fields with a low input of agrochemicals (i.e. no insecticide use and limited use of fungicides), as occurs in some rice-growing areas of Uruguay (Tseng et al. 2020). This particular manner of cultivation (i.e. cultivation under irrigation during part of the crop cycle and limited use of agrochemicals) with respect to other crops and the location in wetland regions presents Uruguayan rice fields as potential biodiversity reservoirs for aquatic organisms (Maltchik et al. 2011, Luo et al. 2014). In addition, depending on the crop management, conservation of several groups of aquatic organisms could be promoted (Settle 1996, Stenert et al. 2009, Gurr et al. 2012). In the same way, natural and semi-natural habitats can provide over-wintering sites, refuges after disturbance and sources of alternative prey for natural enemies in the landscape where crop fields are located (Blitzer et al. 2012). As different groups have different degrees of water pollution tolerance, the composition of macroinvertebrate communities would change amongst water sources and could generate important information about water quality (Rizo-Patron et al. 2013, Melo et al. 2015). Research related to aquatic macroinvertebrates in Uruguay refers to natural environments (Morelli and Verdi 2014), but information related to agroecosystems in the region is scarce (Rizo-Patron et al. 2013, Vilaboa 2012, Gonzalez et al. 2014), while recent research highlights the need for this information to advance the definition of environmental indicators for rice cultivation (Tseng et al. 2020). Furthermore, taking into account the extended area and the production under irrigation, the rice agroecosystems in Uruguay could have the potential to provide ecosystem services and contribute to conservation of regional wetlands biodiversity, whenever insecticide use remains in a low percentage of total productive area (Bao and Martinez 2018). Thus, obtaining information on the current state of the biodiversity of these systems is highly desirable for current decision-making in rice management tending to valuate and preserve the potential mentioned above. This does not replace the need to preserve wetland areas, since beyond the registered biodiversity, differences can be registered between the macroinvertebrate communities associated with paddy fields and wetlands given the dynamics of each of these systems (Pires et al. 2015). This work represents the first approach to the study of aquatic organisms in rice crops and for the evaluation of the role of Uruguayan rice agroecosystems in_ biodiversity conservation of aquatic organisms. The information generated could help to define groups relevant for the ecosystemic services they bring and propose management measures that contribute to their conservation. In addition, group composition in different samples could represent an approximation to water quality knowledge of water sources associated with rice crop production. 4 Bao L et al Material and methods Study area and sampling Sampling was undertaken in February 2015 in three rice crops from Treinta y Tres Department in north-eastern Uruguay (Fig. 2). This Region represents the main area of rice production in Uruguay and belongs to Laguna Merin Basin, where the main production activities are rice and livestock farming. In each locality (site), each crop was sampled in three locations considering water circulation (water source), the points where water entrance (E) to the crop and water outlet (O) occurred from the crop and a control zone (C). Each C zone was a natural wetland consisting of a water body, independent from the crop water source, located at a distance of at least 50 metres from the border of the crop. A LaCharqueada A ElTigre Brazil eae A Julio Maria Sanz Uruguay Treinta yTres A e Figure 1. EESI The location of sampling points for aquatic macroinvertebrates with a Surber net in the rice agroecosystem in Uruguay. Top left corner South America, with Uruguay shaded in grey, left image Uruguayan map indicating localities sampled; right, expanded image of localities sampled. In each crop at each water source (E, O, C), three points were selected and water physicochemical properties were registered: dissolved oxygen (DO), temperature and conductivity (completely randomised design). After this, from each _ point, macroinvertebrates were collected with a Surber net (30 x 30 cm section) for a 1-minute period. Collected material was placed in vials containing 90% alcohol. Collected individuals were counted and organised by taxa with the aid of taxonomic keys, while insect specimens were sorted to morphospecies, the remaining macroinvertebrates were sorted to family (Springer 2006, Springer 2008, Merritt and Cummins 2007, Bentancourt et al. 2009) and functional groups were adjudicated following Ramirez and Gutiérrez-Fonseca (2014). Data analysis Mean values of physicochemical properties for each water source and within each site were calculated and compared by the Tukey Test (p < 0.05). Functional groups and Aquatic macroinvertebrates in Uruguayan rice agroecosystem 5 physicochemical properties were analysed by the Mantel Test to study the correlation between water properties and the composition of functional groups (SAS Institute 2009). Aquatic macroinvertebrate families = Baetidae 100% = Caenidae 90% = Syrphidae Bn = Ceratopogonidae Culicidae 70% = Chironomidae 60% — = Corixidae ® Hydrophillidae = Curculionidae 50% 40% = Noteridae % relative abundance 30% m@ Hydrachnidae 20% = Hydrozetidae* 0% = Acarina unidentified family a Daphniidae Entrance Outlet Control = Cyclopidae Water source ~ Ampullariidae Other aquatics Figure 2. EES] Relative abundances of aquatic macroinvertebrates collected with a Surber net in the rice agroecosystem, grouped by water source. *indicates uncertain identification at this level. Using the insect data counting matrix, richness estimators (Chao 1 and Chao 2, Jackknife 1 and Jackknife 2, Bootstrap and ACE) were calculated and a species accumulation curve was performed with EstimateS 9.1.0 software (Colwell 2013). Species composition was analysed and compared by the ANOSIM Test through Bray-Curtis distances and the SIMPER Test was performed to determine the species that contribute to the major differences between localities and water source (PAST software, Hammer et al. 2001). Using the same software and the Bray-Curtis distances data matrix, a Principal Component Analysis (PCA) with Biplot was undertaken. A heatmap was constructed with the heatmap function using realtive abundances of each taxonomic group (R Core Team 2020). Abundance, Richness and Diversity Indices (Simpson, Shannon, Evenness, Margalef, Equitability and Berger-Parker) were calculated for each sample corresponding to each water source and locality (PAST Software, Hammer et al. 2001) and values were analysed through ANOVA and compared by the Tukey Test (p < 0.05) in Infostat (Di Rienzo et al. 2017). 6 Bao L etal Results Water physicochemical properties varied within each site and also between each water source (Table 1). Dissolved oxygen (DO) values presented statistical differences between sites (F = 5.64, p = 0.011) and between water sources (F = 4.94, p = 0.017). Conductivity values were different only between sites (F = 123.1, p = 0.0001), with higher values in La Charqueada. Table 1. Physicochemical properties of water at different localities and in different water sources (entrance: E, outlet: O and Control zone: C) in rice agroecosystem in Treinta y Tres, Uruguay. Locality (site) Water source Temperature (°C) DO (O2mg/I) Conductivity (yS) La Charqueada (33°11.824'S 53°50.301'W) Entrance 26.57 +0.64a 4.2040.89abc 163.67+10.02ab Outlet 26.47+0.64a 6.2740.89c 139.304+10.02a Control 30.10+0.64b 6.16+0.89c 180.50412.27b El Tigre (33°13.625'S 54°0.387'W) Entrance 27.53+0.57a 9.71+0.89b 89.7343.79a Outlet 29.63+0.57b 9.6340.89b 94.73+3.79a Control 28.23+40.57ab 1.95+0.89a 102.63+3.79a Julio M. Sanz (33°11.201'S 54°1.698'W) Entrance 26.50+0.11a 2.6740.89ab 59.93+2.03a Outlet 27.40+0.11b 5.06+40.89bc 60.95+2.03a Control 27.30+0.11b 2.2140.89a 58.00+2.48a Table 2. Relative abundance of macroinvertebrates (no insects) collected with a Surber net in the rice agroecosystem, in Treinta y Tres, Uruguay. Class Order Family Genus/morpho Relative species abundance (%) Maxillopoda Cyclopoida Cyclopidae Cyclopidae 22.50 morphospecies 1 Branchiopoda Anomopoda Daphniidae Daphnia 18.14 morphospecies 1 Arachnida Trombidiformes Hydrachnidae Hydrachna 8.00 (subclass Acarina) Hydrozetidae“ morphospecies 1 undetermined Hydrozetidae morphospecies 1 Acarina morphospecies 1 Gastropoda Caenogastropoda Ampullariidae Pomacea and other 4.88 (snails) spp. Oligochaeta Tubificida undetermined Tubificida 1.80 morphospecies 1 Aquatic macroinvertebrates in Uruguayan rice agroecosystem 7 Class Order Family Genus/morpho Relative species abundance (%) Malacostraca Amphipoda Hyalellidae Hyalella 1.39 morphospecies 1* Arachnida Araneae Allocosinae, Lyniphiidae, juveniles 0.33 Anyphaenidae Clitellata Hirudinea undetermined unidentified leech 0.29 ‘ indicates uncertain identification at Genus level. A total of 2820 macroinvertebrates in 27 samples were collected from the three localities. Insects represented 41.5% of total macroinvertebrates collected and Maxillopoda (Crustacea), Branchiopoda (Crustacea) and Arachnida were the most abundant macroinvertebrates classes collected after insects (Table 2, Fig. 1). From the total macroinvertebrates collected, 1170 were insects, belonging to 10 orders, 39 families and 50 species/morphospecies (Table 3). The more abundant insect orders were Diptera (59.9%), Hemiptera (16.3%) and Ephemeroptera (14.0%). Chironomidae and Syrphidae were the main dipteran families with 43.9% and 13.8% of the total insect abundance, respectively. Corixidae was the most abundant hemipteran (15.0%), while Ephemeroptera was represented by Baetidae and Caenidae with 7.8% and 6.2% of the total insect abundance, respectively. The more diverse insect groups were Coleoptera with 12 families and 16 morphospecies and Diptera with seven families and 10 morphospecies. 70 Species richness NO oO 012 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 Number of samples —Sobs -Singletons ~Doubletons -Uniques Duplicates Chao 1 «Chao 2 Jack 1 «Jack 2 Bootstrap *ACE Figure 3. EES Cumulative Species Curve for aquatic insects, indicating observed richness (Sobs), singletons, doubletons, duplicates, uniques and richness estimators (500 aleatorisations). 8 Bao L etal Table 3. Insects collected with a Surber net in the rice agroecosystem at the entrance, outlet and control zones of water in Treinta y Tres. Uruguay. Order Family Genus/Species/morphospecies Ephemeroptera Baetidae Americabaetis morphospecies 1% Caenidae Caenis morphospecies 1 Odonata Libellulidae Libellulidae morphospecies 1 Coenagrionidae Acanthagrion morphospecies 1 * Trichoptera Hydroptylidae Oxyethira morphospecies 1 Hemiptera Belostomatidae Belostoma morphospecies 1 Corixidae Sigara chrostowskii Corixidae morphospecies 1 Notonectidae Notonectidae morphospecies 1 Nepidae Ranatra morphospecies 1 Aphididae Rhopalosiphum padi* Cicadellidae Cicadellidae morphospecies 1 (nymph)* Pentatomidae Pentatomidae morphospecies 1 (nymph)* Lygaeidae Nysius simulans* Cercopidae Cercopidae morphospecies 1 (nymph)* Hymenoptera Formicidae Solenopsis morphospecies 1* Coleoptera Hydrophylidae Berosus morphospecies 1 Dytiscidae Curculionidae Scarabaeidae Tropisternus morphospecies 1 Enochrus morphospecies 1 Laccophillus morphospecies 1 Rhantus morphospecies 1 Desmopachria morphospecies 1* Oryzophagus oryzae Hypselus ater Ataenius morphospecies 1 Aphodius morphospecies 1 Haliplidae Haliplus morphospecies 1 Carabidae Meotachys morphospecies 1** Elmidae Elmidae morphospecies 1 Ptinidae Scydmaenus morphospecies 1 * Chrysomelidae Mycetophagidae Staphylinidae Epitrix morphospecies 1* Mycetophagidae morphospecies 1* Staphylinidae morphospecies 1* Aquatic macroinvertebrates in Uruguayan rice agroecosystem 9 Order Family Genus/Species/morphospecies Noteridae Suphisellus morphospecies 1 Diptera Chironomidae Chironomidae morphospecies 1 Chironomidae morphospecies 2 Chironomidae morphospecies 3 Chironomidae morphospecies 4 Chironomidae morphospecies 5 Syrphidae Eristalinae morphospecies 1 Ceratopogonidae Bezzia morphospecies 1 Culicidae Anopheles morphospecies 1 Simuliidae Simuliidae morphospecies 1 Empididae Empididae morphospecies 1 Lepidoptera Piralidae Parapoynx morphospecies 1 Hesperidae Hylephila morphospecies 1 Orthoptera Tettigoniidae Conocephalus morphospecies 1* Acridoidea subfamily Acridinae morphospecies 1* Thysanoptera Thripidae Thripidae morphospecies 1* Psocoptera Psocoptera morphospecies 1* * terrestrial groups not included in analysis. ‘ indicates uncertain identification at this level. The Cumulative Species Curve for insects shows that there are still species pending collection. Richness estimators show that at least 58.8% of expected species were sampled (Jackknife 2: 58.8%, Chao 1: 59.8%, Chao 2: 60.4%, ACE: 62.1%, Jackknife 1: 66.4%, Bootstrap: 74.4%, corresponding to an estimated number of species of 49, 48, 48, 47, 44 and 39, respectively) (Fig. 3). The physicochemical properties of water were correlated with functional groups according to the Mantel Test through Euclidean distances (r = 0.48 p = 0.019), but due to the number of samples, it was not possible to correlate each property to a particular functional group. The Kruskal-Wallis Test for total arthropods number showed the higher values from the water outlet and control zone, with a lower count for the water entrance (H = 18.47, p = 0.0178). Following the same analysis and comparing at family level, there were also differences found between water sources for Syrphidae larvae, with higher abundances in the water outlet (H = 20.34, p = 0.006) (Table 4). Taxa distributions by site and water source are shown as a heat map (Fig. 4). Despite all groups being registered in all the environments, some taxa were more abundant in certain water sources or sites. For example, Ephemeroptera were more abundant at the water entrances, while Syrphidae were more abundant in water outlet samples and Corixidae were more abundant in the control zones. 10 Bao L etal Table 4. Abundance of Syrphidae larvae and total number of arthropods according to water source (Entrance, Outlet and Control) in the rice agroecosystem in Treinta y Tres. Uruguay. Water source Syrphidae larvae Total number of arthropods Entrance 0.67 + 1.32a 40.78 + 6.87a Outlet 15.56 + 11.54b 89.11 + 16.61b Control 1.22 + 1.30a 145.78 + 18.56c The same letter within each column means no significant differences LSD Fisher Test (p * 0.05). ETE CHE CHC w 2 3 YMSE Ww 2 CHO © = ETO ETC JMSC JMSO Co Da Ch Ep Oa Am Hy Sy Mi Ga Cx Taxa Figure 4. EES] Heat map of taxon relative abundance and hierarchical clustering of samples. References: Sites: CH: La Charqueada, JMS: Julio Maria Sanz, ET: El Tigre; Water source: E: entrance, O: outlet, C: control; Taxonomic group: Co: copepods, Da: Daphnia, Ch: Chironomidae, Ep: Ephemeroptera, Oa: Other aquatics, Am: Amphipoda, Hy: Hydrophilidae, Sy: Syrphidae, Mi: Mites, Ga: Gastropoda, Cx: Corixidae. Values are scaled by taxon relative abundances across each site and water source. The ANOSIM Test showed differences according to locality and water source (R = 0.388, p = 0.0013 and R = 0.615, p = 0.0001, respectively). According to the SIMPER Test through Bray-Curtis distances, 80% of the differences between localities was explained by seven species, while differences between water sources were explained by six species (Table 5). Aquatic macroinvertebrates in Uruguayan rice agroecosystem 11 Table 5. SIMPER Test with species accomplishing 80% dissimilarity between localities and between water circulation in the rice agroecosystem in Treinta y Tres, Uruguay. SIMPER Test by LOCALITY Taxon Dissimilarity % % Mean by site Contribution Cumulative e JM EI Charqueada Sanz Tigre Chironomidae morphospecies 2. 22.83 27.94 27.94 13.60 15.20 7.67 Sigara chrostowskii 13.64 16.70 44.64 14.80 0.00 3.67 Eristalinae morphospecies 1 8.08 9.88 54.52 7.89 2.11 1.89 Caenis morphospecies 1 6.47 7.92 62.44 5.33 0.11 2.22 Americabaetis morphospecies 1 6.12 7.50 69.93 3.11 1.11 4.33 A Chironomidae morphospecies 3 6.07 7.43 77.36 1.89 3.11 1.00 Oryzophagus oryzae 2.02 2.47 79.83 0.56 0.22 0.67 SIMPER Test by WATER SOURCE Taxon Dissimilarity % % Mean by water source Contribution |= Cumulative ennence Gute’. | Cantrat Chironomidae morphospecies 2. 26.06 31.02 31.02 5.11 25.10 6.22 Sigara chrostowskii 13.30 15.83 46.85 2.56 0.11 15.80 Eristalinae morphospecies 1 8.34 9.93 56.79 1.78 9.11 1.00 Caenis morphospecies 1 6.76 8.05 64.84 7.56 0.11 0.00 Americabaetis morphospecies 6.40 7.62 72.46 4.11 4.11 0.33 4a Chironomidae morphospecies 3 6.26 7.45 79.91 3.33 2.33 0.33 ‘ indicates uncertain identification at Genus level. The Principal Component Analysis (PCA) shows samples grouping according to the water source (Fig. 5). Chironomidae morphospecies 2 and Syrphidae (Diptera) were associated with water outlet samples, while Caenis (Ephemeroptera) was associated with water entrance samples and corixids (Sigara chrostowskil) to the water control zone. Richness and Diversity Indices are shown in Table 6. Values for control zones did not show differences between sites, while, in the case of the water inlet, higher values of agreement were registered for El Tigre. Water outlet indices showed differences in some values, with the highest specific diversity in El Tigre (mean of 13 species). 12 Component 2 Figure 5. EES gara chrostowskii Bao L etal 204 Corixidae -Smorphospecies PURHHES -20° Component 1 PCA test with Biplot from variance-covariance matrix of Bray-Curtis distances for each locality and water source. (Control: C in green. Entrance: E in sky-blue and Outlet: O in blue). Component 1: 57.74% and Component 2: 21.50% of variance, respectively. Table 6. Richness and Diversity Indices for insects + E.E, by locality (site) in water samples. Diversity Richness (S) Shannon H Simpson 1-D Equitability J Richness (S) Shannon H Simpson 1-D Equitability J Richness (S) Shannon H Simpson 1-D Equitability J EI Tigre Entrance 6.33 + 1.53aA 1.55 + 0.12aA 0.76 + 0.03aA 0.81 + 0.03abA Outlet 13.33 + 0.58bB 2.02 + 0.07bB 0.62 + 0.04aA 0.70 + 0.07aA Control 11.00 + 2.65aB 1.40 + 0.15aA 0.48 + 0.08aA 0.56 + 0.09aA J.M. Sanz 9.00 + 1.73aB 1.54 +0.12aA 0.74 + 0.03aB 0.85 + 0.03bB 4.67 +1.15aA 1.144 0.07aA 0.81 + 0.04bB 0.78 + 0.07aB 7.00 + 1.73aAB 1.19+0.15aA 0.62 + 0.08aA 0.59 + 0.09aA La Charqueada 7.67 + 1.53aA 1.64+0.12aA 0.76 + 0.03aA 0.700 + 0.03aA 5.33 + 1.53aA 1.14+0.07aA 0.61 + 0.04aA 0.77 + 0.07aA 7.33 + 2.31aA 1.05 + 0.15aA 0.58 + 0.08aA 0.62 + 0.09aA The same capital letter within each column and same lower-case letter within each row means differences in LSD Fisher Test (p ” 0.05). no significant Aquatic macroinvertebrates in Uruguayan rice agroecosystem 13 Comparing values within each field according to water source, the major differences were registered for El Tigre. In J.M. Sanz, differences were registered only for richness with the higher value at the water entrance. On the other hand, for Charqueada, no differences were registered between water sources. Discussion This study represents the first description of the aquatic macroinvertebrate community associated with rice crops in Uruguay and establishes a baseline for present and future studies on productivity and ecological sustainability of Uruguayan rice farms. In relation to water physicochemical properties, the greatest differences were registered for dissolved oxygen with lower values in control zones from El Tigre and J.M. Sanz. Control zones are water bodies located in crop surroundings and with less area and without water flux compared with rice areas. In the case of La Charqueada, the control zone belongs to a more extensive water body. Less oxygen content in El Tigre and J. M. Sanz could be due to water body characteristics of these localities compared to La Charqueada control zone, such as the size and sun radiation or vegetable coverage, for example (Ito 1987). The Mantel Test demonstrated an effect of physicochemical properties on the guild composition in different water sources. Richness estimators show that there are still many species remaining to be collected and this work recorded at least 58% of the species, depending on the estimator considered. However, even under these conditions, from the 50 morphospecies collected, 36 were aquatic insects, showing a greater diversity than reported previously in tropical-subtropical rice fields of Brazil (84 morphospecies, Molozzi et al. 2007), but less diversity than that found in Costa Rica (44 insect morphospecies, Rizo-Patron et al. 2013) and Indonesia (45 insect morphospecies, Wakhid et al. 2020). However, our data are the results of sampling only at the grain-filling stage, thus biodiversity could be underestimated. It is important to emphasise that Uruguayan rice crops are mainly managed in a rotation system with pastures, which could bring more stability to the agroecosystem and sustain greater diversity (Norton et al. 2010). In addition, this production system has low inputs of insecticides which could also contribute to retaining biodiversity and natural pest control (Ali et al. 2019). Insecta was the most represented Class in the samples, with Diptera, Hemiptera and Ephemeroptera as the most represented Orders. The most abundant aquatic insects found in the present work were Chironomidae (Diptera), the most ubiquitous family of aquatic flies with a worldwide distribution (Ferrington 2008, Dijkstra et al. 2014). Chironomids include species which are detritivores and can function as alternative food sources for generalist predators like Arachnids and Odonata (Ferrington 2008; Wakhid et al. 2020). Moreover, further work on chironomid species identification would allow them to be assessed for water quality monitoring (Takahashi et al. 2008). The ANOSIM Test proved that there were differences in morphospecies composition according to locality and water source (R = 0.388, p = 0.0013 and R = 0.615, p = 0.0001, 14 Bao L etal respectively). Some groups showed association with a particular water source, so they could be analysed for their potential as water quality indicators for each source studied (Entrance, Outlet and Control zone). Caenis (Ephemeroptera) individuals are associated with the water entrance. It is well documented that Caenis individuals are sensitive to pollution (Roldan 2003, Ocon and Rodriguez-Capitulo 2012), so the greater amounts observed in water entrances could be indicating better water quality compared to outlet and control zones. On the contrary, a higher abundance of Syrphidae larvae in water outlets could be indicating water deterioration in the path through the crop or soil and plant debris addition from the rice crop due to these organisms being associated with an excess of organic matter (Bonada et al. 2006, Arimoro et al. 2007). However, in control zones, corixids were more abundant and the species varied between sampling localities. Corixids are classified as predators and several species are pioneers that quickly colonise new habitats, occasionally newly-planted rice fields (Wakhid et al. 2020). Coleoptera was the most diverse order (18 morphospecies) and Dytiscidae was the most diverse family (four morphospecies). These results are in accordance with data reported by Gomez Lutz et al. (2015) for Argentinean rice fields that reported Dytiscidae and Hydrophilidae as the most diverse coleopteran families. Taking into account these results, further sampling exercises could probably register an even higher number of coleopteran species for Uruguayan rice fields. The Shannon-Weaver Indices, calculated for locality and water source with aquatic insects species, did not show differences between water sources. However, richness values registered for aquatic insect species are higher than those registered in organic rice crops in Australia (Wilson et al. 2008), but lower than those of Costa Rica and Italy (Rizo-Patron et al. 2013, Lupi et al. 2013Wilson et al. 2008). Considering richness estimators calculated in this work and the fact that sampling was only during the grain-filling stage, we expect to find more species with future research. In relation to similar research from Uruguay, Shannon Indices are similar to those reported by Morelli and Verdi (2014) in watercourses from riverine native vegetation, despite the fact that species composition was different. The higher richness and diversity registered in El Tigre could be due to different factors. From the point of view of the landscape where this crop is located, probably there is a better connectivity with native vegetation patches than in the other two sites (J.M. Sanz and La Charqueada), but this has to be confirmed. It is possible that this environment, with the largest area of native riverine vegetation of the three sampled sites, could be a better reservoir for species, supporting a higher diversity and richness (Molozzi et al. 2007, Gonzalez et al. 2016). Moreover, it should be emphasised that sampling was undertaken at the grain-filling stage, so there is no information about possible diversity fluctuation along the crop cycle, as reported previously (Che Salmah et al. 2017). Under these sampling conditions, it is important to discuss the fact that diversity could have been underestimated and there could be still many macroinvertebrate species associated with this crop that are pending to be sampled. A more exhaustive sampling at different crop stages could detect more Aquatic macroinvertebrates in Uruguayan rice agroecosystem 15 differences in macroinvertebrate community composition and register an even higher number of species for Uruguayan rice ecosystems. The progress in aquatic macroinvertebrate community characterisation and agrochemical residues evaluations could enable the definition of indicators as relatively quick sampling tools for water quality screening, thus enabling the measurement of the impact of crop management activities. Conclusions The physicochemical properties of water had an effect on the composition of the guilds. Different morphospecies compositions were found according to locality and water source. Caenis individuals were associated with the water entrance, while Syrphidae and Chironomidae (morphospecies 2) larvae were associated with the water outlet. Some of these species could be useful indicators of water quality. Arthropods were more abundant in the water of control zones. More exhaustive sampling could register more diversity and potential bio-indicators. Richness and diversity indices registered in this work - even with an inventory of species pending registration - are higher than those reported for rice agroecosystems from Costa Rica, Italy and Australia. Acknowledgements Fellowship ANI] POSNAC_2012_4459 to L. Bao. We thank Maria Caraballo for her collaboration in field sampling, to Héctor Da Fonseca and Raul Servetto for the possibility of sampling in their fields and Marcelo Segovia for his assistance in crop selection and providing field management information. Author contributions Leticia Bao conceptualised the study; Leticia Bao, Monica Urrutia and Lucia Seijas contributed to the methodology, Leticia Bao, Monica Cadenazzi, Enrique Castiglioni and Sebastian Martinez carried out formal analysis and investigation; Leticia Bao, Enrique Castiglioni and Sebastian Martinez worked in writing, reviewing and editing. Leticia Bao supervised the study. Conflicts of interest There are no financial support or relationships that may pose any conflict of interest. 16 Bao L etal References ° Ali MP, Bari MN, Haque SS, Kabir MM, Afrinf S, Norwin F, Islam MS, Landis DA (2019) Establishing next-generation pest control in rice fields: eco-agriculture. Scientific Reports 9: 10180. https://doi.org/10.1038/s41598-019-46688-6 ° Arimoro FO, lkomi RB, lwegbue CM (2007) Water quality changes in relation to Diptera community patterns and diversity measured at an organic effluent impacted stream in the Niger Delta, Nigeria. Ecological Indicators 7 (3): 541-552. https://doi.org/10.1016/ j.ecolind.2006.06.002 ° Balian EV, Segers H, Lvque C, Martens K (2008) The freshwater animal diversity assessment: an overview of the results. Hydrobiologia 595: 627-637. htips://doi.org/ 10.1007/978-1-4020-8259-7 61 ° Bao L, Martinez S (2018) Control quimico de insectos en el cultivo de arroz en Uruguay. Arroz 96: 40-45. [In Spanish]. URL: http://www.ainfo.inia.uy/digital/bitstream/ item/12366/1/Arroz-2018-96-Bao.pdf ° Bentancourt CM, Scatoni IB, Morelli E (2009) Insectos del Uruguay. [Insects of Uruguay]. Universidad de la Republica, Facultad de Agronomia — Facultad de Ciencias, ANII-CSIC, Montevideo, 658 pp. ° Blitzer EJ, Dormann CF, Holzschuh A, Klein AM, Rand TA, Tscharntke T (2012) Spillover of functionally important organisms between managed and natural habitats. Agriculture, Ecosystems & Environment 146 (1): 34-43. https://doi.org/10.1016/j.agee.2011.09.005 ° Bonada N, Prat N, Resh VH, B S (2006) Developments in aquatic insect biomonitoring: A comparative analysis of recent approaches. Annual Review of Entomology 51 (1): 495-523. https://doi.org/10.1146/annurev.ento.51.110104.151124 ° Che Salmah MR, Siregar AZ, Abu Hassan A, Nasution Z (2017) Dynamics of aquatic organisms in a rice field ecosystem: effects of seasons and cultivation phases on abundance and predator-prey interactions. Tropical Ecology 58 (1): 177-191. ° Colwell RK (2013) EstimateS: Statistical estimation of species richness and shared species from samples. Version 9. ° Dijkstra KD, Monaghan MT, Pauls SU (2014) Freshwater biodiversity and aquatic insect diversification. Annual Review of Entomology 59: 143-163. https://doi.org/10.1146/ annurev-ento-011613-161958 ° Di Rienzo JA, Casanoves F, Balzarini MG, Gonzalez L, Tablada M, Robledo CW (2017) InfoStat. Version 2017. Grupo InfoStat, FCA, Universidad Nacional de Cordoba, Argentina. URL: http:/Awww.infostat.com.ar ° Ferrington LC (2008) Global diversity of non-biting midges (Chironomidae; Insecta- Diptera) in freshwater. Hydrobiologia 595: 447-455. httos://doi.org/10.1007/ $10750-007-9130-1 ° Gomez Lutz MC, Kehr Al, Fernndez LA (2015) Abundance, diversity and community characterization of aquatic Coleoptera in a rice field of Northeastern Argentina. Revista de Biologa Tropical 63 (3): 629-638. ° Gonzalez E, Salvo M, Defag MT, Valladares G (2016) A moveable feast: moving at the forest-crop interface are affected by crop phenology and the amount of forest in the landscape. PLOS One 11 (7): 1-19. https://doi.org/10.1371/journal.pone.0158836 Aquatic macroinvertebrates in Uruguayan rice agroecosystem 17 Gonzalez ML, Simone MJ, Morais RM, Da Silva G (2014) Diversidad de insectos depredadores en rea orizcola orgnica y de conservacin en Viamo, RS. Brasil, Revista Colombiana de Entomologa 40 (1): 120-128. Gurr G, Wratten S, Snyder W (2012) Biodiversity and insect pests: Key issues for sustainable management. Wiley-Blackwell Ed, 34 pp. https://doi.org/10.1002/ 9781118231838 Hammer T, Harper DA, Ryan PD (2001) PAST: Paleontological statistics software package for education and data analysis. Palaeontologia Electronica 4 (1). URL: http://palaeo-electronica.org/2001 1/past/issue1 01.htm Ito A (1987) Changes of water temperature, pH, dissolved oxygen, inorganic nitrogen, and phosphorus concentrations in flowing irrigation water on paddy surface. Soil Science in Plant Nutrition 33 (3): 449-459. https://doi.org/10.1080/00380768. 1987.10557591 Luo Y, Fu H, Traore S (2014) Biodiversity conservation in rice paddies in China: Toward ecological sustainability. Sustainability 6: 6017-6024. https://doi.org/10.3390/su6096107 Lupi D, Rocco A, Rossaro B (2013) Benthic macroinvertebrates in Italian rice fields. Journal of Limnology 72 (1): 184-200. https://doi.org/10.4081/jlimnol.2013.e15 Maltchik L, Rolon AS, Stenert C, Machado IF, Rocha O (2011) Can rice field contribute to biodiversity in Southern Brazilian wetlands? Journal of Tropical Biology 59 (4): 1895-1914. Melo S, Stenert C, Dalzochio MS, L M (2015) Development of a multimetric index based on aquatic macroinvertebrate communities to assess water quality of rice fields in southern Brazil. Hydrobiologia 742: 1-14. https://doi.org/10.1007/s10750-014-1957-7 Merritt RW, Cummins KW (2007) Trophic relationships of macroinvertebrates. In: Hauer FR, Lamberti GA (Eds) Methods in stream ecology. Second Edition. Elsevier, 585-601 pp. https://doi.org/10.1016/B978-012332908-0.50036-X Molozzi J, Hepp LU, Dias AdS (2007) Influence of rice crop on the benthic community in Itaja Valley (Santa Catarina, Brazil). Acta Limnologica Brasileira 19 (4): 383-392. Morelli E, Verdi A (2014) Diversidad de macroinvertebrados acuticos en cursos de agua dulce con vegetacin riberea nativa de Uruguay. Revista Mexicana de Biodiversidad 85: 1160-1170. https://doi.org/10.7550/rmb.45419 Norton GW, Heong KL, Johnson D, Savary S (2010) Rice pest management: issues and opportunities. In: Pandey S, Byerlee D, Dawe D, Dobermann A, Mohanty S, Rozelle S, Hardy B (Eds) Rice in the global economy: strategic research and policy issues for food security. International Rice Research Institute, Los Banos; Philippines, 297-332 pp. Ocon C, Rodriguez-Capitulo A (2012) Assessment of water quality in temperate-plain streams (Argentina, South America) using a multiple approach. Ecologa Austral 22: 81-91. Pires MM, Kotzian CB, Spies MR, Baptista VdA (2015) Comparative assessment of aquatic macroinvertebrate diversity in irrigated rice fields and wetlands through different spatial scales: an additive partitioning approach. Marine and Freshwater Research 67: 368-379. [In English]. https://doi.org/10.1071/MF14109 Ramirez A, Gutiérrez-Fonseca PE (2014) Functional feeding groups of aquatic insect families in Latin America: a critical analysis and review of existing literature. Revista de Biologa Tropical 62 (2): 153-165. Ramsar Convention Wetlands (2018) Global Wetland Outlook: State of the Worlds Wetlands and their Services to People. Ramsar Convention Secretariat, 84 pp. 18 Bao L et al R Core Team (2020) R: A language and environment for statistical computing. R Foundation for Statistical Computing. R Core Team. URL: https://www.R-project.org/ Rizo-Patron F, Kumar A, McCoy Colton MB, Springer M, Trama FA (2013) Macroinvertebrate communities as bioindicators of water quality in conventional and organic irrigated rice fields in Guanacaste, Costa Rica. Ecological Indicators 29: 68-78. https://doi.org/10.1016/j.ecolind.2012.12.013 Roldan GA (2003) Bioindicacion de la calidad del agua en Colombia. Primera Edicion, Editorial Universidad de Antioquia, 170 pp. SAS Institute (2009) SAS/STAT 9.2 User's guide. Version 9.2. SAS Publishing. Settle W (1996) Managing tropical rice pests through conservation of generalist natural enemies and alternative prey. Ecology 77 (7): 1975-1988. https://doi.org/ 10.2307/2265694 Springer M (2006) Clave taxonomica para larvas de las familias del orden Trichoptera (Insecta) de Costa Rica. Revista de Biologa Tropical 54 (1): 273-286. Springer M (2008) Aquatic insect diversity of Costa Rica: State of knowledge. Revista de Biologa Tropical 56 (4): 273-295. Stenert C, Bacca R, Maltchik L, Rocha O (2009) Can hydrologic management practices of rice fields contribute to macroinvertebrate conservation in southern Brazil wetlands? Hydrobiologia 635: 339-350. https://doi.org/10.1007/s10750-009-9926-2 Takahashi MA, Higuti J, Bagatini YM, Zviejkowski IP, Velho LF (2008) Composition and biomass of larval chironomid (Insecta. Diptera) as potential indicator of trophic conditions in southern Brazil reservoirs. Acta Limnologica Brasileira 20 (1): 5-13. Tseng MC, Roel A, Deambrosi E, Terra JA, Zorrilla G, Riccetto S, Pittelkow CM (2020) Towards actionable research frameworks for sustainable intensification in high-yielding rice systems. Scientific Reports 10 (9975). https://doi.org/10.1038/s41598-020-63251-w Vilaboa N (2012) Evaluacion de calidad de agua en un sistema lotico bajo impacto lechero mediante el uso de macroinvertebrados acuaticos. Tesina de Grado, Ciencias Bioldgicas Facultad de Ciencias, Universidad de la Republica, Uruguay, Montevideo, 54 pp. Wakhid W, Rauf A, Krisanti M, Sumertajaya M, Maryana N (2020) Species richness and diversity of aquatic insects inhabiting rice fields in Bogor, West Java, Indonesia. Biodiversitas Journal of Biological Diversity 21 (1): 34-42. https://doi.org/10.13057/ biodiv/d210106 Wilson AL, Watts RJ, Stevens MM (2008) Effects of different management regimes on aquatic macroinvertebrate diversity in Australian rice fields. Ecological Research 23: 565-572. https://doi.org/10.1007/s11284-007-0410-z