Phylogenetic signal in diatom ecology: perspectives for aquatic ecosystems biomonitoring

Microbial community structure and diversity attached to the periphyton in different urban aquatic habitats

Periphyton is a complex community composed of diverse prokaryotes and eukaryotes; understanding the characteristics of microbial communities within periphyton becomes crucial for biogeochemical cycles and energy dynamics of aquatic ecosystems. To further elucidate the community characteristics of periphyton across varied aquatic habitats, including unpolluted ecologically restored lakes, aquaculture ponds, and areas adjacent to domestic and industrial wastewater treatment plant outfalls, we explored the composition and diversity of prokaryotic and eukaryotic communities in periphyton by employing Illumina MiSeq sequencing. Our findings indicated that the prokaryotic communities were predominantly composed of Proteobacteria (40.92%), Bacteroidota (21.01%), and Cyanobacteria (10.12%), whereas the eukaryotic communities were primarily characterized by the dominance of Bacillariophyta (24.09%), Chlorophyta (20.83%), and Annelida (15.31%). Notably, Flavobacterium emerged as a widely distributed genus among the prokaryotic community. Unclassified_Tobrilidae exhibited higher abundance in unpolluted ecologically restored lakes. Chaetogaster and Nais were enriched in aquaculture ponds and domestic wastewater treatment plant outfall area, respectively, while Surirella and Gomphonema dominated industrial sewage treatment plant outfall area. The alpha diversity of eukaryotes was higher in unpolluted ecologically restored lakes. pH and nitrogen content (NO 2 - - N ,NO 3 - - N , and TN) significantly explained the variations for prokaryotic and eukaryotic community structures, respectively. Eukaryotic communities exhibited a more pronounced response to habitat variations compared to prokaryotic communities. Moreover, the association networks revealed an intensive positive correlation between dominant Bacillariophyta and Bacteroidota. This study provided useful data for identifying keystone species and understanding their ecological functions.

Linking niche size and phylogenetic signals to predict future soil microbial relative abundances

Bacteria provide ecosystem services (e.g., biogeochemical cycling) that regulate climate, purify water, and produce food and other commodities, yet their distribution and likely responses to change or intervention are difficult to predict. Using bacterial 16S rRNA gene surveys of 1,381 soil samples from the Biomes of Australian Soil Environment (BASE) dataset, we were able to model relative abundances of soil bacterial taxonomic groups and describe bacterial niche space and optima. Hold out sample validated hypothetical causal networks (structural equation models; SEM) were able to predict the relative abundances of bacterial taxa from environmental data and elucidate soil bacterial niche space. By using explanatory SEM properties as indicators of microbial traits, we successfully predicted soil bacterial response, and in turn potential ecosystem service response, to near-term expected changes in the Australian climate. The methods developed enable prediction of continental-scale changes in bacterial relative abundances, and demonstrate their utility in predicting changes in bacterial function and thereby ecosystem services. These capabilities will be strengthened in the future with growing genome-level data.

Temporal pesticide dynamics alter specific eukaryotic taxa in a coastal transition zone

Land use change and anthropogenic forcing can drastically alter the rates and patterns of sediment transport and modify biodiversity and ecosystem functions in coastal transition zones, such as the coastal ecosystems. Molecular studies of sediment extracted DNAs provide information on currently living organisms within the upper layers or buried from various periods of time, but might also provide knowledge on species dynamics, replacement and turnover. In this study, we evaluated the eukaryotic communities of a marine core that present a shift in soil erosion that was linked to glyphosate usage and correlated to chlordecone resurgence since 2000. We show differences in community composition between samples from the second half of the last century and those from the last two decades. Temporal analyses of the relative abundance, alpha diversity, and beta diversity for the two periods demonstrated different temporal dynamics depending on the considered taxonomic group. In particular, Ascomycetes showed a decrease in abundance over the most recent period associated with changes in community membership but not community structure. Two photosynthetic groups, Bacillariophyceae and Prasinophytes clade VII, showed a different pattern with an increase in abundance since the beginning of the 21st century with a decrease in diversity and evenness to form more heterogeneous communities dominated by a few abundant OTUs. Altogether, our data reveal that agricultural usages such as pesticide use can have long-term and species-dependent implications for microeukaryotic coastal communities on a tropical island.

Genus-level, trait-based multimetric diatom indices for assessing the ecological condition of rivers and streams across the conterminous United States

Taxonomic inconsistency in species-level identifications has constrained use of diatoms as biological indicators in aquatic assessments. We addressed this problem by developing diatom multimetric indices (MMIs) of ecological condition using genus-level taxonomy and trait-based autecological information. The MMIs were designed to assess river and stream chemical, physical and biological condition across the conterminous United States. Trait-based approaches have the advantage of using both species-level and genus-level data, which require less effort and expense to acquire than traditional species-based approaches and eliminate the persistent taxonomic biases introduced over vast geographic extents. For large-extent assessment programs that require multiple taxonomic laboratories to process samples, such as the United States Environmental Protection Agency's (U.S. EPA's) National Rivers and Streams Assessment (NRSA), the trait approach can eliminate discrepancies in species-level identification or nomenclature that hinder diatom data interpretation. We developed trait-based MMIs using NRSA data for each of the three large ecoregions across the U.S. - the East, Plains, and West. All three MMIs performed well in discriminating least-disturbed from most-disturbed sites. The MMI for the East had the greatest discrimination ability, followed by MMIs for the Plains and West, respectively. The performance of the MMIs was comparable to that observed in existing NRSA fish and macroinvertebrate MMIs. Our research shows that trait-based diatom indices constructed on genus-level taxonomy can be effective for large-scale assessments, and may also allow programs such as NRSA to assess trends in freshwater condition retrospectively, by revisiting older diatom datasets. Moreover, our genus-based approach facilitates including of diatoms into other assessment programs that have limited monitoring resources.

Wastewater constituents impact biofilm microbial community in receiving streams

Microbial life in natural biofilms is dominated by prokaryotes and microscopic eukaryotes living in dense association. In stream ecosystems, microbial biofilms influence primary production, elemental cycles, food web interactions as well as water quality. Understanding how biofilm communities respond to anthropogenic impacts, such as wastewater treatment plant (WWTP) effluent, is important given the key role of biofilms in stream ecosystem function. Here, we implemented 16S and 18S rRNA gene sequencing of stream biofilms upstream (US) and downstream (DS) of WWTP effluents in four Swiss streams to test how bacterial and eukaryotic communities respond to wastewater constituents. Stream biofilm composition was strongly affected by geographic location - particularly for bacteria. However, the abundance of certain microbial community members was related to micropollutants in the wastewater - among bacteria, micropollutant-associated members were found e.g. in Alphaproteobacteria, and among eukaryotes e.g. in Bacillariophyta (algal diatoms). This study corroborates several previously characterized responses (e.g. as seen in diatoms), but also reveals previously unknown community responses - such as seen in Alphaproteobacteria. This study advances our understanding of the ecological impact of the current wastewater treatment practices and provides information about potential new marker organisms to assess ecological change in stream biofilms.

DNA metabarcoding reveals differences in distribution patterns and ecological preferences among genetic variants within some key freshwater diatom species

Our study evaluates differences in the distribution and ecology of genetic variants within several ecologically important diatom species that are also key for Water Framework Directive monitoring of European rivers: Fistulifera saprophila (FSAP), Achnanthidium minutissimum (ADMI), Nitzschia inconspicua (NINC) and Nitzschia soratensis (NSTS). We used DADA2 to infer amplicon sequence variants (ASVs) of a short rbcL barcode in 531 environmental samples from biomonitoring campaigns in Catalonia and France. ASVs within each species showed different distribution patterns. Threshold Indicator Taxa ANalysis revealed three ecological groupings of ASVs in both ADMI and FSAP. Two of these in each species were separated by opposite responses to calcium and conductivity. Boosted regression trees additionally showed that both variables greatly influenced the occurrence of these groupings. A third grouping in FSAP was characterized by a negative response to total organic carbon and hence was better represented in waters with higher ecological status than the other FSAP ASVs, contrasting with what is generally assumed for the species. In the two Nitzschia species, our analyses confirmed earlier studies: NINC preferred higher levels of calcium and conductivity. Our findings suggest that the broad ecological tolerance of some diatom species results from overlapping preferences among genetic variants, which individually show much more restricted preferences and distributions. This work shows the importance of studying the ecological preferences of genetic variants within species complexes, now possible with DNA metabarcoding. The results will help reveal and understand biogeographical distributions and facilitate the development of more accurate biological indexes for biomonitoring programmes.

Dumpster diving for diatom plastid 16S rRNA genes

High throughput sequencing is improving the efficiency of monitoring diatoms, which inhabit and support aquatic ecosystems across the globe. In this study, we explored the potential of a standard V4 515F-806RB primer pair in recovering diatom plastid 16S rRNA sequences. We used PhytoREF to classify the 16S reads from our freshwater biofilm field sampling from three stream segments across two streams in south-eastern Australia and retrieved diatom community data from other, publicly deposited, Australian 16S amplicon datasets. When these diatom operational taxonomic units (OTUs) were traced using the default RDPII and NCBI databases, 68% were characterized as uncultured cyanobacteria. We analysed the 16S rRNA sequences from 72 stream biofilm samples, separated the chloroplast OTUs, and classified them using the PhytoREF database. After filtering the reads attributed to Bacillariophyta (relative abundance >1%), 71 diatom OTUs comprising more than 90% of the diatom reads in each stream biofilm sample were identified. Beta-diversity analyses demonstrated significantly different diatom assemblages and discrimination among river segments. To further test the approach, the diatom OTUs from our biofilm sampling were used as reference sequences to identify diatom reads from other Australian 16S rRNA datasets in the NCBI-SRA database. Across the three selected public datasets, 67 of our 71 diatom OTUs were detected in other Australian ecosystems. Our results show that diatom plastid 16S rRNA genes are readily amplified with existing 515F-806RB primer sets. Therefore, the volume of existing 16S rRNA amplicon datasets initially generated for microbial community profiling can also be used to detect, characterize, and map diatom distribution to inform phylogeny and ecological health assessments, and can be extended into a range of ecological and industrial applications. To our knowledge, this study represents the first attempt to classify freshwater samples using this approach and the first application of PhytoREF in Australia.

Ciliates as bioindicators of water quality: A case study in the neotropical region and evidence of phylogenetic signals (18S-rDNA)

The aim of our study was to evaluate the water quality of an urban stream in southeastern Brazil by analyzing epibenthic ciliates, and to investigate the existence of phylogenetic signal for saprobity in ciliates. However, before conducting this type of phylogenetic study, it is necessary to evaluate if the saprobic classification used frequently in the northern Hemisphere is suitable for neotropical ecosystems. Sediment samples were collected from five sampling stations: two in rural areas and three in urban areas. During the one-year study, with monthly collections, 39 ciliates species were found, of which 32 are included in the saprobic system. Physical, chemical and biological parameters of water confirm the spatial heterogeneity of the sampling stations, with a clear influence of organic pollution on the composition and structure of ciliates taxocenosis. The saprobic index and the saprobic valence index were used to evaluate the water quality of the sampling stations and demonstrated clear heterogeneity between the stations and high degree of pollution of the urban area. These sampling stations were dominated by ciliates indicators of polysaprobric environments. Since we were able to successfully use the saprobic index in a limnic ecosystem in Brazil, we applied the phylogenetic signal validation as a tool for saprobity prediction of the limnic ciliate species not yet analyzed. A phylogenetic tree containing only 18S-rDNA nominal sequences of freshwater ciliates was estimated and used to explore the existence of the phylogenetic signal, which showed that the sensitivity/tolerance of ciliates to organic pollution reflected evolutionary divergence. The results confirm the existence of phylogenetic signal for the saprobrity in Ciliophora. Also, our results suggest that evolutionary analysis is a potential method to predict lineages of ciliates not yet classified for saprobity.

Trait-based diatom functional diversity as an appropriate tool for understanding the effects of environmental changes in soda pans

Saline lakes, among the most seriously endangered ecosystems, are threatened due to climate change and human activities. One valuable feature of these environments is that they constitute areas of high biodiversity. Ecologists are, therefore, under great pressure to improve their understanding of the effects of natural and anthropogenic disturbances on the biodiversity of saline lakes. In this study, a total of 257 samples from 32 soda pans in Central Europe between 2006 and 2015 were examined. The effects of environmental variables and of geographical and limnoecological factors on functional diversity were analyzed. Furthermore, the explanatory power of the trait-based approach was assessed, and the applicability of the indices for biomonitoring purposes was determined. It was found that low habitat heterogeneity and harsh environments lead to the selection of a small number of suitable traits, and consequently, to a naturally low level of functional diversity. Anthropogenic activities enhance diversity at functional level due to the shift toward freshwater characteristics. On the regional scale, the effects of the region and status (natural, degraded, reconstructed) on diatom functional diversity were significant and more pronounced than that of the environmental and other limnoecological factors. The degree of variance found in functional diversity ascribed to environmental variables is five times greater in the case of the application of a trait-based approach, than when a taxonomic one is employed in the literature. Each of the tested functional diversity indices was sensitive to the most important environmental variables. Furthermore, these were type-specific and proved to be more complex indicators than taxonomic metrics. It is possible to suggest four functional diversity indices (FGR, FRic, FDis, and FDiv) which emphasize their independence from substrate and seasonal variations for ecological status assessment and conservation planning.

Community phylogenetic structure reveals the imprint of dispersal-related dynamics and environmental filtering by nutrient availability in freshwater diatoms

Despite important progress, uncertainty persists regarding the ecological forces driving microbial community assembly. Here, we present the first study to use phylogenetic information to interpret the structure and diversity of diatom communities. We examined local phylogenetic divergence and beta- phylogenetic diversity in a large dataset of 595 freshwater benthic diatom communities and we investigated how this diversity is influenced by gradients in nutrients, pH, organic matter and catchment size. Overall, we found that diatom communities were phylogenetically clustered, i.e. species within communities were more closely related than expected by chance. Phylogenetic clustering was stronger in nutrient-poor environments and in sites with a small catchment area. The variation of the phylogenetic beta-diversity index was much better explained by space and environment than the variation of the taxonomic index was. Both approaches detected a significant effect of environment and space on diatom community turnover. Our results support the view that diatom communities are primarily shaped by environmental filtering, in particular by nutrient availability. Moreover, they highlight the importance of considering dispersal-related processes and the depth of phylogenetic signal in functional traits when interpreting patterns of diversity.

Exploring taxonomic and phylogenetic relationships to predict radiocaesium transfer to marine biota

One potentially useful approach to fill data gaps for concentration ratios, CRs, is based upon the hypothesis that an underlying taxonomic and/or phylogenetic relationship exists for radionuclide transfer. The objective of this study was to explore whether these relationships could be used to explain variation in the transfer of radiocaesium to a wide range of marine organisms. CR data for ¹³⁷Cs were classified in relation to taxonomic family, order, class and phylum. A Residual Maximum Likelihood (REML) mixed-model regression modelling approach was adopted. The existence of any patterns were then explored using phylogenetic trees constructed with mitochondrial COI gene sequences from various biota groups and mapping the REML residual means onto these trees. A comparison of the predictions made using REML with blind datasets allowed the efficacy of the procedure to be tested. The only significant correlation between predicted and measured activity concentrations was revealed at the taxonomic level of order when comparing REML analysis output with data from the Barents Sea Region. For this single case a correlation 0.80 (Spearman rank) was derived which was significant at the 0.01 level (1-tailed test) although this was not the case once a (Bonferroni) correction was applied. The application of the REML approach to marine datasets has met with limited success, and the phylogenetic trees illustrate complications of using predictions based on values from different levels of taxonomic organization, where predicted values for the order level can mask the values at lower taxonomic levels. Any influence of taxonomy and phylogeny on transfer is not immediately conspicuous and categorizing marine organisms in this way is limited in providing a potentially robust prognostic extrapolation tool. Other factors may plausibly affect transfer to a much greater degree in marine systems, such as quite diverse life histories and different diets, which may confound any phylogenetic pattern.

Sensitivities of seven algal species to triclosan, fluoxetine and their mixtures

Increasing release of pharmaceuticals and personal care products (PPCPs) into aquatic ecosystems is a growing environmental concern. Triclosan and fluoxetine are two widely used PPCPs and frequently detected in aquatic ecosystems. In this study, the sensitivities of 7 algal species from 4 genera to triclosan, fluoxetine and their mixture were evaluated. The results showed that the inhibitory effect on algal growth (EC₅₀-96h) of triclosan varied with 50 times differences among the 7 algal species. Chlorella ellipsoidea was the least susceptible species and Dunaliella parva was the most sensitive species to triclosan. The inhibitory effect of fluoxetine was less variable than triclosan. Slightly higher toxicity of fluoxetine than triclosan was shown in the 7 tested algal species. No consistent pattern of the effects from mixture of triclosan and fluoxetine was observed among the 7 algal species and among the 4 genera. Additive effects of the mixture occured in 4 species and antagonistic effects in the other 3 species but no synergistic effect was detected. The algal species might show some sign of phylogenetic response to triclosan, as evidenced by the wide range of differences in their sensitivity at the genus level. This study provides important data which could be beneficial for biomonitoring programs on the ecological risk (algal species diversity) of these two chemicals.

The future of biotic indices in the ecogenomic era: Integrating (e)DNA metabarcoding in biological assessment of aquatic ecosystems

The bioassessment of aquatic ecosystems is currently based on various biotic indices that use the occurrence and/or abundance of selected taxonomic groups to define ecological status. These conventional indices have some limitations, often related to difficulties in morphological identification of bioindicator taxa. Recent development of DNA barcoding and metabarcoding could potentially alleviate some of these limitations, by using DNA sequences instead of morphology to identify organisms and to characterize a given ecosystem. In this paper, we review the structure of conventional biotic indices, and we present the results of pilot metabarcoding studies using environmental DNA to infer biotic indices. We discuss the main advantages and pitfalls of metabarcoding approaches to assess parameters such as richness, abundance, taxonomic composition and species ecological values, to be used for calculation of biotic indices. We present some future developments to fully exploit the potential of metabarcoding data and improve the accuracy and precision of their analysis. We also propose some recommendations for the future integration of DNA metabarcoding to routine biomonitoring programs.

Boosting DNA metabarcoding for biomonitoring with phylogenetic estimation of operational taxonomic units' ecological profiles

DNA metabarcoding has been introduced as a revolutionary way to identify organisms and monitor ecosystems. However, the potential of this approach for biomonitoring remains partially unfulfilled because a significant part of the sampled DNA cannot be affiliated to species due to incomplete reference libraries. Thus, biotic indices, which are based on the estimated abundances of species in a community and their ecological profiles, can be inaccurate. We propose to compute biotic indices using phylogenetic imputation of operational taxonomic units (OTUs') ecological profiles (OTU-PITI approach). First, OTUs sequences are inserted within a reference phylogeny. Second, OTUs' ecological profiles are estimated on the basis of their phylogenetic relationships with reference species whose ecology is known. Based on these ecological profiles, biotic indices can be computed using all available OTUs. Using freshwater diatoms as a case study, we show that short DNA barcodes can be placed accurately within a phylogeny and their ecological preferences estimated with a satisfactory level of precision. In the light of these results, we tested the approach with a data set of 139 environmental samples of benthic river diatoms for which the same biotic index (specific sensitivity index) was calculated using (a) traditional microscopy, (b) OTUs with taxonomic assignment approach, (c) OTUs with phylogenetic estimation of ecological profiles (OTU-PITI) and (d) OTU with taxonomic assignment completed by the phylogenetic approach (OTU-PITI) for unclassified OTUs. Using traditional microscopy as a reference, we found that the combination of the OTUs' taxonomic assignment completed by the phylogenetic method performed satisfactorily and substantially better than the other methods tested.

DNA metabarcoding and microscopic analyses of sea turtles biofilms: Complementary to understand turtle behavior

Sea turtles are distributed in tropical and subtropical seas worldwide. They play several ecological roles and are considered important indicators of the health of marine ecosystems. Studying epibiotic diatoms living on turtle shells suggestively has great potential in the study of turtle behavior because diatoms are always there. However, diatom identification at the species level is time consuming, requires well-trained specialists, and there is a high probability of finding new taxa growing on turtle shells, which makes identification tricky. An alternative approach based on DNA barcoding and high throughput sequencing (HTS), metabarcoding, has been developed in recent years to identify species at the community level by using a DNA reference library. The suitabilities of morphological and molecular approaches were compared. Diatom assemblages were sampled from seven juvenile green turtles (Chelonia mydas) from Mayotte Island, France. The structures of the epibiotic diatom assemblages differed between both approaches. This resulted in different clustering of the turtles based on their diatom communities. Metabarcoding allowed better discrimination between turtles based on their epibiotic diatom assemblages and put into evidence the presence of a cryptic diatom diversity. Microscopy, for its part, provided more ecological information of sea turtles based on historical bibliographical data and the abundances of ecological guilds of the diatom species present in the samples. This study shows the complementary nature of these two methods for studying turtle behavior.

Species complexes within epiphytic diatoms and their relevance for the bioindication of trophic status

The popularity of aquatic bioassessments has increased in Europe and worldwide, with a considerable number of methods being based on benthic diatoms. Recent evidence from molecular data and mating experiments has shown that some traditional diatom morphospecies represent species complexes, containing several to many cryptic species. This case study is based on epiphytic diatom and environmental data from shallow fishponds, investigating whether the recognition and use of fine taxonomic resolution (cryptic species) can improve assessment of community response to environmental drivers and increase sharpness of classification, compared to coarse taxonomic resolution (genus level and species level with unresolved species complexes). Secondly, trophy bioindication based on a species matrix divided into two compartments (species complexes and remaining species) was evaluated against the expectation that species complexes would be poor trophy indicators, due to their expected wide ecological amplitude. Finally, the response of species complexes and their members (cryptic species) to a trophic gradient (phosphorus) were compared. Multivariate analyses showed similar efficiency of all three taxonomic resolutions in depicting community patterns and their environmental correlates, suggesting that even genus level resolution is sufficient for routine bioassessment of shallow fishponds with a wide trophic range. However, after controlling for coarse taxonomic matrices, fine taxonomic resolution (with resolved cryptic species) still showed sufficient variance related to the environmental variable (habitat groups), and increased the sharpness of classification, number of indicator species for habitat categories, and gave better separation of habitat categories in the ordination space. Regression analysis of trophic bioindication and phosphorus concentration showed a weak relationship for species complexes but a close relationship for the remaining taxa. GLM models also showed that no species complex responded to phosphorus concentration. It follows that the studied species complexes have wide tolerances to, and no apparent optima for, phosphorus concentrations. In contrast, various responses (linear, unimodal, or no response) of cryptic species within species complexes were found to total phosphorus concentration. In some cases, fine taxonomic resolution to species level including cryptic species has the potential to improve data interpretation and extrapolation, supporting recent views of species surrogacy.

Can we predict diatoms herbicide sensitivities with phylogeny? Influence of intraspecific and interspecific variability

Diatoms are used as indicators of freshwater ecosystems integrity. Developing diatom-based tools to assess impact of herbicide pollution is expected by water managers. But, defining sensitivities of all species to multiple herbicides would be unattainable. The existence of a phylogenetic signal of herbicide sensitivity was shown among diatoms and should enable prediction of new species sensitivity. However, diatoms present a cryptic diversity that may lead to variation in their sensitivity to herbicides that would need to be taken into account. Using bioassays, the sensitivity to four herbicides (Atrazine, Terbutryn, Diuron, Isoproturon) was evaluated for 11 freshwater diatom taxa and intraspecific variability was assessed for two of them (Nitzschia palea and Achnanthidium spp.). Intraspecific variability of herbicide sensitivity was always smaller than interspecific variability, but intraspecific variability was more important in N. palea than in Achnanthidium spp. Indeed, one species showed no intraspecific phylogenetic signal (N. palea) whereas the other did (Achnanthidium spp.). On one hand, species boundaries are not set properly for Achnanthidium spp. which encompass several taxa. On the other hand, there is a higher phenotypic plasticity for N. palea. Finally, a phylogenetic signal of herbicide sensitivity was measured at the interspecific level, opening up prospects for setting up reliable biomonitoring tools based on sensitivity prediction, insofar as species boundaries are correctly defined.

Integrative behavioral ecotoxicology: bringing together fields to establish new insight to behavioral ecology, toxicology, and conservation

The fields of behavioral ecology, conservation science, and environmental toxicology individually aim to protect and manage the conservation of wildlife in response to anthropogenic stressors, including widespread anthropogenic pollution. Although great emphasis in the field of toxicology has been placed on understanding how single pollutants affect survival, a comprehensive, interdisciplinary approach that includes behavioral ecology is essential to address how anthropogenic compounds are a risk for the survival of species and populations in an increasingly polluted world. We provide an integrative framework for behavioral ecotoxicology using Tinbergen’s four postulates (causation and mechanism, development and ontogeny, function and fitness, and evolutionary history and phylogenetic patterns). The aims of this review are: 1) to promote an integrative view and re-define the field of integrative behavioral ecotoxicology; 2) to demonstrate how studying ecotoxicology can promote behavior research; and 3) to identify areas of behavioral ecotoxicology that require further attention to promote the integration and growth of the field.

phylosignal: an R package to measure, test, and explore the phylogenetic signal

Phylogenetic signal is the tendency for closely related species to display similar trait values as a consequence of their phylogenetic proximity. Ecologists and evolutionary biologists are becoming increasingly interested in studying the phylogenetic signal and the processes which drive patterns of trait values in the phylogeny. Here, we present a new R package, phylosignal which provides a collection of tools to explore the phylogenetic signal for continuous biological traits. These tools are mainly based on the concept of autocorrelation and have been first developed in the field of spatial statistics. To illustrate the use of the package, we analyze the phylogenetic signal in pollution sensitivity for 17 species of diatoms.