Integrative phylogenetic, phylogeographic and morphological characterisation of the Unio crassus species complex reveals cryptic diversity with important conservation implications

The global decline of freshwater mussels and their crucial ecological services highlight the need to understand their phylogeny, phylogeography and patterns of genetic diversity to guide conservation efforts. Such knowledge is urgently needed for Unio crassus, a highly imperilled species originally widespread throughout Europe and southwest Asia. Recent studies have resurrected several species from synonymy based on mitochondrial data, revealing U. crassus to be a complex of cryptic species. To address long-standing taxonomic uncertainties hindering effective conservation, we integrate morphometric, phylogenetic, and phylogeographic analyses to examine species diversity within the U. crassus complex across its entire range. Phylogenetic analyses were performed using cytochrome c oxidase subunit I (815 specimens from 182 populations) and, for selected specimens, whole mitogenome sequences and Anchored Hybrid Enrichment (AHE) data on ∼ 600 nuclear loci. Mito-nuclear discordance was detected, consistent with mitochondrial DNA gene flow between some species during the Pliocene and Pleistocene. Fossil-calibrated phylogenies based on AHE data support a Mediterranean origin for the U. crassus complex in the Early Miocene. The results of our integrative approach support 12 species in the group: the previously recognised Unio bruguierianus, Unio carneus, Unio crassus, Unio damascensis, Unio ionicus, Unio sesirmensis, and Unio tumidiformis, and the reinstatement of five nominal taxa: Unio desectusstat. rev., Unio gontieriistat. rev., Unio mardinensisstat. rev., Unio nanusstat. rev., and Unio vicariusstat. rev. Morphometric analyses of shell contours reveal important morphospace overlaps among these species, highlighting cryptic, but geographically structured, diversity. The distribution, taxonomy, phylogeography, and conservation of each species are succinctly described.

Deconstructing Difflugia: The tangled evolution of lobose testate amoebae shells (Amoebozoa: Arcellinida) illustrates the importance of convergent evolution in protist phylogeny

Protists, the micro-eukaryotes that are neither plants, animals nor fungi build up the greatest part of eukaryotic diversity on Earth. Yet, their evolutionary histories and patterns are still mostly ignored, and their complexity overlooked. Protists are often assumed to keep stable morphologies for long periods of time (morphological stasis). In this work, we test this paradigm taking Arcellinida testate amoebae as a model. We build a taxon-rich phylogeny based on two mitochondrial (COI and NADH) and one nuclear (SSU) gene, and reconstruct morphological evolution among clades. In addition, we prove the existence of mitochondrial mRNA editing for the COI gene. The trees show a lack of conservatism of shell outlines within the main clades, as well as a widespread occurrence of morphological convergences between far-related taxa. Our results refute, therefore, a widespread morphological stasis, which may be an artefact resulting from low taxon coverage. As a corollary, we also revise the groups systematics, notably by emending the large and highly polyphyletic genus Difflugia. These results lead, amongst others, to the erection of a new infraorder Cylindrothecina, as well as two new genera Cylindrifflugia and Golemanskia.

Anthropogenic Stressors in Upland Rivers: Aquatic Macrophyte Responses. A Case Study from Bulgaria

Upland rivers across Europe still exhibit undisturbed conditions and represent a treasure that we cannot afford to lose. We hypothesize that the combination of pristine and modified conditions could demonstrate biological responses along the stressor gradients. Thus, the response of aquatic macrophyte communities to anthropogenic stressors along upland rivers in Bulgaria was studied. Six stressors were selected out of 36 parameters grouped into hydromorphological, chemical variables and combined drivers (catchment land use). The stressors strongly affected species richness on the basis of biological type (bryophytes vs. vascular plants) and ecomorphological type (hydrophytes vs. helophytes). Hydrological alteration expressed by the change of the river's base flow and altered riparian habitats has led to a suppression of bryophytes and a dominance of riverbank plant communities. Seventy-five percent of mountain sites were lacking bryophytes, and the vegetation at semi-mountainous sites was dominated by vascular plants. It can be concluded that hydropeaking, organic and inorganic pollution, and discontinuous urban structures caused important modifications in the aquatic macrophyte assemblages. Macrophyte abundance and the biological and ecomorphological type of aquatic macrophytes reflect multi-stressor effects in upland rivers.

Invasion of Corbicula fluminea (Müller, 1774) (Bivalvia: Corbiculidae) in water bodies from the East Aegean River Basin in Bulgaria

Two decades after the first record of the invasive mussel Corbicula fluminea (Müller, 1774) in Bulgaria in 2001, it has been found in many surface lotic and standing waters in the Bulgarian Danube River Basin. The aim of the current study is to report the finding of C. fluminea in lotic environment of the East Aegean water basin in Bulgaria. Six new localities of the mussel, all situated in the lower stretches of the Tundzha (two) and Maritsa (four) River were registered. The distributional patterns of the mussel and the abundance of its populations in relation to the environmental conditions were analyzed. The adaptive character of the Corbicula, the possible path for invasion as well as the probability of rapid spread of the mussel in the largest rivers and their tributary systems in south-eastern part of Balkan Peninsula was discussed.

Multiple convergences in the evolutionary history of the testate amoeba family Arcellidae (Amoebozoa: Arcellinida: Sphaerothecina): when the ecology rules the morphology

Protists are probably the most species-rich eukaryotes, yet their systematics are inaccurate, leading to an underestimation of their actual diversity. Arcellinida (= lobose testate amoebae) are amoebozoans that build a test (a hard shell) whose shape and composition are taxonomically informative. One of the most successful groups is Arcellidae, a family found worldwide in many freshwater and terrestrial environments where they are indicators of environmental quality. However, the systematics of the family is based on works published nearly a century ago. We re-evaluated the systematics based on single-cell barcoding, morphological and ecological data. Overall, test shape appears to be more related to environmental characteristics than to the species’ phylogenetic position. We show several convergences in organisms with similar ecology, some traditionally described species being paraphyletic. Based on conservative traits, we review the synapomorphies of the infraorder Sphaerothecina, compile a list of synonyms and describe a new genus Galeripora, with five new combinations. Seven new species: Arcella guadarramensis sp. nov., Galeripora balari sp. nov., Galeripora bufonipellita sp. nov., Galeripora galeriformis sp. nov., Galeripora naiadis sp. nov., Galeripora sitiens sp. nov. andGaleripora succelli sp. nov. are also described here.

Checklist of Sphagnum-dwelling testate amoebae in Bulgaria

BACKGROUND

Until now, a complete checklist of Sphagnum-dwelling testate amoebae in Bulgaria has never been published. Records for species diversity and distribution in the country were scattered in many faunistic and ecological publications. The aim of the present study is to summarise all data for the species distribution at the level of country by reviewing the existing literature and by additional data obtained in our research over the past two years.

NEW INFORMATION

The checklist comprises 171 species, classified into 43 genera, 20 families, three orders, three classes and three phyla. We present data for 16 new Sphagnum-dwelling testate amoebae in Bulgaria and new distribution data for 134 species. Of them, 99 species are recorded from Stara Planina Mt., for which there was no available data to date. Additionally are recorded 69 new species for Pirin Mt., 21 for Vitosha Mt. and 18 for Rila Mt. Thirty six species are synonymised according to the latest taxonomic changes. Two misidentified taxa (Euglypha brachiata Penard, 1902 and Difflugia compressa var. africana Gauthier-Lièvre et Thomas, 1958) are transferred into valid species E. acanthophora and Zivkovicia compressa, respectively. Three of the recorded species have not been included in the checklist, because they are currently not refering to testate amoebae (Cochliopodium bilimbosum (Auerbach 1856) and Cochliopodium echinatum Korotneef, 1879 are gymnamoebae (naked amoebae) and Microgromia elegantula (Penard 1904) = Paralieberkuehnia elegantula (Penard 1904) is freshwater foraminifera).

Hamiltonian models for the propagation of irrotational surface gravity waves over a variable bottom

A single incompressible, inviscid, irrotational fluid medium bounded by a free surface and varying bottom is considered. The Hamiltonian of the system is expressed in terms of the so-called Dirichlet-Neumann operators. The equations for the surface waves are presented in Hamiltonian form. Specific scaling of the variables is selected which leads to approximations of Boussinesq and Korteweg-de Vries (KdV) types, taking into account the effect of the slowly varying bottom. The arising KdV equation with variable coefficients is studied numerically when the initial condition is in the form of the one-soliton solution for the initial depth.This article is part of the theme issue 'Nonlinear water waves'.

A contribution to the phylogeny of agglutinating Arcellinida (Amoebozoa) based on SSU rRNA gene sequences

Arcellinid testate amoebae include a wide variety of amoeboid organisms whose test (shell) varies in shape, composition and size. A decade ago, we initiated molecular phylogenetic analyses based on SSU rRNA gene sequences and a taxonomic revision of Arcellinida. However, many lineages within Arcellinida still lack molecular data, and the phylogeny of this group is largely incomplete. In this study, we obtained SSU rRNA gene sequences from seven taxa, of which six have agglutinated shell (Difflugia oblonga, D. labiosa, D. gramen, Mediolus corona, Netzelia wailesi, and N. tuberculata), and one has an entirely proteinaceous shell (Arcella intermedia). All species but Difflugia oblonga branched within the recently erected suborder Sphaerothecina, confirming the synapomorphic value of an oviform or discoid shell. Thus, we propose that species with an oviform or discoid shell currently classified within genus Difflugia must be transferred to other genera, thus continuing the process of taxonomic revision of genus Difflugia, the largest Arcellinida genus. We therefore transferred the current and the previously sequenced oviform Difflugia spp. to Netzelia spp., based on the shared globular/oviform shell shape and their monophyly. Another species, D. labiosa, formed an independent lineage that branched as a sister clade to Arcella spp.; based on the shell morphology and their phylogenetic position, we considered D. labiosa as incertae sedis.

On the morphology, biometry and biogeography of Lamtopyxiscallistoma (Amoebozoa: Arcellinida)

The ultra-structure of the shell and the morphometric variability of soil inhabiting testate amoeba Lamtopyxiscallistoma from Madagascar were studied by using light- and scanning electron microscopy. The biometrical characteristic of the species was made on the basis of 75 specimens measured. In addition to the diameter of the shell, six other shell characters were described biometrically for the first time. The analysis of the variation coefficients shows that the studied population of L.callistoma is comparatively homogeneous and almost all measured characters are weakly to moderate variable (CV less than 10%). Scanning electron microscopy (SEM) studies on the shell ultra-morphology show that it has a smooth apertural surface with a thick layer of porous and fibrous organic cement and a rough dorsal surface composed of bigger and angular pieces of quartz. The shell wall has a thickness of about 5-6 µm and is composed of three layers. Unlike the previously accepted opinion that species is characterized by the presence of four teeth, this study shows that population of L.callistoma from Madagascar is comprised of both, specimens with four teeth and specimens with three teeth, in ratio of about 60% to 40%. Taking into account the restricted geographical distribution, large sizes and characteristic apertural morphology of L.callistoma it is assumed that this species, like some bryophilic ‘Nebelas’ with circumaustral distribution (e.g. Apoderavas, Alocoderacockayni, Certesellacertesi, Certesellamartiali, etc.), can be used as an example that in free-living microbial eukaryotes ‘not everything is everywhere’.

Morphological and molecular diversification of Asian endemic Difflugia tuberspinifera (Amoebozoa, Arcellinida): a case of fast morphological evolution in protists?

Planktonic arcellinid testate amoebae exhibit a broad-range of morphological variability but it is currently unclear to what extent this variability represents phenotypic plasticity or if it is genetically determined. We investigated the morphology and phylogenetic relationships of three endemic east-asian Difflugia taxa 1) the vase-shaped D. mulanensis, 2) and a spinose and a spineless morphotypes of D. tuberspinifera using scanning electron microscopy and two ribosomal genetic markers (SSU rDNA and ITS sequences). Our phylogenetic analyses shows that all three taxa are genetically distinct and closely related to D. achlora and Netzelia oviformis. The genetic variations between the spineless and spinose morphotypes of D. tuberspinifera were low at the SSU rRNA level (0.4%), but ten times higher at the ITS level (4.5-6%). Our data suggest that the two forms of D. tuberspinifera are sufficiently differentiated in terms of morphology and genetic characteristics to constitute two separate entities and that the presence of spines does not result from phenotypic plasticity due to environmental selective pressure. However further observational and experimental data are needed to determine if these two forms constitute different biological species.

Simulation of chemical metabolism for fate and hazard assessment. V. Mammalian hazard assessment

Animals and humans are exposed to a wide array of xenobiotics and have developed complex enzymatic mechanisms to detoxify these chemicals. Detoxification pathways involve a number of biotransformations, such as oxidation, reduction, hydrolysis and conjugation reactions. The intermediate substances created during the detoxification process can be extremely toxic compared with the original toxins, hence metabolism should be accounted for when hazard effects of chemicals are assessed. Alternatively, metabolic transformations could detoxify chemicals that are toxic as parents. The aim of the present paper is to describe specificity of eukaryotic metabolism and its simulation and incorporation in models for predicting skin sensitization, mutagenicity, chromosomal aberration, micronuclei formation and estrogen receptor binding affinity implemented in the TIMES software platform. The current progress in model refinement, data used to parameterize models, logic of simulating metabolism, applicability domain and interpretation of predictions are discussed. Examples illustrating the model predictions are also provided.

Androgen receptor binding affinity: a QSAR evaluation

The multiparameter formulation of the COmmon REactivity PAttern (COREPA) approach has been used to describe the structural requirements for eliciting rat androgen receptor (AR) binding affinity, accounting for molecular flexibility. Chemical affinity for AR binding was related to the distances between nucleophilic sites and structural features describing electronic and hydrophobic interactions between the receptor and ligands. Categorical models were derived for each binding affinity range in terms of specific distances, local (maximal donor delocalizability associated with the oxygen atom of the A ring), global nucleophilicity (partial positive surface areas and energy of the highest occupied molecular orbital) and hydrophobicity (log Kow) of the molecules. An integral screening tool for predicting binding affinity to AR was constructed as a battery of models, each associated with different activity bins. The quality of the screening battery of models was assessed using a high value (0.9) of the Pearson contingency coefficient. The predictability of the model was assessed by testing the model performance on external validation sets. A recently developed technique for selection of potential androgenically active chemicals was used to test the performance of the model in its applicability domain. Some of the selected chemicals were tested for AR transcriptional activation. The experimental results confirmed the theoretical predictions.

Can mutagenicity information be useful in an Integrated Testing Strategy (ITS) for skin sensitization?

Our previous work has investigated the utility of mutagenicity data in the development and application of Integrated Testing Strategies (ITS) for skin sensitization by focusing on the chemical mechanisms at play and substantiating these with experimental data where available. The hybrid expert system TIMES (Tissue Metabolism Simulator) was applied in the identification of the chemical mechanisms since it encodes a comprehensive set of established structure-activity relationships for both skin sensitization and mutagenicity. Based on the evaluation, the experimental determination of mutagenicity was thought to be potentially helpful in the evaluation of skin sensitization potential. This study has evaluated the dataset reported by Wolfreys and Basketter (Cutan. Ocul. Toxicol. 23 (2004), pp. 197-205). Upon an update of the experimental data, the original reported concordance of 68% was found to increase to 88%. There were several compounds that were 'outliers' in the two experimental evaluations which are discussed from a mechanistic basis. The discrepancies were found to be mainly associated with the differences between skin and liver metabolism. Mutagenicity information can play a significant role in evaluating sensitization potential as part of an ITS though careful attention needs to be made to ensure that any information is interpreted in the appropriate context.

Comparing potential COI and SSU rDNA barcodes for assessing the diversity and phylogenetic relationships of cyphoderiid testate amoebae (Rhizaria: Euglyphida)

The mitochondrial Cytochrome Oxidase Subunit 1 gene (COI) has been promoted as an ideal "DNA barcode" for animal species and other groups of eukaryotes. However, the utility of the COI marker for species level discrimination and for phylogenetic analyses has yet to be tested within the Rhizaria. Accordingly, we analysed mitochondrial COI gene sequences and nuclear small subunit rDNA (SSU) sequences from several morphospecies of euglyphid testate amoebae (Cercozoa, Rhizaria) in order to evaluate the utility of these DNA markers for species discrimination and phylogenetic reconstructions. Sequences were obtained from eleven populations belonging to sixCyphoderiamorphospecies that were isolated from field samples in North America and Europe. Mean inter-population COI sequence dissimilarities were on average 2.9 times greater than in the SSU, while the intra-population sequence dissimilarities were higher in the SSU (0-0.95%) than in the COI (0%); this suggests that the COI fragment is valuable for discriminating Cyphoderiidae isolates. Our study also demonstrated that COI sequences are useful for inferring phylogenetic relationships among Cyphoderiidae isolates. COI and SSU tree topologies were very similar even though the COI fragment used in these analyses (500bp) was much shorter than the SSU sequences (1600bp). Altogether, these results demonstrate the utility of the COI as a potential taxonomic DNA barcode for assessing cyphoderiid species diversity and for inferring phylogenetic relationships within the group.

Strategic selection of chemicals for testing. Part I. Functionalities and performance of basic selection methods

To develop quantitative structure-activity relationships (QSAR) models capable of predicting adverse effects for large chemical inventories and diverse structures, an interactive approach is presented that includes testing of strategically selected chemicals to expand the scope of a preliminary model to cover a target inventory. The goal of chemical selection in this context is to make the testing more effective in terms of adding maximal new structural information to the predictive model with minimal testing. The aim of this paper is to describe a set of algorithmic solutions and modelling techniques that can be used to efficiently select chemicals for testing to achieve a variety of goals. One purpose of chemical selection is to refine the model thus extending our knowledge about the modelled subject. Alternatively, the selection of chemicals for testing could be targeted at achieving a more adequate structural representation of a specific universe of untested chemicals to extend the model applicability domain on each subsequent step of model development. The chemical selection tools are collectively provided in a software package referred to as ChemPick. The system also allows the visualisation of chemical inventories and training sets in multidimensional (two and three dimensions) descriptor space. The software environment allows one or more datasets to be simultaneously loaded in a three-dimensional (or N-dimensional) chart where each point represents a combination of values for the descriptors for a given conformation of a chemical. The application of the chemical selection tools to select chemicals to expand a preliminary model of human oestrogen receptor (hER) ligand binding to more adequately cover a diverse chemical inventory is presented to demonstrate the application of these tools.

QSAR and mechanistic interpretation of estrogen receptor binding

A multi-dimensional formulation of the COmmon REactivity PAttern (COREPA) modeling approach has been used to investigate chemical binding to the human estrogen receptor (hER). A training set of 645 chemicals included 497 steroid and environmental chemicals (database of the Chemical Evaluation and Research Institute, Japan - CERI) and 148 chemicals to further explore hER-structure interactions (selected J. Katzenellenbogen references). Upgrades of modeling approaches were introduced for multivariate COREPA analysis, optimal conformational generation and description of the local hydrophobicity of chemicals. Analysis of reactivity patterns based on the distance between nucleophilic sites resulted in identification of distinct interaction types: a steroid-like A-B type described by frontier orbital energies and distance between nucleophilic sites with specific charge requirements; an A-C type where local hydrophobic effects are combined with electronic interactions to modulate binding; and mixed A-B-C (AD) type. Chemicals were grouped by type, then COREPA models were developed for within specific relative binding affinity ranges of >10%, 10 > RBA > or = 0.1%, and 0.1 > RBA > 0.0%. The derived models for each interaction type and affinity range combined specific prefiltering requirements (interatomic distances) and a COREPA classification node using no more than 2 discriminating parameters. The interaction types are becoming less distinct in the lowest activity range for each chemicals of each type; here, the modeling was performed within chemical classes (phenols, phthalates, etc.). The ultimate model was organized as a battery of local models associated to interaction type and mechanism.

Identification of the structural requirements for mutagencitiy, by incorporating molecular flexibility and metabolic activation of chemicals. II. General Ames mutagenicity model

The tissue metabolic simulator (TIMES) modeling approach is a hybrid expert system that couples a metabolic simulator together with structure toxicity rules, underpinned by structural alerts, to predict interaction of chemicals or their metabolites with target macromolecules. Some of the structural alerts representing the reactivity pattern-causing effect could interact directly with the target whereas others necessitated a combination with two- or three-dimensional quantitative structure-activity relationship models describing the firing of the alerts from the rest of the molecules. Recently, TIMES has been used to model bacterial mutagenicity [Mekenyan, O., Dimitrov, S., Serafimova, R., Thompson, E., Kotov, S., Dimitrova, N., and Walker, J. (2004) Identification of the structural requirements for mutagenicity by incorporating molecular flexibility and metabolic activation of chemicals I: TA100 model. Chem. Res. Toxicol. 17 (6), 753-766]. The original model was derived for a single tester strain, Salmonella typhimurium (TA100), using the Ames test by the National Toxicology Program (NTP). The model correctly identified 82% of the primary acting mutagens, 94% of the nonmutagens, and 77% of the metabolically activated chemicals in a training set. The identified high correlation between activities across different strains changed the initial strategic direction to look at the other strains in the next modeling developments. In this respect, the focus of the present work was to build a general mutagenicity model predicting mutagenicity with respect to any of the Ames tester strains. The use of all reactivity alerts in the model was justified by their interaction mechanisms with DNA, found in the literature. The alerts identified for the current model were analyzed by comparison with other established alerts derived from human experts. In the new model, the original NTP training set with 1341 structures was expanded by 1626 proprietary chemicals provided by BASF AG. Eventually, the training set consisted of 435 chemicals, which are mutagenic as parents, 397 chemicals that are mutagenic after S9 metabolic activation, and 2012 nonmutagenic chemicals. The general mutagenicity model was found to have 82% sensitivity, 89% specificity, and 88% concordance for training set chemicals. The model applicability domain was introduced accounting for similarity (structural, mechanistic, etc.) between predicted chemicals and training set chemicals for which the model performs correctly.