Increase in egg resistance to desiccation in springtails correlates with blastodermal cuticle formation: Eco-evolutionary implications for insect terrestrialization.

Land colonization was a major event in the history of life. Among animals, insects exerted a staggering terrestrialization success, due to traits usually associated with postembryonic life stages, while the egg stage has been largely overlooked in comparative studies. In many insects, after blastoderm differentiation, the extraembryonic serosal tissue wraps the embryo and synthesizes the serosal cuticle, an extracellular matrix that lies beneath the eggshell and protects the egg against water loss. In contrast, in noninsect hexapods such as springtails (Collembola) the early blastodermal cells synthesize a blastodermal cuticle. Here, we investigate the relationship between blastodermal cuticle formation and egg resistance to desiccation in the springtails Orchesella cincta and Folsomia candida, two species with different oviposition environments and developmental rates. The blastodermal cuticle becomes externally visible in O. cincta and F. candida at 22% and 29% of embryogenesis, respectively. To contextualize, we describe the stages of springtail embryogenesis, exemplified by F. candida. Our physiological assays then showed that blastodermal cuticle formation coincides with an increase in egg viability in a dry environment, significantly contributing to hatching success. However, protection differs between species: while O. cincta eggs survive at least 2 hr outside a humid environment, the survival period recorded for F. candida eggs is only 15 min, which correlates with this species' requirement for humid microhabitats. We suggest that the formation of this cuticle protects the eggs, constituting an ancestral trait among hexapods that predated and facilitated the process of terrestrialization that occurred during insect evolution.

Multi-faceted analysis provides little evidence for recurrent whole-genome duplications during hexapod evolution.

BACKGROUND: Gene duplication events play an important role in the evolution and adaptation of organisms. Duplicated genes can arise through different mechanisms, including whole-genome duplications (WGDs). Recently, WGD was suggested to be an important driver of evolution, also in hexapod animals. RESULTS: Here, we analyzed 20 high-quality hexapod genomes using whole-paranome distributions of estimated synonymous distances (KS), patterns of within-genome co-linearity, and phylogenomic gene tree-species tree reconciliation methods. We observe an abundance of gene duplicates in the majority of these hexapod genomes, yet we find little evidence for WGD. The majority of gene duplicates seem to have originated through small-scale gene duplication processes. We did detect segmental duplications in six genomes, but these lacked the within-genome co-linearity signature typically associated with WGD, and the age of these duplications did not coincide with particular peaks in KS distributions. Furthermore, statistical gene tree-species tree reconciliation failed to support all but one of the previously hypothesized WGDs. CONCLUSIONS: Our analyses therefore provide very limited evidence for WGD having played a significant role in the evolution of hexapods and suggest that alternative mechanisms drive gene duplication events in this group of animals. For instance, we propose that, along with small-scale gene duplication events, episodes of increased transposable element activity could have been an important source for gene duplicates in hexapods.

Transgenerational and multigenerational stress gene responses to the insecticide etofenprox in Folsomia candida (Collembola).

Insecticide exposure may cause both transgenerational and multigenerational effects on populations, but the molecular mechanisms of these changes remain largely unclear. Many studies have focused on either transgenerational or multigenerational mechanisms but did neglect the comparative aspects. This study assessed whether the pyrethroid insecticide etofenprox (formulation Trebon® 30 EC) shows transgenerational and/or multigenerational effects on the survival and reproduction of Folsomia candida (Collembola). The activation of stress-related genes was studied to detect whether etofenprox modifies the expression of reproduction-associated genes in trans- and multigenerational treatments. A laboratory study was carried out for three generations with five insecticide concentrations in LUFA 2.2 soil. In the transgenerational treatment, only the parent generation (P) was exposed, but the subsequent generations were not. In the multigenerational treatment, all three generations were exposed to the insecticide in the same manner. Multigenerational exposure resulted in reduced reproduction effects over generations, suggesting that F. candida is capable of acclimating to enhanced concentration levels of etofenprox during prolonged exposure over multiple generations. In the transgenerational treatment, the heat shock protein 70 was up-regulated and cytochrome oxidase 6N4v1 expression down-regulated in a dose-dependent manner in the F2 generation. This finding raises the possibility of the epigenetic inheritance of insecticide impacts on parents. Furthermore, CYP6N4v1 expression was oppositely regulated in the trans- and multigenerational treatments. Our results draw attention to the differences in molecular level responses of F. candida to trans- and multigenerational etofenprox exposure.

Metagenomic analysis of bacterial community structure and diversity of lignocellulolytic bacteria in Vietnamese native goat rumen.

OBJECTIVE: In a previous study, analysis of Illumina sequenced metagenomic DNA data of bacteria in Vietnamese goats' rumen showed a high diversity of putative lignocellulolytic genes. In this study, taxonomy speculation of microbial community and lignocellulolytic bacteria population in the rumen was conducted to elucidate a role of bacterial structure for effective degradation of plant materials. METHODS: The metagenomic data had been subjected into Basic Local Alignment Search Tool (BLASTX) algorithm and the National Center for Biotechnology Information non-redundant sequence database. Here the BLASTX hits were further processed by the Metagenome Analyzer program to statistically analyze the abundance of taxa. RESULTS: Microbial community in the rumen is defined by dominance of Bacteroidetes compared to Firmicutes. The ratio of Firmicutes versus Bacteroidetes was 0.36:1. An abundance of Synergistetes was uniquely identified in the goat microbiome may be formed by host genotype. With regard to bacterial lignocellulose degraders, the ratio of lignocellulolytic genes affiliated with Firmicutes compared to the genes linked to Bacteroidetes was 0.11:1, in which the genes encoding putative hemicellulases, carbohydrate esterases, polysaccharide lyases originated from Bacteroidetes were 14 to 20 times higher than from Firmicutes. Firmicutes seem to possess more cellulose hydrolysis capacity showing a Firmicutes/Bacteroidetes ratio of 0.35:1. Analysis of lignocellulolytic potential degraders shows that four species belonged to Bacteroidetes phylum, while two species belonged to Firmicutes phylum harbouring at least 12 different catalytic domains for all lignocellulose pretreatment, cellulose, as well as hemicellulose saccharification. CONCLUSION: Based on these findings, we speculate that increasing the members of Bacteroidetes to keep a low ratio of Firmicutes versus Bacteroidetes in goat rumen has resulted most likely in an increased lignocellulose digestion.

Coping with living in the soil: the genome of the parthenogenetic springtail Folsomia candida.

BACKGROUND: Folsomia candida is a model in soil biology, belonging to the family of Isotomidae, subclass Collembola. It reproduces parthenogenetically in the presence of Wolbachia, and exhibits remarkable physiological adaptations to stress. To better understand these features and adaptations to life in the soil, we studied its genome in the context of its parthenogenetic lifestyle. RESULTS: We applied Pacific Bioscience sequencing and assembly to generate a reference genome for F. candida of 221.7 Mbp, comprising only 162 scaffolds. The complete genome of its endosymbiont Wolbachia, was also assembled and turned out to be the largest strain identified so far. Substantial gene family expansions and lineage-specific gene clusters were linked to stress response. A large number of genes (809) were acquired by horizontal gene transfer. A substantial fraction of these genes are involved in lignocellulose degradation. Also, the presence of genes involved in antibiotic biosynthesis was confirmed. Intra-genomic rearrangements of collinear gene clusters were observed, of which 11 were organized as palindromes. The Hox gene cluster of F. candida showed major rearrangements compared to arthropod consensus cluster, resulting in a disorganized cluster. CONCLUSIONS: The expansion of stress response gene families suggests that stress defense was important to facilitate colonization of soils. The large number of HGT genes related to lignocellulose degradation could be beneficial to unlock carbohydrate sources in soil, especially those contained in decaying plant and fungal organic matter. Intra- as well as inter-scaffold duplications of gene clusters may be a consequence of its parthenogenetic lifestyle. This high quality genome will be instrumental for evolutionary biologists investigating deep phylogenetic lineages among arthropods and will provide the basis for a more mechanistic understanding in soil ecology and ecotoxicology.

Evolutionary ecology of beta-lactam gene clusters in animals.

Beta-lactam biosynthesis was thought to occur only in fungi and bacteria, but we recently reported the presence of isopenicillin N synthase in a soil-dwelling animal, Folsomia candida. However, it has remained unclear whether this gene is part of a larger beta-lactam biosynthesis pathway and how widespread the occurrence of penicillin biosynthesis is among animals. Here, we analysed the distribution of beta-lactam biosynthesis genes throughout the animal kingdom and identified a beta-lactam gene cluster in the genome of F. candida (Collembola), consisting of isopenicillin N synthase (IPNS), δ-(L-α-aminoadipoyl)-L-cysteinyl-D-valine synthetase (ACVS), and two cephamycin C genes (cmcI and cmcJ) on a genomic scaffold of 0.76 Mb. All genes are transcriptionally active and are inducible by stress (heat shock). A beta-lactam compound was detected in vivo using an ELISA beta-lactam assay. The gene cluster also contains an ABC transporter which is coregulated with IPNS and ACVS after heat shock. Furthermore, we show that different combinations of beta-lactam biosynthesis genes are present in over 60% of springtail families, but they are absent from genome- and transcript libraries of other animals including close relatives of springtails (Protura, Diplura and insects). The presence of beta-lactam genes is strongly correlated with an euedaphic (soil-living) lifestyle. Beta-lactam genes IPNS and ACVS each form a phylogenetic clade in between bacteria and fungi, while cmcI and cmcJ genes cluster within bacteria. This suggests a single horizontal gene transfer event most probably from a bacterial host, followed by differential loss in more recently evolving species.

Evolutionary history of a dispersal-associated locus across sympatric and allopatric divergent populations of a wing-polymorphic beetle across Atlantic Europe.

Studying the evolutionary history of trait divergence, in particular those related to dispersal capacity, is of major interest for the process of local adaptation and metapopulation dynamics. Here, we reconstruct the evolution of different alleles at the nuclear-encoded mitochondrial NADP(+)-dependent isocitrate dehydrogenase (mtIdh) locus of the ground beetle Pogonus chalceus that are differentially and repeatedly selected in short- and long-winged populations in response to different hydrological regimes at both allopatric and sympatric scales along the Atlantic European coasts. We sequenced 2788 bp of the mtIdh locus spanning a ~7-kb genome region and compared its variation with that of two supposedly neutral genes. mtIdh sequences show (i) monophyletic clustering of the short-winged associated mtIDH-DE haplotypes within the long-winged associated mtIDH-AB haplotypes, (ii) a more than tenfold lower haplotype diversity associated with the mtIDH-DE alleles compared to the mtIDH-AB alleles and (iii) a high number of fixed nucleotide differences between both mtIDH haplotype clusters. Coalescent simulations suggest that this observed sequence variation in the mtIdh locus is most consistent with a singular origin in a partially isolated subpopulation, followed by a relatively recent spread of the mtIDH-DE allele in short-winged populations along the Atlantic coast. These results demonstrate that even traits associated with decreased dispersal capacity can rapidly spread and that reuse of adaptive alleles plays an important role in the adaptive potential within this sympatric mosaic of P. chalceus populations.

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.

Transcriptome assembly and microarray construction for Enchytraeus crypticus, a model oligochaete to assess stress response mechanisms derived from soil conditions.

BACKGROUND: The soil worm Enchytraeus crypticus (Oligochaeta) is an ecotoxicology model species that, until now, was without genome or transcriptome sequence information. The present research aims at studying the transcriptome of Enchytraeus crypticus, sampled from multiple test conditions, and the construction of a high-density microarray for functional genomic studies. RESULTS: Over 1.5 million cDNA sequence reads were obtained representing 645 million nucleotides. After assembly, 27,296 contigs and 87,686 singletons were obtained, from which 44% and 25% are annotated as protein-coding genes, respectively, sharing homology with other animal proteomes. Concerning assembly quality, 84% of the contig sequences contain an open reading frame with a start codon while E. crypticus homologs were identified for 92% of the core eukaryotic genes. Moreover, 65% and 77% of the singletons and contigs without known homologs, respectively, were shown to be transcribed in an independent microarray experiment. An Agilent 180 K microarray platform was designed and validated by hybridizing cDNA from 4 day zinc- exposed E. crypticus to the concentration corresponding to 50% reduction in reproduction after three weeks (EC50). Overall, 70% of all probes signaled expression above background levels (mean signal + 1x standard deviation). More specifically, the probes derived from contigs showed a wider range of average intensities when compared to probes derived from singletons. In total, 522 significantly differentially regulated transcripts were identified upon zinc exposure. Several significantly regulated genes exerted predicted functions (e.g. zinc efflux, zinc transport) associated with zinc stress. Unexpectedly, the microarray data suggest that zinc exposure alters retro transposon activity in the E. crypticus genome. CONCLUSION: An initial investigation of the E. crypticus transcriptome including an associated microarray platform for future studies proves to be a valuable resource to investigate functional genomics mechanisms of toxicity in soil environments and to annotate a potentially large number of lineage specific genes that are responsive to environmental stress conditions.

A functional isopenicillin N synthase in an animal genome.

Horizontal transfer of genes is widespread among prokaryotes, but is less common between microorganisms and animals. Here, we present evidence for the presence of a gene encoding functional isopenicillin N synthase, an enzyme in the ß-lactam antibiotics biosynthesis pathway, in the genome of the soil-living collembolan species, Folsomia candida (FcIPNS). At present, this gene is only known from bacteria and fungi, as is the capacity to produce ß-lactam antibiotics. The FcIPNS gene was located on two genomic contigs, was physically linked to a predicted insect ATP-binding cassette transporter gene, and contained three introns each flanked by eukaryotic splicing recognition sites (GT/AG). Homology searches revealed no similarity between these introns and the FcIPNS regions of bacteria or fungi. All amino acids conserved across bacteria and fungi were also conserved in F. candida. Recombinant FcIPNS was able to convert its substrate amino δ-(l-α-aminoadipyl)-l-cysteinyl-d-valine into isopenicillin N, providing strong evidence that FcIPNS is functional. Phylogenetic analysis clustered FcIPNS outside the bacterial IPNS clade, and also outside the fungal IPNS clade, suggesting an ancient gene transfer followed by divergence in the F. candida genome. In conclusion, the data suggest that the soil-living collembolan F. candida has assimilated the capacity for antibacterial activity by horizontal gene transfer, which may be an important adaptive trait in the microbe-dominated soil ecosystem.

Gapless assembly of complete human and plant chromosomes using only nanopore sequencing.

The combination of ultra-long Oxford Nanopore (ONT) sequencing reads with long, accurate PacBio HiFi reads has enabled the completion of a human genome and spurred similar efforts to complete the genomes of many other species. However, this approach for complete, "telomere-to-telomere" genome assembly relies on multiple sequencing platforms, limiting its accessibility. ONT "Duplex" sequencing reads, where both strands of the DNA are read to improve quality, promise high per-base accuracy. To evaluate this new data type, we generated ONT Duplex data for three widely-studied genomes: human HG002, Solanum lycopersicum Heinz 1706 (tomato), and Zea mays B73 (maize). For the diploid, heterozygous HG002 genome, we also used "Pore-C" chromatin contact mapping to completely phase the haplotypes. We found the accuracy of Duplex data to be similar to HiFi sequencing, but with read lengths tens of kilobases longer, and the Pore-C data to be compatible with existing diploid assembly algorithms. This combination of read length and accuracy enables the construction of a high-quality initial assembly, which can then be further resolved using the ultra-long reads, and finally phased into chromosome-scale haplotypes with Pore-C. The resulting assemblies have a base accuracy exceeding 99.999% (Q50) and near-perfect continuity, with most chromosomes assembled as single contigs. We conclude that ONT sequencing is a viable alternative to HiFi sequencing for de novo genome assembly, and has the potential to provide a single-instrument solution for the reconstruction of complete genomes.

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.

Transcriptional responses indicate attenuated oxidative stress in the springtail Folsomia candida exposed to mixtures of cadmium and phenanthrene.

Since the 'omics revolution', the assessment of toxic chemical mixtures has incorporated approaches where phenotypic endpoints are connected to a mechanistic understanding of toxicity. In this study we determined the effect of binary mixtures of cadmium and phenanthrene on the reproduction of Folsomia candida and investigated the cellular mechanisms underlying this response. Mixture toxicity modeling showed an antagonistic deviation from concentration addition for reproduction effects of the mixtures. Subsequent transcriptional response analysis was done using five mixtures at the modeled 50 % effect level for reproduction. The transcription profiles of 86 high throughput RT-qPCR assays were studied by means of partial least squares regression analysis. The first and second principal components (PCs) were correlated with global responses to cadmium and phenanthrene, while correlations with the mixture treatments were found in the higher PCs. Specifically associated with the mixture treatments were a biotransformation phase II gene, four mitochondrial related genes and a gene involved in the biosynthesis of antioxidant selenoproteins. Membrane integrity related gene inductions were correlated with the single phenanthrene treatment but not with the mixtures. Immune and inflammatory response assays did not correlate with any of the mixtures. These results suggest moderated oxidative stress, a higher mitochondrial maintenance and less compromised membrane function in the mixture exposed samples compared to the separate cadmium or phenanthrene exposures. The antagonism found for inhibition of reproduction may partially originate from these differences. Mechanistic studies on mixture toxicity can ultimately aid risk assessment by defining relevant toxicity pathways in organisms exposed to real-world mixture exposures present in the field.

Validation of biomarkers for neonicotinoid exposure in Folsomia candida under mutual exposure to diethyl maleate.

Neonicotinoid insecticides are harmful to non-target soil invertebrates, which are crucial for sustainable agriculture. Gene expression biomarkers could provide economic and high-throughput metrics of neonicotinoid exposure and toxicity to non-target invertebrates. Thereby, biomarkers can help guide remediation efforts or policy enforcement. Gene expression of Glutathione S-Transferase 3 (GST3) has previously been proposed as a biomarker for the neonicotinoid imidacloprid in the soil ecotoxicological model species Folsomia candida (Collembola). However, it remains unclear how reliably gene expression of neonicotinoid biomarkers, such as GST3, can indicate the exposure to the broader neonicotinoid family under putative GST enzymatic inhibition. In this work, we exposed springtails to two neonicotinoids, thiacloprid and imidacloprid, alongside diethyl maleate (DEM), a known GST metabolic inhibitor that imposes oxidative stress. First, we determined the influence of DEM on neonicotinoid toxicity to springtail fecundity. Second, we surveyed the gene expression of four biomarkers, including GST3, under mutual exposure to neonicotinoids and DEM. We observed no effect of DEM on springtail fecundity. Moreover, the expression of GST3 was only influenced by DEM under mutual exposure with thiacloprid but not with imidacloprid. The results indicate that GST3 is not a robust indicator of neonicotinoid exposure and that probable GST enzymatic inhibition mediates the toxicity of imidacloprid and thiacloprid differentially. Future research should investigate biomarker reliability under shifting metabolic conditions such as provided by DEM exposure.

Nonribosomal Peptide Synthetases in Animals.

Nonribosomal peptide synthetases (NRPSs) are a class of cytosolic enzymes that synthesize a range of bio-active secondary metabolites including antibiotics and siderophores. They are widespread among both prokaryotes and eukaryotes but are considered rare among animals. Recently, several novel NRPS genes have been described in nematodes, schistosomes, and arthropods, which led us to investigate how prevalent NRPS genes are in the animal kingdom. We screened 1059 sequenced animal genomes and showed that NRPSs were present in 7 out of the 19 phyla analyzed. A phylogenetic analysis showed that the identified NRPSs form clades distinct from other adenylate-forming enzymes that contain similar domains such as fatty acid synthases. NRPSs show a remarkably scattered distribution over the animal kingdom. They are especially abundant in rotifers and nematodes. In rotifers, we found a large variety of domain architectures and predicted substrates. In the nematode Plectus sambesii, we identified the beta-lactam biosynthesis genes L-δ-(α-aminoadipoyl)-L-cysteinyl-D-valine synthetase, isopenicillin N synthase, and deacetoxycephalosporin C synthase that catalyze the formation of beta-lactam antibiotics in fungi and bacteria. These genes are also present in several species of Collembola, but not in other hexapods analyzed so far. In conclusion, our survey showed that NRPS genes are more abundant and widespread in animals than previously known.

Maternal exposure to polystyrene nanoparticles retarded fetal growth and triggered metabolic disorders of placenta and fetus in mice.

Microplastics can enter the human body via direct body contact or the food chain, increasing the likelihood of adverse impacts on pregnancy and fetal development. We investigated the potential effects and modes of action of polystyrene nanoplastics (PS-NPs) in placenta and fetus using mice as a model species. Maternal PS-NP exposure (100 nm; 1 and 10 mg/L) via drinking water induced a significant decline in fetal weights at the higher exposure concentration. Abnormal morphologies of cells in the placenta and fetus were observed after exposure. For the placenta, transcriptomic analyses indicated that PS-NPs significantly disturbed cholesterol metabolism and complement and coagulation cascades pathways. Metabolomics showed appreciable metabolic disorders, particularly affecting sucrose and daidzein concentrations. For the fetal skeletal muscle, transcriptomics identified many significantly regulated genes, involving muscle tissue development, lipid metabolism, and skin formation. Transcriptomic analysis of the placenta and fetal skeletal muscle at the high PS-NP concentration showed that APOA4 and its transcriptional factors, facilitating cholesterol transportation, were significantly regulated in both tissues. Our study revealed that PS-NPs caused fetal growth restriction and significantly disturbed cholesterol metabolism in both placenta and fetus, offering new insights into the mechanisms underlying the placental and fetal effects in mice exposed to PS-NPs.

Combining time-resolved transcriptomics and proteomics data for Adverse Outcome Pathway refinement in ecotoxicology.

Conventional Environmental Risk Assessment (ERA) of pesticide pollution is based on soil concentrations and apical endpoints, such as the reproduction of test organisms, but has traditionally disregarded information along the organismal response cascade leading to an adverse outcome. The Adverse Outcome Pathway (AOP) framework includes response information at any level of biological organization, providing opportunities to use intermediate responses as a predictive read-out for adverse outcomes instead. Transcriptomic and proteomic data can provide thousands of data points on the response to toxic exposure. Combining multiple omics data types is necessary for a comprehensive overview of the response cascade and, therefore, AOP development. However, it is unclear if transcript and protein responses are synchronized in time or time lagged. To understand if analysis of multi-omics data obtained at the same timepoint reveal one synchronized response cascade, we studied time-resolved shifts in gene transcript and protein abundance in the springtail Folsomia candida, a soil ecotoxicological model, after exposure to the neonicotinoid insecticide imidacloprid. We analyzed transcriptome and proteome data every 12 h up to 72 h after onset of exposure. The most pronounced shift in both transcript and protein abundances was observed after 48 h exposure. Moreover, cross-correlation analyses indicate that most genes displayed the highest correlation between transcript and protein abundances without a time-lag. This demonstrates that a combined analysis of transcriptomic and proteomic data from the same time-point can be used for AOP improvement. This data will promote the development of biomarkers for the presence of neonicotinoid insecticides or chemicals with a similar mechanism of action in soils.

Daphnia magna Multigeneration Exposure to Carbendazim: Gene Transcription Responses.

The world population is experiencing colossal growth and thus demand for food, leading to an increase in the use of pesticides. Persistent pesticide contamination, such as carbendazim, remains a pressing environmental concern, with potentially long-term impacts on aquatic ecosystems. In the present study, Daphnia magna was exposed to carbendazim (5 µg L-1) for 12 generations, with the aim of assessing gene transcription alterations induced by carbendazim (using a D. magna custom microarray). The results showed that carbendazim caused changes in genes involved in the response to stress, DNA replication/repair, neurotransmission, ATP production, and lipid and carbohydrate metabolism at concentrations already found in the environment. These outcomes support the results of previous studies, in which carbendazim induced genotoxic effects and reproduction impairment (increasing the number of aborted eggs with the decreasing number of neonates produced). The exposure of daphnids to carbendazim did not cause a stable change in gene transcription between generations, with more genes being differentially expressed in the F0 generation than in the F12 generation. This could show some possible daphnid acclimation after 12 generations and is aligned with previous multigenerational studies where few ecotoxicological effects at the individual and populational levels and other subcellular level effects (e.g., biochemical biomarkers) were found.

The influence of long-term copper contaminated agricultural soil at different pH levels on microbial communities and springtail transcriptional regulation.

Copper has long been applied for agricultural practises. Like other metals, copper is highly persistent in the environment and biologically active long after its use has ceased. Here we present a unique study on the long-term effects (27 years) of copper and pH on soil microbial communities and on the springtail Folsomia candida an important representative of the soil macrofauna, in an experiment with a full factorial, random block design. Bacterial communities were mostly affected by pH. These effects were prominent in Acidobacteria, while Actinobacteria and Gammaroteobacteria communities were affected by original and bioavailable copper. Reproduction and survival of the collembolan F. candida was not affected by the studied copper concentrations. However, the transcriptomic responses to copper reflected a mechanism of copper transport and detoxification, while pH exerted effects on nucleotide and protein metabolism and (acute) inflammatory response. We conclude that microbial community structure reflected the history of copper contamination, while gene expression analysis of F. candida is associated with the current level of bioavailable copper. The study is a first step in the development of a molecular strategy aiming at a more comprehensive assessment of various aspects of soil quality and ecotoxicology.

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

Natural toxins, such as isothiocyanate (ITC), are harmful secondary metabolites produced by plants. Many natural toxins occur in commercial crops, yet their possible negative repercussions on especially non-target soil organisms are largely unknown. This study examined life-history and gene transcriptional responses to 2-phenylethyl ITC on two soil arthropod species: Folsomia candida and Protaphorura fimata. To that end the standardized ISO guideline for ecotoxicological tests and a microarray for F. candida were used. The dissipation of 2-phenylethyl ITC in natural soil was investigated using GC-MS/MS for quantification. Half-lives, tested at four concentration levels in natural soil, were on average 16 h with biodegradation as the plausible main removal process. Regardless, toxic effects on reproduction were shown for F. candida and P. fimata, with EC50 values of around 11.5 nmol/g soil illustrating the toxic character of this compound. Gene expression profiles revealed the importance of fatty acid metabolism at low exposure concentrations (EC10), which is associated with the lipophilic nature of 2-phenylethyl ITC. At higher concentrations (EC50) gene expression became more ubiquitous with over-expression of especially stress-related genes and sugar metabolism. The regulation of a gene encoding a precursor of follistatin, furthermore, implied the inhibition of reproduction and may be an important molecular target that can be linked to the observed adverse effect of life-history traits.