DOI: 10.5061/dryad.7s11m96
title: Data from: Differential response to heat-stress among evolutionary lineages of an aquatic invertebrate species complex

description: Under global warming scenarios, rising temperatures can constitute heat
 stress to which species may differentially respond. Within a described
 species, knowledge on cryptic diversity is of further relevance, as
 different lineages/cryptic species may differentially respond to
 environmental change. The Brachionus calyciflorus species complex
 (Rotifera), which was recently described using integrative taxonomy, is an
 essential component of aquatic ecosystems. Here, we tested the hypothesis
 that these (formerly cryptic) species differ in their heat tolerance. We
 assigned 47 clones with nuITS1 and mtCOI markers to evolutionary lineages,
 now named Brachionus calyciflorus sensu stricto and Brachionus fernandoi.
 We selected 15 representative clones, and assessed their heat tolerance as
 a bi-dimensional phenotypic trait affected by both intensity and duration
 of heat-stress. We found two distinct groups, with B. calyciflorus sensu
 stricto clones having higher heat tolerance than the novel species B.
 fernandoi. This apparent temperature specialization among former cryptic
 species underscores the necessity of a sound species delimitation and
 assignment, when organismal responses to environmental change are
 investigated.
description: alignmentscontains all the alignments of data presented either in the main text or the supplementary materialtreesInput and output files for BEAST v.1.8.1Thermal_Death_Time_Curves_pipeline_REV2Contains an excel file with the pipeline needed to create the Thermal Death Time curves and extrapolate the CTmax values. The pipeline is separated in 7 steps. Step1 contains the observed data (time in minutes and percentage (%) of incapacitated individuals). Step 2 contains the log10 transformed data of step 1. Step 3 represents the fitting of linear regressions to the data of step 2. Step 4 contains the estimation of time needed for LT50 from the regressions presented in step 3. In Step 5 we plotted the time needed for LT50 (in logarithmic scale) against the experimental temperature to create the TDT curves for each clone separately. In step 6 we estimated the CTmax and the z value from the intercept and slope of the TDT curves. In step 7 we used the estimated values in step 6 to create the thermal tolerance landscape.Generalized Linear Mixed ModelsContains the data input in R software we used to run GLMM models to test the effect of temperature and species on the response variable, i.e., the time when 50% of the experimental population was incapacitated (LT50). Lineage c represents the species Brachionus calyciflorus s.s. and lineage d represents the species Brachionus fernandoi. Also contains the code for running GLMM in R using lmer4 and car packages and a file with the console output.GLMM.rart-test for CTmax valuesContains the CTmax input data used to test for normal distribution and t-test. Also contains the R script used for both Shapiro-Wilk and t-test as well as the output of the previous analyses on the console. CTmax represents Critical Themal Maximun values estimated as an extrapolation from the Thermat Death Time (TDT) curves. Lineage c represents species Brachinous calyciforus s.s. and lineage d represents the species Brachionus fernandoi.t-test.rar