Molecular and life-history effects of a natural toxin on herbivorous and non-target soil arthropods

Assessment of the reproduction of six collembolan species in tropical soils naturally rich in potentially toxic elements

Laboratory ecotoxicological tests are important tools for the management of environmental changes derived from anthropogenic activities. Folsomia candida is usually the model species used in some procedures. However, this species may not be sufficiently representative of the sensitivity of the other collembolan species. This study aimed to evaluate (i) the effects of soils naturally rich in potentially toxic elements (PTE) and soil characteristics on the reproduction and survival of different collembolan species, (ii) whether the habitat function of these soils is compromised, and (iii) to what extent F. candida is representative of the other collembolan species. For this, reproduction tests with six collembolan species were conducted in 14 different samples of soils. In general, collembolan reproduction was not completely inhibited in none of the natural tested soils. Even soils with high pollution load index values did not negatively affect collembolan reproduction for most of the species. In contrast, the lowest collembolan reproduction rates were found in a visually dense soil (lowest volume/weight ratio), highlighting that soil attributes other than total PTE concentration also interfere in the reproduction of collembolans. Our results support the idea that the F. candida species might not be representative of other collembolan species and that laboratory tests to assess soil contaminations should be conducted using diverse collembolan species.

Study on molecular level toxicity of Sb(V) to soil springtails: using a combination of transcriptomics and metabolomics

To date, numerous studies have focused on the toxicity of antimony (Sb) to soil-dwelling organisms at the individual level. However, little is known about Sb-caused molecular level toxicity. Here, an integrated transcriptomics and metabolomics approach was used to better reveal toxicity of Sb(V) to springtails Folsomia candida considering environmentally relevant speciation of Sb. No significant effects of Sb(V) on survival, reproduction and growth of springtails were observed using the ISO standard test. Transcriptomics analysis identified 1015 and 3367 differentially expressed genes (DEGs) after 2 and 7 d of exposure, indicating an increasing transcriptomal changes with time. Significantly enriched top GO (Gene Ontology) terms (chitin metabolic process, chitin binding and extracellular region) were shared between the two time exposure groups. However, no enriched KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway was shared, with fatty acid metabolism and apoptosis-fly being the most significant pathway, respectively. Metabolomics analysis identified 155 differential changed metabolites (DCMs) in springtails after 7 d of exposure. Antifolate resistance was the most significantly enriched pathway, in which dihydrofolic acid was up-regulated and three purine nucleotides (adenosine 5'-monophosphate, inosine 5'-monophosphate, guanosine 5'-monophosphate) were down-regulated. This indicated obvious repression of DNA replication, which was also observed by transcriptomics. Additionally, metabolites level related to chitin, oxidative stress, and protein metabolism significantly changed, and these metabolites could also support and confirm main transcriptomic results. Thus, the combination of multiomics facilitated better understanding of the molecular level of toxicity of Sb(V) in Collembola.

Bringing ecology into toxicology: Life-cycle toxicity of two neonicotinoids to four different species of springtails in LUFA 2.2 natural soil

Collembolans comprise one of the most abundant groups of soil invertebrates within the arthropods. The parthenogenetic species Folsomia candida (Willem, 1902) is the most well-studied representative, being used since the beginning of the 1960s as a model organism for assessing toxicity of chemicals in soil. In this paper we aimed at answering three questions by exposing four different species of springtails (F. candida, Folsomia fimetaria, Sinella curviseta and Heteromurus nitidus) to the neonicotinoids imidacloprid and thiacloprid: i) How representative as a model organism is F. candida for species of springtails that reproduce sexually? (ii) How suitable are other species of springtails to be used as model organisms for ecotoxicological testing? (iii) Is it possible to use the life history of these species to extrapolate the impact of neonicotinoids on the population level? Our results showed that F. candida is a good model organism, despite being the most sensitive species tested, when analysing both endpoints - survival and reproduction. The tests performed with S. curviseta and H. nitidus showed that they could be used as surrogates in ecotoxicity tests, and also to predict how their population might be affected after being exposed to chemicals. The adjustments made to the test performed with F. candida: introducing adults (20-22 days old) into the test jars and exposing them for 21 days instead of 28 days, proved to be as efficient as the standardized test guideline (OECD 232, 2009).

Agricultural trichothecene mycotoxin contamination affects the life-history and reduced glutathione content of Folsomia candida Willem (Collembola)

There is limited data available concerning the effect of T-2/HT-2 toxin or deoxynivalenol (DON) on invertebrates such as springtails, and no data on their life history and oxidative stress. Control maize and DON or T-2 toxin contaminated maize were fed to Folsomia candida with a toxin content of 16324 mg DON kg–1 or 671 mg T-2 kg–1 maize. Ten to twelve days old animals were investigated in a life-history test and a stress protein test.T-2 toxin did not affect Folsomia candida in any measured parameters. The DON exposed group showed decreased growth and reproduction, and a higher survival rate. DON treatment resulted in lower protein content, while reduced glutathione content was higher than in control. It suggests that DON activated the glutathione-related detoxification pathway, which possibly causes a higher survival rate. The results also suggest that the oral toxicity of DON or T-2 is lower than through physical contact.For that reason, DON and T-2 toxin contaminated maize is not suggested to be used as green manure in the native state. Alternative solutions could be using mycotoxin contaminated maize for biogas production, or after decontamination by bacterial strains, it can be used as organic fertilizer.

Biofumigation for Fighting Replant Disease- A Review

Replant disease is a soil (micro-) biome-based, harmfully-disturbed physiological and morphological reaction of plants to replanting similar cultures on the same sites by demonstrating growth retardation and leading to economic losses especially in Rosaceae plant production. Commonly, replant disease is overcome by soil fumigation with toxic chemicals. With chemical soil fumigation being restricted in many countries, other strategies are needed. Biofumigation, which is characterized by the incorporation of Brassicaceae plant materials into soil, is a promising method. We review the potential of biofumigation in the fight against replant disease. Biofumigation using optimized Brassicaceae seed meal compositions in combination with replant disease tolerant plant genotypes shows promising results, but the efficacy is still soil and site-dependent. Therefore, future studies should address the optimal timing as well as amount and type of incorporated plant material and environmental conditions during incubation in dependence of the soil physical and chemical characteristics.

Transcriptome-level effects of the model organic pollutant phenanthrene and its solvent acetone in three amphipod species

Polyaromatic hydrocarbons (PAH) are common pollutants of water ecosystems originating from incineration processes and contamination with mineral oil. Water solubility of PAHs is generally low; for toxicity tests with aquatic organisms, they are therefore usually dissolved in organic solvents. Here we examined the effects of a typical model PAH, phenanthrene, and a solvent, acetone, on amphipods as relevant aquatic invertebrate models. Two of these species, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, are common endemics of the oligotrophic and pristine Lake Baikal, while one, Gammarus lacustris, is widespread throughout the Holarctic and inhabits smaller and more eutrophic water bodies in the Baikal area. Neither solvent nor phenanthrene caused mortality at the applied concentrations, but both substances affected gene expression in all species. Differential gene expression was more profound in the species from Lake Baikal than in the Holarctic species. Moreover, in one of the Baikal species, E. cyaneus, we found that many known components of the cellular xenobiotic detoxification system reacted to the treatments. Finally, we detected a negative relationship between changes in transcript abundances in response to the solvent and phenanthrene. This mixture effect, weaker than the impact by a single mixture component, needs further exploration.

Population Genomics Insights into Adaptive Evolution and Ecological Differentiation in Streptomycetes

Deciphering the genomic variation that represents microevolutionary processes toward species divergence is key to understanding microbial speciation, which has long been under debate. Streptomycetes are filamentous bacteria that are ubiquitous in nature and the richest source of antibiotics; however, their speciation processes remain unknown. To tackle this issue, we performed a comprehensive population genomics analysis on Streptomyces albidoflavus residing in different habitats and with a worldwide distribution and identified and characterized the foundational changes within the species. We detected three well-defined phylogenomic clades, of which clades I and III mainly contained free-living (soil/marine) and insect-associated strains, respectively, and clade II had a mixed origin. By performing genome-wide association studies (GWAS), we identified a number of genetic variants associated with free-living or entomic (denoting or relating to insects) habitats in both the accessory and core genomes. These variants contributed collectively to the population structure and had annotated or confirmed functions that likely facilitate differential adaptation of the species. In addition, we detected higher levels of homologous recombination within each clade and in the free-living group than within the whole species and in the entomic group. A subset of the insect-associated strains (clade III) showed a relatively independent evolutionary trajectory with more symbiosis-favorable genes but little genetic interchange with the other lineages. Our results demonstrate that ecological adaptation promotes genetic differentiation in S. albidoflavus, suggesting a model of ecological speciation with gene flow in streptomycetes.IMPORTANCE Species are the fundamental units of ecology and evolution, and speciation leads to the astounding diversity of life on Earth. Studying speciation is thus of great significance to understand, protect, and exploit biodiversity, but it is a challenge in the microbial world. In this study, using population genomics, we placed Streptomyces albidoflavus strains in a spectrum of speciation and showed that the genetic differences between phylogenomic clusters evolved mainly by environmental selection and gene-specific sweeps. These findings highlight the role of ecology in structuring recombining bacterial species, making a step toward a deeper understanding of microbial speciation. Our results also raise concerns of an underrated microbial diversity at the intraspecies level, which can be utilized for mining of ecologically relevant natural products.

Optimizing isothiocyanate formation during enzymatic glucosinolate breakdown by adjusting pH value, temperature and dilution in Brassica vegetables and Arabidopsis thaliana

Consumption of glucosinolate-rich Brassicales vegetables is associated with a decreased risk of cancer with enzymatic hydrolysis of glucosinolates playing a key role. However, formation of health-promoting isothiocyanates is inhibited by the epithiospecifier protein in favour of nitriles and epithionitriles. Domestic processing conditions, such as changes in pH value, temperature or dilution, might also affect isothiocyanate formation. Therefore, the influences of these three factors were evaluated in accessions of Brassica rapa, Brassica oleracea, and Arabidopsis thaliana. Mathematical modelling was performed to determine optimal isothiocyanate formation conditions and to obtain knowledge on the kinetics of the reactions. At 22 °C and endogenous plant pH, nearly all investigated plants formed nitriles and epithionitriles instead of health-promoting isothiocyanates. Response surface models, however, clearly demonstrated that upon change in pH to domestic acidic (pH 4) or basic pH values (pH 8), isothiocyanate formation considerably increases. While temperature also affects this process, the pH value has the greatest impact. Further, a kinetic model showed that isothiocyanate formation strongly increases due to dilution. Finally, the results show that isothiocyanate intake can be strongly increased by optimizing the conditions of preparation of Brassicales vegetables.

Combined Transcriptomics Analysis for Classification of Adverse Effects As a Potential End Point in Effect Based Screening

Environmental risk assessment relies on the use of bioassays to assess the environmental impact of chemicals. Gene expression is gaining acceptance as a valuable mechanistic end point in bioassays and effect-based screening. Data analysis and its results, however, are complex and often not directly applicable in risk assessment. Classifier analysis is a promising method to turn complex gene expression analysis results into answers suitable for risk assessment. We have assembled a large gene expression data set assembled from multiple studies and experiments in the springtail Folsomia candida, with the aim of selecting a set of genes that can be trained to classify general toxic stress. By performing differential expression analysis prior to classifier training, we were able to select a set of 135 genes which was enriched in stress related processes. Classifier models from this set were used to classify two test sets comprised of chemical spiked, polluted, and clean soils and compared to another, more traditional classifier feature selection. The gene set presented here outperformed the more traditionally selected gene set. This gene set has the potential to be used as a biomarker to test for adverse effects caused by chemicals in springtails to provide end points in environmental risk assessment.

Verticillium longisporum infection induces organ-specific glucosinolate degradation in Arabidopsis thaliana

The species Verticillium represents a group of highly destructive fungal pathogens, responsible for vascular wilt in a number of crops. The host response to infection by Verticillium longisporum at the level of secondary plant metabolites has not been well explored. Natural variation in the glucosinolate (GLS) composition of four Arabidopsis thaliana accessions was characterized: the accessions Bur-0 and Hi-0 accumulated alkenyl GLS, while 3-hydroxypropyl GLS predominated in Kn-0 and Ler-0. With respect to GLS degradation products, Hi-0 and Kn-0 generated mainly isothiocyanates, whereas Bur-0 released epithionitriles and Ler-0 nitriles. An analysis of the effect on the composition of both GLS and its breakdown products in the leaf and root following the plants' exposure to V. longisporum revealed a number of organ- and accession-specific alterations. In the less disease susceptible accessions Bur-0 and Ler-0, colonization depressed the accumulation of GLS in the rosette leaves but accentuated it in the roots. In contrast, in the root, the level of GLS breakdown products in three of the four accessions fell, suggestive of their conjugation or binding to a fungal target molecule(s). The plant-pathogen interaction influenced both the organ- and accession-specific formation of GLS degradation products.

Effects of a natural toxin on life history and gene expression of Eisenia andrei

Earthworms perform key functions for a healthy soil ecosystem, such as bioturbation. The soil ecosystem can be challenged by natural toxins such as isothiocyanates (ITCs), produced by many commercial crops. Therefore, the effects of 2-phenylethyl ITC were investigated on the earthworm Eisenia andrei using an ecotoxicogenomics approach. Exposure to 2-phenylethyl ITC reduced both survival and reproduction of E. andrei in a dose-dependent manner (median effective concentration [EC50] = 556 nmol/g). Cross-species comparative genomic hybridization validated the applicability of an existing 4 × 44,000 Eisenia fetida microarray to E. andrei. Gene expression profiles revealed the importance of metallothionein (MT) as an early warning signal when E. andrei was exposed to low concentrations of 2-phenylethyl ITC. Alignment of these MT genes with the MT-2 gene of Lumbricus rubellus showed that at least 2 MT gene clusters are present in the Eisenia sp. genome. At high-exposure concentrations, gene expression was mainly affected by inhibiting chitinase activity, inducing an oxidative stress response, and stimulating energy metabolism. Furthermore, analysis by Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway implied that the high concentration may have caused impaired light sensitivity, angiogenesis, olfactory perception, learning, and memory. Increased levels of ITCs may be found in the field in the near future. The results presented call for a careful investigation to quantify the risk of such compounds before allowing them to enter the soil on a large scale.

Genetic variation in parthenogenetic collembolans is associated with differences in fitness and cadmium-induced transcriptome responses

Ecotoxicological tests may be biased by the use of laboratory strains that usually contain very limited genetic diversity. It is therefore essential to study how genetic variation influences stress tolerance relevant for toxicity outcomes. To that end we studied sensitivity to cadmium in two distinct genotypes of the parthogenetic soil ecotoxicological model organism Folsomia candida. Clonal lines of both genotypes (TO1 and TO2) showed divergent fitness responses to cadmium exposure; TO2 reproduction was 20% less affected by cadmium. Statistical analyses revealed significant differences between the cadmium-affected transcriptomes: i) the number of genes affected by cadmium in TO2 was only minor (~22%) compared to TO1; ii) 97 genes showed a genotype × cadmium interaction and their response to cadmium showed globally larger fold changes in TO1 when compared to TO2; iii) the interaction genes showed a concerted manner of expression in TO1, while a less coordinated pattern was observed in TO2. We conclude that (1) there is genetic variation in parthenogenetic populations of F. candida, and (2) this variation affects life-history and molecular end points relative to cadmium toxicity. This sheds new light on the sources of biological variability in test results, even when the test organisms are thought to be genetically homogeneous because of their parthenogenetic reproduction.

Time-related survival effects of two gluconasturtiin hydrolysis products on the terrestrial isopod Porcellio scaber

Glucosinolates are compounds produced by commercial crops which can hydrolyse in a range of natural toxins that may exert detrimental effects on beneficial soil organisms. This study examined the effects of 2-phenylethyl isothiocyanate and 3-phenylpropionitrile on the survival and growth of the woodlouse Porcellio scaber exposed for 28 d. 2-Phenylethyl isothiocyanate dissipated from the soil with half-lives ranging from 19 to 96 h. Exposure through soil showed toxic effects only on survival. The LC50s after 28 d were significantly different at 65.3 mg kg(-1) for 2-phenylethyl isothiocyanate and 155 mg kg(-1) for 3-phenylpropionitrile. A toxicokinetic-toxicodynamic (TKTD) approach, however, revealed that both compounds in fact have very similar effect patterns. The TKTD model was better suited to interpret the survival data than descriptive dose-response analysis (LC(x)), accounting for the fast dissipation of the compounds in the soil. Found effects were within environmentally relevant concentrations. Care should therefore be taken before allowing these natural toxins to enter soil ecosystems in large quantities.

Genes involved in the evolution of herbivory by a leaf-mining, Drosophilid fly

Herbivorous insects are among the most successful radiations of life. However, we know little about the processes underpinning the evolution of herbivory. We examined the evolution of herbivory in the fly, Scaptomyza flava, whose larvae are leaf miners on species of Brassicaceae, including the widely studied reference plant, Arabidopsis thaliana (Arabidopsis). Scaptomyza flava is phylogenetically nested within the paraphyletic genus Drosophila, and the whole genome sequences available for 12 species of Drosophila facilitated phylogenetic analysis and assembly of a transcriptome for S. flava. A time-calibrated phylogeny indicated that leaf mining in Scaptomyza evolved between 6 and 16 million years ago. Feeding assays showed that biosynthesis of glucosinolates, the major class of antiherbivore chemical defense compounds in mustard leaves, was upregulated by S. flava larval feeding. The presence of glucosinolates in wild-type (WT) Arabidopsis plants reduced S. flava larval weight gain and increased egg-adult development time relative to flies reared in glucosinolate knockout (GKO) plants. An analysis of gene expression differences in 5-day-old larvae reared on WT versus GKO plants showed a total of 341 transcripts that were differentially regulated by glucosinolate uptake in larval S. flava. Of these, approximately a third corresponded to homologs of Drosophila melanogaster genes associated with starvation, dietary toxin-, heat-, oxidation-, and aging-related stress. The upregulated transcripts exhibited elevated rates of protein evolution compared with unregulated transcripts. The remaining differentially regulated transcripts also contained a higher proportion of novel genes than the unregulated transcripts. Thus, the transition to herbivory in Scaptomyza appears to be coupled with the evolution of novel genes and the co-option of conserved stress-related genes.