Geographical origin determines responses to salinity of Mediterranean caddisflies.

Many freshwater ecosystems worldwide, and particularly Mediterranean ones, show increasing levels of salinity. These changes in water conditions could affect abundance and distribution of inhabiting species as well as the provision of ecosystem services. In this study we conduct laboratory experiments using the macroinvertebrate Smicridea annulicornis as a model organism. Our factorial experiments were designed to evaluate the effects of geographical origin of organisms and salinity levels on survival and behavioral responses of caddisflies. The experimental organisms were captured from rivers belonging to three hydrological basins along a 450 Km latitudinal gradient in the Mediterranean region of Chile. Animals were exposed to three conductivity levels, from 180 to 1400 μS/cm, close to the historical averages of the source rivers. We measured the behavioral responses to experimental stimuli and the survival time. Our results showed that geographical origin shaped the behavioral and survival responses to salinity. In particular, survival and activity decreased more strongly with increasing salinity in organisms coming from more dilute waters. This suggests local adaptation to be determinant for salinity responses in this benthic invertebrate species. In the current scenario of fast temporal and spatial changes in water levels and salt concentration, the conservation of geographic intra-specific variation of aquatic species is crucial for lowering the risk of salinity-driven biodiversity loss.

Introduction

Increases in water salinity affect freshwater ecosystems in many regions of the globe [1, 2]. These changes in ionic concentration and composition impair water quality, which disrupts the abundance and distribution of aquatic species [1, 3–8]. This disruption threatens ecosystem functioning and the provision of services of crucial importance for our society such as drinking and irrigation water, climate regulation, recreation and food [9]. Mediterranean regions are specially vulnerable to this stress due to their marked seasonality in environmental temperature and evaporation.

To date, our understanding about the physiological mechanisms by which freshwater organisms deal with salinity stress is limited [2, 10]. Moreover, our knowledge of the consequences of salinity changes for organisms performance is relatively poor, considering the attention given to other stressors in freshwaters, such as metals and pesticides [11]. Salinity changes may exert direct lethal effects if those changes are strong enough as to exceed the tolerance limits of organisms [10]. Nevertheless, if salinity changes are more subtle, sublethal effects may be expressed such as alterations in development, physiology or behavior [12-17]. In particular if behavior is affected, this would imply a dysfunction in the short-term behavioral capabilities of reacting to environmental signals and stimuli. In benthic freshwater invertebrates, behavioral responses allow them to perform habitat selection, through deciding either to stay in the same place, to perform small—scale movements or to release from the substrate to derive and move away quickly [7, 18]. These behavioral responses to environmental risk are vital for organism performance and population persistence.

Freshwaters from Mediterranean regions are of special concern for biodiversity maintenance and ecosystem functioning. This is due to that Mediterranean ecosystems, on one hand, exhibit high levels of biodiversity and endemicity and, on the other hand, are subjected to high levels of biological invasions and are expected to be highly affected by climate change [19]. Rivers from Mediterranean Chile are characterized by strong environmental gradients [20-22], where both altitude and latitude are strong predictors of abiotic conditions.

Besides this, hydrological basins in this region have been threatened by economic activities mainly derived from mining and agriculture industry during the last century. These actions have disrupted physical, chemical and biological properties of these ecosystems [4, 23]. Thus, both natural –physical and geological– features as well as anthropogenic factors have contributed to create contrasting salinity levels in rivers belonging to the Mediterranean region of central Chile that extends from 31°S to 35°S along the Andean hydrological basins.

Physiological adaptation to different levels of water salinity has been reported for aquatic freshwater organisms, mainly linked with mechanisms of ion regulation [11, 24]. Thus, organisms from different local populations of a species may exhibit different degrees of salinity tolerance, which could also be expressed into diverse abilities to respond behaviorally to external stimuli. Nevertheless, due to our limited understanding of the ecological consequences of osmotic stress, it is unclear how the phenotypic responses to different salinity levels are shaped by local adaptation.

In this study we evaluate experimentally the effects of geographical origin and local salinity levels on survival and behavioral responses of a model benthic freshwater species inhabiting the Mediterranean region of Chile, the caddisfly Smicridea annulicornis. Our general prediction is that survival and behavioral responses to different salinity levels will be function of the geographical origin of the organisms. However, we have specific expectations according to two competing hypotheses. Habitats with lower water salinity represent a challenging environment for benthic invertebrate species, due to difficulties for osmotic regulation in a medium with a lower ionic concentration than internal fluids [25]. Thus, it has been suggested that the ability of some invertebrates to tolerate dilute water may allow them to inhabit a wide range of water salinity. However, physiological mechanisms to overcome low water salinity likely evolve at the expense of other traits that are beneficial at higher water salinity conditions [11, 26, 27]. According to this hypothesis, we expect that organisms adapted to low salinity will maintain similar levels of survival through different experimental salinity levels, but will reduce their behavioral performance at higher salinity as an expression of the costs of low-salinity adaptation. Also, we expect that organisms adapted to higher salinity levels will reduce their survival under low salinity. On the other hand, several studies indicate that freshwater organisms are negatively affected by increased salinity [28, 29]. According to this hypothesis, we would expect that low-salinity adapted organisms will be more affected by higher salinity levels both in behavioral responses and survival, as compared to organisms adapted to higher salinity levels. Under this hypothesis, low-salinity adapted organisms should be more sensitive to experimental changes in salinity.

Our analyses shed light into how benthic freshwater fauna may respond to current salinity changes by evaluating an array of lethal (mortality) and nonlethal (behavioral responses) attributes.

Materials and methods

Study organisms

We conducted this study using as a model species the Neotropical trichopteran Smicridea annulicornis (Blanchard 1851), which has a broad geographic distribution in latitude and altitude [5]. The larval stages of this species inhabit the rhithral areas of rapids and pools where they are abundant and play a functional role as filter—feeders/collectors, consuming coarse organic particles and algae. During spring, large numbers of adults emerge and fly away [13, 30]. However genetic analyses indicated a low vagility of this species [31]. We collected the animals from three hydrological basins along a 450 km latitudinal gradient in the Mediterranean zone of central Chile: Choapa (S31° 52’; W70° 57’), Maipo (S33° 22’; W70° 27’) and Maule (S35° 33’; W7‘1° 13’) (Fig 1). Biological samples were collected from free access riparian areas in Chilean rivers (DL N° 1.939-1977). In these areas, no permits are required for sampling invertebrate animals. In each of these localities, larvae of S. annulicornis were collected by hand from the riverbed in the spring-summer of 2016, and placed together in plastic jars of 1L to be transported to the laboratory. All animals were kept in water from the same place where the animals were collected and they were acclimated to the experimental conditions during one week before performing the measurements. We used electric conductivity as proxy of salinity. To determine the levels of conductivity to be used in the experiments we considered the range of variation of water conductivity at each of the hydrological basins reported for the last 50 years (Dirección General de Aguas, Chile; www.snia.dga.cl). We choose three experimental conductivity levels: low (180 μS/cm), medium (500 μS/cm) and high (1400 μS/cm). These levels were close to the ones measured by us in situ and approximate the conductivity values found historically at Choapa, Maipo and Maule basins respectively. Then, organisms from each population were randomly transferred individually into 40 mL glass jars to each experimental conductivity level for conducting the factorial experiment. Different salinity levels of the experimental media were prepared with aquarium sea salt (Sera®) dissolved in distilled water up to reach the appropriate conductivity values, measured with a desktop electrical conductivity meter (Bante Instruments®). The media were renewed every two days during the experimental period to avoid oxygen depletion and changes in pH. These parameters were recorded in parallel during the experiment.

Fig 1

Main features of the study sites.

A: Geographical location of the hydrological basins Choapa, Maipo and Maule in central Chile. B: Values of electric conductivities recorded for the three hydrological basins by the DGA (Dirección General de Agua, Chile; http://snia.dga.cl/BNAConsultas/reportes) at the monitoring water quality stations during the last 50 years. Horizontal dashed lines indicate the three experimental conductivity levels used in our study.

Measurement of response variables

During the experiment, all animals were kept in dechlorinated water, maintained at 20°C with a photoperiod of 16L:8D. Animals were acclimated during 48 hours prior experimentation under the same ambient conditions used during measurements. To avoid further variability due to food conditions and to strengthen the responses to experimental factors, animals were kept unfed during the experiment. To conduct our experiment, we used 85 larvae from Choapa basin, 81 from Maipo basin and 48 larvae from Maule basin, distributed uniformly over the three conductivity levels (low, medium, and high. The body size (length from the top of the head to the end of the tail) of all experimental organisms was measured from digital images obtained with a digital camera attached to a stereo—microscope. Individual survival in each treatment was measured as the time from the beginning of the experiment to the death of the animal, after which the organisms were kept individually in 70% ethanol to measure body size. The behavioral traits were defined as the response to direct stimuli, reflecting the capability of responding to an external disturbance. To trigger the behavioral responses we used a pencil brush, with which we gently touched each animal when it was in a resting state. We defined four response traits from this procedure: i) “swimming” shown as an active movement through the water column; ii) “sheltering”, characterized by assuming a spherical protection posture; iii) “walking” characterized by a regular displacement on the bottom of the glass using the front legs; and iv) “pushing-up” is characterized by a strong dorsal contraction of the body, forming an arched posture, at the same time that the animal sinks in the water column. We recorded the occurrence of each behavioral trait within a time window of 30 seconds of observations per day during the entire experiment (≈ two weeks). Because we allowed the animals to perform any of these behaviors in each time window of observation, we considered the expression of each behavior responses as independent.

Statistics

A factorial design was used to analyze the recorded responses. Our design included three levels of river basin and three levels of water conductivity. Survival analysis was conducted using a parametric model based on the Weibell distribution, with body size used as a covariate. For each trait, repeatability was evaluated by intra-class correlation coefficient [32] by calculating the between– and within—individual variance component, from models which consider the biased distribution of counts and proportional data [33]. We conducted the repeatability analysis using the rptR library in R software. The exhibition of any given behavior, here referred as “activity” was analyzed using a model which included body size, the time at measurement as well as the kind of behavioral trait adopted in each activity event recorded at each experimental condition. Because the behaviors were recorded as the presence or absence of the responses in discrete time (count data), our variable responses swimming, sheltering, pushing-up and walking were biased. Thus, to analyze the behavioral response variables, we fit Generalized Linear Models (GLMs) with error family according to this characteristic (quasibinomial). Body size and measurement times were used as covariates into the models to avoid under—sub estimating the behavioral responses to the combination of stress conditions. All the analyses were performed using R software [34].

Results

Survival time

Survival of the caddisflies Smicridea annulicornis differed significantly among organisms from different river basin. In addition there were nearly significant effects of conductivity (Table 1). The interaction between basin and conductivity did not exert significant effects on survival. The median survival time (time to 50% survival) for low, medium and high conductivity levels respectively, was 349 ± 45, 385 ± 24 and 197 ± 47 hours in caddisflies from Maule, 235 ± 26, 317 ± 22 and 214 ± 22 hours in caddisflies from Choapa, and 347 ± 29, 329 ± 33 and 349 ± 33 hours in caddisflies from Maipo (Fig 2). A visual analysis of the survival curves (Fig 2) suggests that conductivity level exerted: a) stronger effects on survival in organisms coming from Maule basin, with marked lower survival at high conductivity, b) lighter effects in caddisflies from Choapa, with higher survival at medium conductivity, and c) negligible effects in caddisflies from Maipo. These differences in organisms’ sensitivity to conductivity is more evident when considering the variance of median survival time among the three conductivity levels within each basin: 9957.3 (Maule), 2962 (Choapa), and 121.3 (Maipo). The patterns describe above indicate intrinsic phenotypic differences among the tested populations.

Table 1

ANOVA results from survival parametric analyses with Weibull distribution to three caddisfly populations (Choapa, Maipo and Maule) exposed to low (180 μS/cm), medium (500 μS/cm) and high (1400 μS/cm) conductivity levels, body size was used as a covariate.

Factor	−2* Log Likelihood	χ—value(df)	P-value
Null model	2760.77
Body size (mm)	2758.85	1.92(1,207)	0.165
Basin	2751.66	7.19(2,205)	0.027
Conductivity (μS/cm)	2746.67	4.98(2,203)	0.082
Basin x Conductivity	2741.12	5.55(4,199)	0.235

Fig 2

Survival curves.

Kaplan-Meir survival curves of S. annulicornis from populations belonging to Choapa, Maipo and Maule river basins subjected to conductivity treatments. Doted black line indicate the 50% of survival.

Behavioral responses

Activity (i.e. the rate of behavioral responses) was significantly affected by experimental time and body size (Table 2). In addition, organisms from different river basins showed significant differences in their activity levels in dependence on the specific behavior displayed (Table 2, Fig 3). In average, activity was lower in organisms from Choapa than Maule and it was the highest in organisms from Maipo. However, organisms from Maule basin exhibited lower levels of swimming behavior than organisms from both Maipo basin (Tukey test, P < 0.001) and Choapa basin (Tukey test, P < 0.05). Conversely, organisms from Maule basin exhibited higher levels of pushing-up than organisms from both Maipo (Tukey test, P < 0.001), and Choapa (Tukey test, P < 0.001). Pairwise comparison between Choapa and Maipo basins revealed lower levels of pushing-up in organisms from Maipo (Tukey test, z-ratio = 2.84, P < 0.05). Specific comparisons within Maule basin revealed that total activity as well as sheltering, walking and pushing-up behaviors diminished at high as compared to both low and medium conductivity levels, although these differences were not significant (Tukey test, P > 0.05). In Choapa basin, all behavioral traits exhibited higher activity rates at the medium conductivity level. These differences were statistically significant (P < 0.05) for total activity. In Maipo basin the behavioral responses expressed at similar levels irrespective of conductivity, with the exception of swimming behavior, which was marginally lower (P = 0.04) at low conductivity than at medium conductivity. The behavioral traits evaluated showed significant repeatability (Intraclass Correlation Coefficient—ICC), thought the magnitude was different among them. Swimming and pushing-up behaviors showed low repeatability as compared to walking and sheltering (S1 Table).

Table 2

Summary of results of ANOVA analyses with general linear model (GLMs) for total activity recorded from three caddisfly population (Choapa, Maipo and Maule) exposed to low (180 μS/cm), medium (500 μS/cm) and high (1400 μS/cm) conductivity levels.

Body mass and experimental time were used as covariates. Behavior categories were used as a factor to test if the different behaviors (swimming, sheltering, walking and pushing-up) differed in proportion among experimental treatments.

Factor	F-value(df)	P-value
Experimental time (days)	990(1,826)	<0.001
Body Size (mm)	22.2(1,825)	<0.001
Behavior Category (BC)	159(3,822)	<0.001
Basin	3.78(2,820)	0.021
Conductivity (μS/cm)	0.07(2,818)	0.923
BC x Basin	15.5(6,812)	<0.001
BC x Conductivity	0.87(2,806)	0.513
Basin x Conductivity	1.62(1,802)	0.166
BC x Basin x Conductivity	1.49(2,790)	0.123

Fig 3

Behavioral responses.

Mean total activity and its decomposition into specific behavioral responses exhibited by S. annulicornis. Values of activity levels are shown in parenthesis. Organisms came from three river basins: Choapa, Maipo and Maule. Experimental conductivity levels were: low (180 μS/cm), medium (500 μS/cm) and high (1400 μS/cm).

Discussion

Survival and behavioral responses of Smicridea annulicornis showed to be sensitive to water salinity in our laboratory treatments. However, the responses were stronger in organisms whose geographic origin was the river with lower historical levels of water salinity. This result support our general prediction that geographical origin shape the responses to water salinity. Although, this result contrasts with those of Kefford et al. [35], where organisms collected from different rivers having different conductivity ranges did not exhibit differences in their salinity tolerance. Differences in tolerance of invertebrates to salinity increments has been found associated to taxonomic affiliation [36] or habitat type (e.g. lentic vs lotic, [37]). However, there is little evidence on differences in salinity responses attributable to geographic origin [38]. This is surprising considering the current relevance of salinity pollution for the conservation of freshwater biodiversity [26]. Freshwater ecosystems suffering from anthropogenic salinization have shown effects not only on community structure but also on the phenotypic expression of reproductive, morphological and physiological traits [39].

Regarding our competing hypotheses about the way in which experimental organisms respond according to their geographical origin, our results supported our second hypothesis. This hypothesis proposes that low-salinity adapted organisms should be more sensitive to experimental changes in salinity, as compared to organisms adapted to higher salinity levels. This supports that high salinity levels, within the observed ranges in Mediterranean rivers, act as a stressing factor to benthic macroinvertebrate organisms [40]. These effects could be modulated by differential osmoregulatory capabilities due to local adaptation [25]. Many species facing osmotic stress have developed adaptations to deal with salinity variation in their environment [41]. Furthermore, there is evidence of rapid evolution of osmoregulatory adaptation in aquatic organisms, which indicates that local adaptation could occur within few generations in the wild [42]. Recent findings showed these adaptation being present in benthic macroinvertebrates facing sublethal salinity levels [43, 44]. The mechanisms involve adjustment (i.e. up and down regulation) of the abundance of specialized cells (ionocytes) as well as of the abundance of transporter proteins in their membranes. Thus, the conjecture of local adaptation to explain the observed effect of geographical origin in our study organisms appears plausible, regarding the current understanding of physiological mechanisms of osmotic regulation in a variety of aquatic organisms [45-47]. Also, our experimental results are in line with field information that indicates salinity increment in rivers having a detrimental effect on biodiversity [26]. Furthermore, increased salinity has also shown to induce changes in community structure [48], and affect the phenotypic variation of the freshwater species. In this vein, it is expected that selection for tolerant variants triggers biodiversity loss, both at intraspecific and interspecific levels [25, 26].

This study aids to expand our knowledge about how the organisms respond to environmental variation, specifically to salinity, and which is the role of the historical context in shaping phenotypic responses to ecological conditions. In particular, we showed that organisms from different geographic locations responded differently to salinity, both in survival and behavior. The significant interaction effect between behavioral category and basin (Table 2) could be attributed to basin-specific environmental parameters not considered in our experimental setup, or to different microevolutionary trajectories of the populations. This highlights the potential role of local adaptation processes in shaping phenotypic and genetic variation in freshwater populations. In Chilean rivers, S. annulinicornis is recognized as a conspicuous species, exhibiting a wide geographical distribution with a clear genetic structure among habitats, studied at the level of allele and haplotype frequencies [31]. This suggests that genetic differentiation may explain the variation of behavioral responses observed in our study organisms across basins. This does not discard the possibility of having a set of closely related species or subspecies still unrecognized. The observed geographically—driven diversity of responses to the instantaneous level of water salinity, reflects that the maintenance of intraspecific diversity –which influence extinction risk in changing environments– rests on the conservation of geographic variation. This appeals for focusing on a landscape—level scale for planning biological freshwater conservation actions. However, our behavioral results should be interpreted with caution. As usual in laboratory experimentation, it is difficult to assess the consistency between the observed behavioral responses under controlled conditions and the animals’ behavior in the field. Although we analyzed the responses of juvenile developmental stages of our model species S. annulicornis, it would be relevant to extend this kind of studies to assess effects also on adult stages, where reproductive outcomes could be measured. Furthermore, and since emerging adult stages are typical prey of several riparian carnivores, studying how water conditions affect adult stages of these invertebrates also allows understanding how the effects of environmental perturbations propagate between adjacent habitats that traditionally have been considered functionally apart from each other, such as river and terrestrial systems. This represents a topic of increasing interest [49-52].

Repeatability and responses of each behavioral trait.

(DOCX)

Click here for additional data file.

8 Nov 2019

PONE-D-19-19484

Geographical origin determines responses to salinity of Mediterranean caddisflies

PLOS ONE

Dear Dr. Ramos-Jiliberto,

Thank you for submitting your manuscript to PLOS ONE. After careful consideration, we feel that it has merit but does not fully meet PLOS ONE’s publication criteria as it currently stands. Therefore, we invite you to submit a minor revision of the manuscript that addresses the points raised during the review process.

In your re-submission, please carefully address each of the reviewer comments. The significant difference in responses between basins was striking, raising the question about the nature of population differentiation. Please address this point in revisions, as well as more minor comments regarding clarification of details.

We would appreciate receiving your revised manuscript by Dec 23 2019 11:59PM. When you are ready to submit your revision, log on to https://www.editorialmanager.com/pone/ and select the 'Submissions Needing Revision' folder to locate your manuscript file.

If you would like to make changes to your financial disclosure, please include your updated statement in your cover letter.

To enhance the reproducibility of your results, we recommend that if applicable you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see: http://journals.plos.org/plosone/s/submission-guidelines#loc-laboratory-protocols

Please include the following items when submitting your revised manuscript:

A rebuttal letter that responds to each point raised by the academic editor and reviewer(s). This letter should be uploaded as separate file and labeled 'Response to Reviewers'.

A marked-up copy of your manuscript that highlights changes made to the original version. This file should be uploaded as separate file and labeled 'Revised Manuscript with Track Changes'.

An unmarked version of your revised paper without tracked changes. This file should be uploaded as separate file and labeled 'Manuscript'.

Please note while forming your response, if your article is accepted, you may have the opportunity to make the peer review history publicly available. The record will include editor decision letters (with reviews) and your responses to reviewer comments. If eligible, we will contact you to opt in or out.

We look forward to receiving your revised manuscript.

Kind regards,

Rachel A Hovel

Academic Editor

PLOS ONE

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. The PLOS ONE style templates can be found at

http://www.journals.plos.org/plosone/s/file?id=wjVg/PLOSOne_formatting_sample_main_body.pdf and http://www.journals.plos.org/plosone/s/file?id=ba62/PLOSOne_formatting_sample_title_authors_affiliations.pdf

1. In your Methods section, please provide additional location information of the study areas, including geographic coordinates for the data set if available.

2. In your Methods section, please provide additional information regarding the permits you obtained for the work. Please ensure you have included the full name of the authority that approved the field site access and, if no permits were required, a brief statement explaining why.

3. We note that you have stated that you will provide repository information for your data at acceptance. Should your manuscript be accepted for publication, we will hold it until you provide the relevant accession numbers or DOIs necessary to access your data. If you wish to make changes to your Data Availability statement, please describe these changes in your cover letter and we will update your Data Availability statement to reflect the information you provide.

[Note: HTML markup is below. Please do not edit.]

Reviewers' comments:

Reviewer's Responses to Questions

Comments to the Author

1. Is the manuscript technically sound, and do the data support the conclusions?

The manuscript must describe a technically sound piece of scientific research with data that supports the conclusions. Experiments must have been conducted rigorously, with appropriate controls, replication, and sample sizes. The conclusions must be drawn appropriately based on the data presented.

Reviewer #1: Yes

Reviewer #2: Yes

**********

2. Has the statistical analysis been performed appropriately and rigorously?

Reviewer #1: Yes

Reviewer #2: Yes

**********

3. Have the authors made all data underlying the findings in their manuscript fully available?

The PLOS Data policy requires authors to make all data underlying the findings described in their manuscript fully available without restriction, with rare exception (please refer to the Data Availability Statement in the manuscript PDF file). The data should be provided as part of the manuscript or its supporting information, or deposited to a public repository. For example, in addition to summary statistics, the data points behind means, medians and variance measures should be available. If there are restrictions on publicly sharing data—e.g. participant privacy or use of data from a third party—those must be specified.

Reviewer #1: No

Reviewer #2: Yes

**********

4. Is the manuscript presented in an intelligible fashion and written in standard English?

PLOS ONE does not copyedit accepted manuscripts, so the language in submitted articles must be clear, correct, and unambiguous. Any typographical or grammatical errors should be corrected at revision, so please note any specific errors here.

Reviewer #1: Yes

Reviewer #2: Yes

**********

5. Review Comments to the Author

Please use the space provided to explain your answers to the questions above. You may also include additional comments for the author, including concerns about dual publication, research ethics, or publication ethics. (Please upload your review as an attachment if it exceeds 20,000 characters)

Reviewer #1: According to the instructions to authors, the authors either need to place the data in a public repository, where URL, accession numbers or DOIs can be cited, or can be contained in Supporting Information files. The current statement says that the data will be available upon request.

The paper needs a little editing, but, per instructions, I did not note specific edits, except that in the second sentence of the Introduction, the verbs should be "disrupts" and "threatens".

The authors might consider discussing how a species might adapt to differences in salinity across its range. At higher taxonomic levels, adaptation to different salinities can involve small mutations to the genes for the ion transporter genes that affect influx and efflux rates. In anadromous or catadromous species, such as salmonids and eels, adaptation may involve up-regulation or down-regulation of these ion transporters or of entire ionocytes.

Reviewer #2: Carter, Flores, Ramos-Jiliberto

Geographical origin determines responses to salinity of Mediterranean caddisflies

This study looks the importance of source (river basin of origin) when assessing behavioral and survivorship responses of larval caddisflies (Smicridea annulicornis) when exposed to elevated salinity. This species appears to be widely distributed across Chile and Argentina, and expresses significant genetic differentiation within a basin, and across basins. Larvae were collected from three basins: Choapa (low conductivity), Maipo (medium conductivity) and Maule (high conductivity), and based on ambient conductivities, they tested three experimental conductivity levels: low (180 μS/cm), medium (500 μS/cm) and high (1400 μS/cm) different salinity levels were prepared with aquarium sea salt dissolved in distilled water. Response variables measured were time to death, and behavior (swimming, sheltering, walking, pushing up). This assessment of behavioral responses (i.e., non-lethal) was an interesting addition to the standard survival study. The experiments found significant differences between basins – big differences. This is most unusual and suggests measurable local adaptation or acclimation.

One factor that is only lightly addressed is this species in this region exhibits significant genetic differentiation among basins (Sabando et al. 2011) – suggesting either local adaptation with little or no gene exchange, or possibly some closely related species or subspecies that are unrecognized and responding to a strong habitat template.

Specific comments;

I assume the experiments were conducted at 20 C, same as the acclimation temperature.

Survival time – Table 1 indicates that basin matters, and that the effect of conductivity was not significant (p=0.08), but the results present a bunch of antidotal assessements of conductivity effects. Table 1 results and the Survival results need to be in agreement – at least justify the details by calling the factorial result “nearly significant” and not “marginal”.

In the behavioral tests, why do you think the BC x Basin was significant? Why would behavioral responses differ between basins within a species?

In the behavioral tests, do you believe behavioral responses are in some what adaptive with regard to salinity?

The results presented here make one question standard bioassay results – where was the original test population from? Have laboratory conditions changed the test organisms?

**********

6. PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: No

Reviewer #2: No

[NOTE: If reviewer comments were submitted as an attachment file, they will be attached to this email and accessible via the submission site. Please log into your account, locate the manuscript record, and check for the action link "View Attachments". If this link does not appear, there are no attachment files to be viewed.]

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool, https://pacev2.apexcovantage.com/. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Registration is free. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org. Please note that Supporting Information files do not need this step.

28 Nov 2019

PONE-D-19-19484

Geographical origin determines responses to salinity of Mediterranean caddisflies

PLOS ONE

RESPONSE TO REVIEWERS (POINT-BY-POINT)

Journal Requirements:

When submitting your revision, we need you to address these additional requirements.

Please ensure that your manuscript meets PLOS ONE's style requirements, including those for file naming. Response: Checked

1. In your Methods section, please provide additional location information of the study areas, including geographic coordinates for the data set if available.

Response: Done as requested. L 79-80

2. In your Methods section, please provide additional information regarding the permits you obtained for the work. Please ensure you have included the full name of the authority that approved the field site access and, if no permits were required, a brief statement explaining why.

Response: No permits are required, according to national regulations. We included this information in the text. L 80-82

3. We note that you have stated that you will provide repository information for your data at acceptance. Should your manuscript be accepted for publication, we will hold it until you provide the relevant accession numbers or DOIs necessary to access your data. If you wish to make changes to your Data Availability statement, please describe these changes in your cover letter and we will update your Data Availability statement to reflect the information you provide.

Response: We uploaded the data to the dryad repository. We included this information in the cover letter as requested.

Review Comments to the Author

Reviewer #1:

According to the instructions to authors, the authors either need to place the data in a public repository, where URL, accession numbers or DOIs can be cited, or can be contained in Supporting Information files. The current statement says that the data will be available upon request.

Response: We uploaded the data to the dryad repository. We included this information in the cover letter as requested.

The paper needs a little editing, but, per instructions, I did not note specific edits, except that in the second sentence of the Introduction, the verbs should be "disrupts" and "threatens".

Response: Corrected. L 4-5

The authors might consider discussing how a species might adapt to differences in salinity across its range. At higher taxonomic levels, adaptation to different salinities can involve small mutations to the genes for the ion transporter genes that affect influx and efflux rates. In anadromous or catadromous species, such as salmonids and eels, adaptation may involve up-regulation or down-regulation of these ion transporters or of entire ionocytes.

Response: Corrected. In the 2nd paragraph of the Discussion section we now talk about rapid evolution and mechanisms of ion regulation in aquatic organisms, including macroinvertebrates. L 210-227

Reviewer #2:

One factor that is only lightly addressed is this species in this region exhibits significant genetic differentiation among basins (Sabando et al. 2011) – suggesting either local adaptation with little or no gene exchange, or possibly some closely related species or subspecies that are unrecognized and responding to a strong habitat template.

Response: We extend our discussion about the empirical evidence of genetic structure in our study species on the study sites, and also include the reviewer’s idea of unrecognized species or subspecies. L 237-242

Specific comments;

I assume the experiments were conducted at 20 C, same as the acclimation temperature.

Response: Yes. This information was included in the revised version. L 102-104

Survival time – Table 1 indicates that basin matters, and that the effect of conductivity was not significant (p=0.08), but the results present a bunch of antidotal assessments of conductivity effects. Table 1 results and the Survival results need to be in agreement – at least justify the details by calling the factorial result “nearly significant” and not “marginal”.

Response: We corrected the text following the reviewer's suggestions. L 148

In the behavioral tests, why do you think the BC x Basin was significant? Why would behavioral responses differ between basins within a species?

Response: Done. We added a phrase in which we –with caution-- mention possible causes. L 232-235

In the behavioral tests, do you believe behavioral responses are in some what adaptive with regard to salinity?

Response: Our results do not allow inferring that behavioral responses are adaptive respect to salinity. We guess that the previous response covers this aspect too.

The results presented here make one question standard bioassay results – where was the original test population from? Have laboratory conditions changed the test organisms?

Response: Corrected. We commented on the limitation of lab experimentation. 250-253

Submitted filename: Response to Reviewers.docx

Click here for additional data file.

13 Dec 2019

Geographical origin determines responses to salinity of Mediterranean caddisflies

PONE-D-19-19484R1

Dear Dr. Ramos-Jiliberto,

We are pleased to inform you that your manuscript has been judged scientifically suitable for publication and will be formally accepted for publication once it complies with all outstanding technical requirements.

Within one week, you will receive an e-mail containing information on the amendments required prior to publication. When all required modifications have been addressed, you will receive a formal acceptance letter and your manuscript will proceed to our production department and be scheduled for publication.

Shortly after the formal acceptance letter is sent, an invoice for payment will follow. To ensure an efficient production and billing process, please log into Editorial Manager at https://www.editorialmanager.com/pone/, click the "Update My Information" link at the top of the page, and update your user information. If you have any billing related questions, please contact our Author Billing department directly at authorbilling@plos.org.

If your institution or institutions have a press office, please notify them about your upcoming paper to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, you must inform our press team as soon as possible and no later than 48 hours after receiving the formal acceptance. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information, please contact onepress@plos.org.

With kind regards,

Rachel A Hovel

Academic Editor

PLOS ONE

Additional Editor Comments (optional):

This revision has sufficiently addressed reviewer comments and outstanding concerns on this manuscript, and I am pleased to accept it. Please make sure to thoroughly proofread the manuscript before submission. For example, Line 71 now reads: "We conducted this study usin as a model speciesgthe...". It will be your responsibility to identify and remedy typographical or other errors.

Reviewers' comments:

19 Dec 2019

PONE-D-19-19484R1

Geographical origin determines responses to salinity of Mediterranean caddisflies

Dear Dr. Ramos-Jiliberto:

I am pleased to inform you that your manuscript has been deemed suitable for publication in PLOS ONE. Congratulations! Your manuscript is now with our production department.

If your institution or institutions have a press office, please notify them about your upcoming paper at this point, to enable them to help maximize its impact. If they will be preparing press materials for this manuscript, please inform our press team within the next 48 hours. Your manuscript will remain under strict press embargo until 2 pm Eastern Time on the date of publication. For more information please contact onepress@plos.org.

For any other questions or concerns, please email plosone@plos.org.

Thank you for submitting your work to PLOS ONE.

With kind regards,

PLOS ONE Editorial Office Staff

on behalf of

Dr. Rachel A Hovel

Academic Editor

PLOS ONE